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The  objective  of  the  present  study  is  the  systematic  experi- 
mental investigation  of  interfacial  and  bed  shear  stresses  together  with 
the  shear  flow  structure  of  an  arrested  saline  wedge.  A special 
experimental  apparatus  was  used.  This  apparatus  was  designed  at  the 
Aristoteles  University  of  Thessaloniki  in  Greece  and  at  the  University 
of  Florida  in  Gainesville  and  was  constructed  at  the  Coastal  and  Ocean- 
ographic Engineering  Laboratory  of  the  latter  university.  The  apparatus 
consists  of  an  open  tilting  flume  45.7  cm  wide,  45.7  cm  deep , and  20. 24  m 
long,  supplemented  with  all  the  pertinent  units  for  the  formation  of 
arrested  saline  wedges  over  a large  range  of  flow  and  salinity  conditions. 
The  apparatus  was  equipped  with  a constant  temperature  hot  film  anemometer, 
a high  resolution  electronic  conductivity  meter,  and  some  other  instru- 
ments for  operation  and  measurements. 
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A number  of  74  experiments  have  been  carried  out.  The  analysis 

of  the  experimental  results  led  to  the  following  conclusions:  (1)  The 

best  correlation  of  the  interfacial  and  bed  friction  coefficients  with 

the  flow  parameters  and  fluid  properties  has  been  obtained  when  these 

2 

coefficients  were  plotted  against  the  dimensionless  parameter  ReFr  , 
where  Re  is  the  Reynolds  number  and  Fr  is  the  nondensimetric  Froude 
number,  with  the  relative  density  difference  Ap/(p  + Ap)  as  an  inde- 
pendent parameter.  (2)  The  ratio  of  interfacial  and  bed  shear  stresses 
is  independent  of  Ap/(p  + Ap).  (3)  Within  a core  of  about  4 mm  thickness, 
in  the  interfacial  zone,  the  nondimensional  velocity  distribution  is 
linear.  (4)  Outside  the  linear  zone  the  velocity  distribution  law  for 
the  two  layers  is  logarithmic.  In  these  zones  we  distinguish  the  law 
of  interface  and  the  law  of  the  outer  region,  corresponding  to  wall  law 
and  outer  region  law  for  fixed  boundaries,  respectively. 
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CHAPTER  I 
INTRODUCTION 

Density  stratification  occurs  in  the  atmosphere,  oceans,  lakes, 
estuaries,  rivers,  and  reservoirs  as  a result  of  compressibility,  heat, 
dissolved  salts,  and  suspended  matter.  It  can  be  natural  or  induced  by 
man-made  activities,  such  as  discharge  of  heated  water  and  other 
effluents.  The  induced  density  differences  may  strongly  affect  the 
dynamics  of  the  system  as  well  as  its  dispersive  and  transport  proper- 
ties; therefore,  they  play  a dominant  role  in  the  environmental  quality 
control . 

Density  stratification  can  be  either  continuous  or  discontinuous. 
In  the  latter  case  the  system  is  composed  of  a finite  number  of  distinct 
layers  of  constant  density  separated  by  relatively  thin  zones  of  sharp 
density  gradients.  For  modeling  purposes  these  zones  are  normally  simu- 
lated by  surfaces  through  which  densities  and  sometimes  velocities  are 
assumed  to  vary  discontinuously.  The  simplest  two-layered  system 
has  been  extensively  studied  for  two  reasons:  first,  a large  number  of 

natural  systems  and  problems  in  the  domain  of  hydraulic,  coastal,  and 
environmental  engineering  can  be  simulated  and  treated  as  two-layered 
systems;  second,  its  study  can  often  lead  to  an  understanding  and  mathe- 
matical formulation  of  more  complicated  cases.  Some  of  the  most  common 
occurrences  of  stratified  flows  are  the  warm  and  salinity  wedges,  the 
formation  and  evolution  of  the  thermocline  in  oceans  and  lakes,  the  lock 
exchange  flow,  the  density  and  turbidity  currents  in  estuaries,  oceans, 
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and  lakes  and  the  motion  of  fluid  fronts.  In  this  study  the  term 
"stratified  flow"  refers  to  the  two-layered  system  and  the  interface 
is  considered  as  a thin  zone  of  finite  width. 

In  the  modeling  of  stratified  flow  systems  and  in  the  analysis 
of  the  associated  problems,  the  interfacial  shear  stresses,  x . , as  well 
as  the  associated  coefficient,  f . , has  to  be  introduced.  For  this 
reason,  x.  and  f^  have  been  studied  extensively  during  the  last  thirty 
years  through  theoretical  analysis,  laboratory  research,  and  field 
investigations.  Some  of  these  studies  aimed  at  a mere  correlation  of 
f.  with  the  pertinent  gross  flow  parameters  and  fluid  properties,  while 
some  others  concentrated  in  the  detailed  flow  structure  on  the  neighbor- 
hood and  within  the  interfacial  zone.  The  first  group  of  studies 
resulted  in  a number  of  semi -empirical  laws  for  f . . The  application 
of  these  laws  to  laboratory  and  field  data  obtained  by  other  investiga- 
tions leads,  sometimes,  to  values  of  f.  which  may  differ  by  orders  of 
magnitude.  The  reason  for  these  discrepancies  lies  in  the  basic  assump- 
tions, the  methodology  of  approach,  and  the  range  of  variation  of  the 
controlling  parameters.  The  reliability  of  the  up  to  1975  semi- 
empirical  laws  for  f.  and  the  derivation  of  a relationship  of  general 
validity  linking  the  interfacial  friction  coefficient,  f . , to  the  per- 
tinent gross  flow  and  density  parameters  were  the  main  objectives  of  the 
author's  doctoral  dissertation  [11],*  which  for  the  first  time  led  to 
satisfactory  generalized  correlations  for  f . . 

Stratified  flows  can  be  divided  into  the  following  three  general 

types: 


*Numbers  in  brackets  indicate  references  at  the  end  of  the  dissertation. 
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a.  Flow  of  lighter  fluid  over  a quasi-stagnant  denser 
fluid,  as  in  the  case  of  arrested  saline  wedges  in 
stratified  estuaries 

b.  Flow  of  heavier  fluid  under  a quasi-stagnant  lighter 
fluid,  as  in  arrested  thermal  wedges  in  rivers  and 
estuaries,  in  density  and  turbidity  currents  in  oceans, 
lakes,  and  reservoirs,  and  in  advancing  density  fronts 

c.  Motion  of  both  layers  in  opposite  directions,  as  in  the 
case  of  lock  exchange  flow  (counter  flow) 

The  hydrodynamic  behavior  of  arrested  saline  wedges  has  been 
theoretically  and  experimentally  studied  in  the  last  thirty  years.  Among 
the  best  known  studies  are  Schijf  and  Schonfeld's  theoretical  analysis 
[37]  and  Keulegan's  theoretical  and  experimental  research  [22,  24].  The 
fundamental  difference  between  these  two  approaches  lies  in  the  introduc- 
tion of  bed  shear  stresses.  Schijf  and  Schonfeld  [37]  defined  the  bed 
shear  stresses  as  a function  of  the  mean  velocity  of  the  lower  layer. 

This  velocity  goes  to  zero  in  each  cross  section  of  the  wedge;  therefore, 
the  bed  shear  stresses  by  definition  are  zero.  Keulegan  [22,  24],  on 
the  other  hand,  through  an  integration  of  the  equations  of  motion  in  each 
layer  and  using  the  actual  velocity  profiles,  derived  expressions  for  the 
bed  and  interfacial  shear  stresses  Tq  and  x. , respectively,  and  the 
associated  friction  coefficients  fg  and  f .. . His  experiments  showed  that 
the  bed  shear  stresses  are  substantial.  The  bed  shear  stresses  distribu- 
tion is  important  since  these  stresses  control  the  overall  flow  dynamics 
and  specifically  the  sediment  erosion  and  deposition  [28,  29].  Therefore, 
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a reliable  estimate  of  this  distribution  would  contribute  to  an  improved 
design  of  navigation  canals  and  to  a better  environmental  control. 

In  the  analysis  of  any  problem  which  is  associated  with  mass, 
energy,  and  momentum  transport,  a knowledge  of  the  vertical  velocity 
distribution  is  of  primary  importance.  From  the  analysis  of  the  first 
experimental  results,  it  became  obvious  that  the  nondimensional  velocity 
distribution  in  the  vicinity  of  the  interface  follows  a logarithmic  law 
as  in  the  case  of  fixed  boundaries.  These  first  results  were  reinforced 
by  Csanady's  [9]  and  Powell's  [34]  data.  In  general,  there  are  limited 
studies  on  the  laws  which  govern  the  nondimensional  velocity  distribu- 
tion in  strongly  stratified  flows.  Among  the  most  interesting  of  them  are 
the  studies  by  Csanady  [9],  Powell  [34],  and  Hino  et  al.  [17]. 

Csanady  [9]  elaborated  on  Lofquist's  experimental  results  [26] 
on  stratified  flows  where  there  was  uniform  motion  of  the  lower  denser 
layer  under  a quasi-stagnant  upper  lighter  layer.  These  experiments 
took  place  before  1960  without  any  use  of  modern  electronic  instrumenta- 
tion. Csanady  [9],  in  his  effort  to  give  a universal  law  for  the  non- 
dimensional  velocity  distribution,  defined  a constant  coefficient  < 
analogous  to  the  von  Karman  universal  constant  k.  Powell  [34],  in  his 
study  on  fresh  water  flow  over  a quasi-stagnant  salt  water  layer,  found 
that  the  velocity  distribution  law  in  the  upper  layer  is  logarithmic 
and  claimed  that  the  distribution  depends  on  an  equivalent  von  Karman 
constant.  His  data,  however,  were  not  sufficient  to  permit  the  deriva- 
tion of  the  pertinent  law.  Hino  et  al . [17]  studied  the  dynamics  of 
arrested  saline  wedges  and  they  considered  a 1 ogarthmic-1 inear  velocity 
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distribution  for  the  upper  layer  according  to  Monin-Obukov  similarity 
theory.  Keulegan  [22,  24]  and  some  other  investigators  gave  a dimen- 
sional or  nondimensional  velocity  distribution  without  making  an  effort 
to  find  the  governing  law. 

The  objectives  of  the  present  study  were  (1)  the  experimental 
investigation  of  interfacial  and  bed  shear  stresses  and  especially  the 
derivation  of  generalized  laws  linking  these  stresses  to  the  pertinent 
gross  flow  parameters  and  fluid  properties  and  (2)  the  study  of  the  flow 
structure  in  the  vicinity  of  the  interface  and  in  the  upper  and  lower 
layer,  for  the  purpose  of  obtaining  a velocity  distribution  law  of 
universal  validity. 

The  experiments  were  conducted  in  an  open  salinity  flume  con- 
structed in  the  Coastal  and  Oceanographic  Laboratory  of  the  University 
of  Florida,  in  Gainesville,  Florida,  U.S.A.  The  apparatus  consists  of  a 
20.24  m long  open  channel  of  a square  cross  section  45.7  cm  x 45.7  cm 
constructed  of  transparent  plexiglass,  and  of  the  accessory  tanks  for 
recirculation  of  fresh  water  and  the  supply  of  salt  water.  The  system 
was  equipped  with  the  necessary  electronic  instrumentation  for  the 
measurements  of  velocities,  salinities,  temperatures,  and  Reynolds 
stresses. 

The  investigation  herein  described  led  to  original  correlations 
for  the  interfacial  and  bed  shear  stress  coefficient  and  to  equally 
original  velocity  distribution  laws  in  the  neighborhood  and  within  the 
interfacial  zone.  Useful  and  interesting  results  have  also  been  obtained 
regarding  the  nature  and  the  validity  of  the  interfacial  friction  coef- 
ficient obtained  by  Schijf  and  Schonfeld's  simplified  model  [37]. 


CHAPTER  II 

BACKGROUND  AND  REVIEW  OF  EARLIER  INVESTIGATIONS 


2.1  Basic  Parameters 

The  simple  stratified  flow  system  of  Figure  2.1,  where  denser 
liquid  flows  over  a spillway,  is  first  analyzed  in  order  to  introduce 
some  basic  concepts  of  stratified  flows  [15],  The  two  layers  are  con- 
sidered frictionless  and  nonmixing.  Therefore,  the  velocity  and  the 
density  distribution  are  discontinuous  at  the  interface.  The  upper  layer 
is  stagnant  and  the  velocity  of  the  lower  layer  is  considered  uniform 
over  any  cross  section. 

Application  of  the  energy  equation  between  points  A and  B at  the 
interface,  first  for  the  stagnant  upper  layer  and  second  for  the  moving 
lower  layer,  gives 

V2  - V? 

B ? A = ^ gAh  (2.1.D 

2 P2 

where  g is  the  gravity  acceleration  with  the  other  parameters  identified 
in  Figure  2.1.  If  the  velocity  VA  is  very  small  compared  with  the 
velocity  Vg,  Equation  (2.1.1)  becomes 


In  the  case  where  Ap  = p2,  as  in  air-water  systems.  Equation  (2.1.2) 
becomes  identical  with  the  free  surface  flow  equation  over  a spillway. 


6 


7 


V,  = 0 

Pi  = P 


Figure  2.1  Gravity  flow  in  a two- layered  system 
Source:  [15,  Figure  26.1] 
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In  the  case  where  Ap  = 0.0,  as  in  homogeneous  fluids,  the  velocity  is 
zero  (no  motion).  Therefore,  the  driving  force  is  the  submerged  weight 
of  the  denser  fluid,  which  corresponds  to  a reduced  gravitational 
acceleration  denoted  by 


g' 


Ap 

p + Ap 


g 


(2.1.3) 


and  defined  as  "densimetric  gravitational  acceleration."  For  small 
values  of  Ap  compared  with  p.  Equation  (2.1.3)  takes  the  form 


(2.1.4) 


Substituting  g by  g1  in  the  definition  of  Froude  number,  a new  nondi- 
mensional  number  is  introduced  as 


Fr'  = — (2.1.5) 

v^Th 

where  V and  h represent  characteristic  velocity  and  length  scales  of  the 
flow  system,  respectively.  This  number  is  called  "densimetric  Froude 
number"  and  it  is  a basic  parameter  in  stratified  flows  for  dynamic 
similarity. 

In  1949,  Keulegan  [21],  in  order  to  study  the  interfacial 
stability  of  viscous  fluids,  introduced  a nondimensional  parameter  0 
defined  by  the  equation 


(2.1.6) 
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where  V is  the  local  average  velocity  of  the  moving  layer  and  v is  the 
kinematic  viscosity  of  the  heavier  fluid.  This  equation  comes  from 
Jeffrey's  criterion  [20]  for  the  generation  of  wind  waves.  Equation 
(2.1.6)  can  also  be  written  in  the  form 

K3  = = 7 (2.1.7) 

Re(Fr') 

where  Re  is  the  Reynolds  number.  K is  known  as  "Keulegan's  number." 

In  the  system  of  Figure  2.1,  the  velocity  and  density  distribu- 
tions are  discontinuous  at  the  interface.  Viscosity  and  mixing,  how- 
ever, lead  to  a continuous  distribution  of  both  variables  around  the 
interface,  albeit  at  sharp  rates.  The  Richardson  number,  Ri , defined  by 
the  relation 

Ri  = - 2.JM4Z)  (2.1.8) 

p (du/dy)2 

is  introduced  as  a stability  parameter  of  flow  where  dp/dy  and  du/dy  are 
the  density  and  velocity  gradients,  respectively.  When  the  flow  is  con- 
sidered as  one-dimensional,  the  Richardson  number  becomes  identical 
with  the  densimetric  Froude  number  [Equation  (2.1.5)]. 

2.2  Mixing  Mechanisms  and  Density  Interface 
The  main  mixing  mechanisms  are  (a)  the  entrainment  and  (b)  the 
turbulent  diffusion.  Entrainment  is  the  vertical  one-way  volume  flux 
through  the  interface  from  the  lower  to  the  upper  layer.  Turbulent 
diffusion  is  the  mass  transport  of  the  dissolved  or  suspended  matter  by 
the  turbulent  velocity  fluctuations.  As  a result  of  the  mixing,  a thin 
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layer,  called  the  interfacial  layer  or  interfacial  zone,  is  formed  with 
high  gradients  of  continuously  varying  density.  In  the  present  study 
the  term  "interface"  indicates  the  surface  at  which  the  density  is  the 
average  of  the  two  fluid  densities. 

2.3  Velocity  Interface 

Experimental  observations  suggest  that  in  the  vicinity  of  the 
interface  there  is  a zone  where  the  velocity  gradients  are  large. 

Linden  [25]  calls  this  zone  "velocity  interface  zone."  Lofquist  [26] 
observed  that  the  point  of  inflection  of  the  velocity  distribution 
curve  is  about  1 mm  to  3 mm  over  the  interface.  In  Powell's  experiments 
[34],  the  above  point  of  inflection  was  found  to  be  at  a pus f.tfon where  the 
velocity  was  the  average  of  the  velocities  U-^  and  of  the  two 
layers  very  far  from  the  velocity  interfacial  zone.  In  this  study  the 
term  "velocity  interface"  designates  a theoretical  surface,  where  the 
velocity  has  the  mean  value  of  and  l^.  A characteristic  of  the 
velocity  interface  is  that  the  viscous  stresses  in  its  vicinity  are 
maximized. 

From  Figure  2.2  we  can  see  the  influence  of  interfacial  zone 
and  velocity  interfacial  zone  widths  on  Richardson  number,  Ri . The  two 
interfaces  usually  do  not  coincide.  The  stratified  flow,  which  is 
studied  herein,  belongs  to  the  second  case  of  Figure  2.2. 

2.4  Review  and. Discussion  of  Earlier  Investigations 
2,4.1  The  Interfacial  Friction  Coefficient 

The  scientific  interest  and  practical  significance  of  the 
interfacial  friction  coefficient,  f.,  motivated  extensive  laboratory 
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(a)  Thick  density  interface 


(b)  Thin  density  interface 


Figure  2.2  The  effect  of  width  of  interfacial  zone  on  the  Richardson 
number 


Source:  [34] 
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research  and  field  investigations  in  the  last  thirty  years  in  order  to 
obtain  a reliable  correlation  of  f..  to  the  gross  flow  parameters  and 
fluid  properties,  which  define  the  flow. 

An  extensive  review  and  discussion  of  the  best  known  investiga- 
tions was  presented  in  Chapter  II  of  the  author's  previous  dissertation 
[11].  For  this  reason  only  a short  review  of  earlier  investigations 
will  be  given  in  this  section  with  some  more  details  on  Keulegan's 
experiments  [20,  24]  and  on  the  author's  previous  research  [11]. 

For  the  arrested  saline  wedge,  only  Keulegan's  experimental 
studies  [22,  24]  give  essentially  bottom  shear  stresses.  His  mathemati- 
cal model,  modified  and  refined,  is  used  in  this  study  for  the  formula- 
tion of  bottom  shear  stresses.  The  Schijf  and  Schonfeld  mathematical 
model  [37]  is  widely  known  for  the  evaluation  of  the  interfacial  fric- 
tion coefficient  f. . Finally,  in  the  author's  earlier  work  [11]  a 
satisfactory  correlation  of  f..  to  the  mean  flow  parameters  and  the  fluid 
properties  has  been  obtained.  This  correlation  can  be  used  as  a first 
approximation  of  f.  in  practical  applications.  Moreover,  in  the 
present  study,  f.  is  also  correlated  with  the  same  parameters  as  in 
[11,  30],  with  satisfactory  results. 

Keulegan  [22]  applied  the  equations  of  motion  in  each  layer  to 
compute  the  interfacial  and  bed  shear  stresses  of  arrested  saline 
wedges.  His  basic  assumptions  were  (a)  negligible  vertical  accelera- 
tions, (b)  hydrostatic  pressure  distribution,  and  (c)  small  density 
difference  between  the  two  layers  as  compared  to  the  average  density. 

Keulegan's  analysis  [22,  24]  led  to  the  following  equations  for 
the  mean  interfacial  and  bed  shear  stresses  T^  and  Tq,  respectively: 
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x=0 


Figure  2.3  Arrested  saline  wedge 
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(2. 4. 1.2) 


dH 

Here  is  the  mean  slope  of  the  piezometric  line, 


a = (1/A.,) 


(u/V.TdA, 


A-,  is  the  cross  sectional  area  of  the  upper  layer,  u is  the  local  time- 
averaged  velocity,  V-,  is  the  mean  velocity  of  the  upper  layer,  fwQ  is 
the  wall  friction  coefficient  upstream  of  the  wedge, 


is  the  densimetric  Froude  number  upstream  of  the  wedge,  Vq,  hg  are  the 
velocity  and  the  depth  upstream  of  the  wedge,  respectively,  Ap  is  the 
density  difference  between  the  two  layers,  p is  the  fresh  water  density, 
g is  the  gravitational  acceleration,  h-,c  is  the  depth  of  the  upper  layer 
at  the  ocean  end  of  the  channel,  hsc  = hg  - h^c,  L is  the  length  of  the 
arrested  saline  wedge,  b is  the  width  of  the  channel.  Now,  introducing 


n as 
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(2. 4. 1.3) 


I , 1^,  Ic,  and  1^  are  given  by  the  following  formulas: 


1=1-  0.485n 

a 


(2. 4. 1.4) 


1.028n 
1 - n 


(2. 4. 1.5) 


1=1+  0.369n  + 0. 195n 
c 


2 


(2. 4. 1.6) 


Id  = 0.496 


(2.4.1 .7) 


Keulegan  [22,  24]  computed  the  stresses  T.  and  Tg  using  Equations 
(2. 4. 1.1)  and  (2. 4. 1.2)  and  data  from  his  experimental  results.  Figures 


upstream  of  the  wedge.  The  scattering  of  the  points  in  Figures  2.4 
and  2.5  is  significant.  This  means  that  there  are  some  other  parameters 
besides  the  Reynolds  number,  Re,  which  influence  the  stresses.  In 


one  order  of  magnitude  for  the  same  value  of  Re.  Keulegan,  on  the  basis 
of  the  mean  values  of  the  experimental  points  of  Figures  2.4  and  2.5 
suggested  the  graphical  relationships  of  Figure  2.6.  This  figure  indi- 
cates that  the  stresses  T.  are  inversely  proportional  to  the  Reynolds 
number.  Re.  This  observation  led  him  to  consider  the  interfacial 
stresses  as  laminar.  This  is  contrary  to  the  experimental  results  of 
several  investigators,  where  significant  Reynolds  stresses  were  detected 
[3,  4,  12,  13,  19,  26,  35,  39]. 


2 

2.4  and  2.5  give  the  values  of  T./pVg  in  terms  of  the  Reynolds  number 


2 2 

particular,  the  extreme  values  of  T^/pVg  or  Tg/pVg  differ  by  more  than 


16 


Figure  2 
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Figure  2.6 
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It  is,  moreover,  observed  that  for  large  Reynolds  numbers,  the 
shear  stresses  at  the  bottom  are  larger  than  the  corresponding  inter- 
facial shear  stresses. 

In  1977,  the  author  completed  a study  on  the  dynamics  of  a 
two-layered  system  which  consisted  of  a fresh  water  layer  flowing  over 
a quasi-stagnant  salt  water  layer.  The  analysis  of  the  experimental 
data  was  based  on  a one-dimensional  steady  uniform  flow  model,  similar 
to  that  developed  by  Schijf  and  Schonfeld  [37]  for  open  channels.  Com- 
parison of  the  data  of  the  latter  study  with  data  obtained  by  other 
investigators  for  the  same  type  of  stratified  flow  is  given  in  Figure 
2.7  [11,  Figure  6.13].  The  scattering  of  the  data  points  does  not  per- 
mit any  reliable  conclusion  on  the  variation  law  for  f. . The  same 
scattering  appears  in  other  types  of  stratified  flows  when  f..  is  plotted 
versus  the  Renyolds  number  only.  Figure  2.8  [11,  Figure  6.14]  gives  an 
example  of  such  scattering  for  underflows.  From  the  comparison  of 
Figures  2.4,  2.7,  and  2.8,  it  can  be  concluded  that  f.  is  not  only  a 
function  of  Re,  but  it  also  depends  on  some  other  parameters.  A satis- 
factory correlation  was  obtained  when  f.  was  expressed  as  a function  of 
2 

ReFr  , where  Fr  is  the  nondensimetri c Froude  number  of  the  flow,  with 

Ap/p  as  an  independent  parameter  (Figure  2.9).  In  Figure  2.10  the  line 

of  Figure  2.9  for  Ap/p  = 0.020  has  been  reproduced,  while  reanalyzed 

data  by  Abraham  and  Eysink  [1]  for  counterflows  have  been  plotted.  The 

latter  fall  quite  closely  on  the  present  curve  and  its  extrapolation. 

Finally,  one  point  corresponding  to  a single  set  of  data  from  the  South 

Pass  of  the  Mississippi  River  has  been  added.  In  spite  of  the  wide  gap 

2 

of  the  parameter  ReFr  , that  point  falls  on  a smooth  extrapolation  of 
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Figure  2.7  Correlation  of  fn-  with  Reynolds  number  for  upper  layer  motion  and  semi-stagnant  lower 
layer 

Source:  [11,  Figure  2.7] 
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Figure  2.8  Correlation  of  f-j  with  Reynolds  number  for  underflows 
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the  Ap/p  = 0.020  line  of  the  author's  study  [11].  This  is  indicative 
that  the  correlation  parameters  have  been  properly  selected. 

In  Figure  2.11,  the  curves  of  Figure  2.9  have  been  extrapolated 
in  a qualitative  manner  following  roughly  the  extrapolation  of  the 
Ap/p  = 0.20  curve.  This  extrapolation  should  be  considered  tentative, 
since  only  one  field  point  was  used.  The  graphs  of  Figures  2.9  and  2.11 
can  be  used  as  a first  approximation  for  the  f.  value  for  any  type  of 
stratified  flow. 

In  1953  Schijf  and  Schonfeld  [37]  studied  a system  of  two  layers 
of  infinite  width  and  particularly  the  arrested  saline  wedge.  For  the 
latter  they  gave  the  following  equation  for  the  mean  value  of  f. : 


f . 


2h, 


5(Fr*)‘ 


2 + 3(Fr' )2/3 


— (Fr ' )4/3 
5^rr0; 


(2. 4. 1.8) 


In  1953  Bata  and  Knezevich  [8]  studied  the  turbidity  currents  in 
reservoirs  as  a density  stratified  underflow.  They  found  f = 0.010  for 
smooth  bed  and  f = 0.020  for  rough  bed,  where  f incorporates  the 
resistance  of  both  the  solid  boundaries  and  the  interface. 

In  1957  Bata  [7]  studied  the  warm  wedge  using  Schijf  and  Schon- 
feld's  [37]  equation. 

In  1960  Lofquist  [26]  studied  the  uniform  flow  of  a heavier 
salt  water  under  a quasi-stagnant  fresh  water  layer  using  the  equations 
of  motion.  Besides  the  detailed  velocity  and  density  profiles  at  the 
interfacial  zone,  he  introduced  the  ratio  of  the  total  stresses  to  the 
laminar  stresses  and,  neglecting  the  significant  scattering  of  the  data 
points,  he  derived  the  relation 
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(2. 4. 1.9) 


where  u is  the  maximum  cross  sectional  velocity  in  the  lower  layer, 
max 

v is  the  mean  kinematic  viscosity  of  the  layers,  (du/dy)  v is  the  maxi- 

Mia  X 

mum  velocity  gradient,  and  Re  is  the  Reynolds  number  of  the  lower  layer 

defined  as  Re  = ( ^2Rh2 ^ ^ ^ v2 ^ ‘ Here  V2’  v2’  Rh2  are  ttie  mean  veloci'ty> 
hydraulic  radius,  and  kinematic  viscosity  of  the  lower  layer,  respec- 
tively. Lofquist  [26]  subsequently  introduced  the  ratio  3mnv  as 
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(2.4.1.10) 


where  A is  the  interfacial  friction  coefficient  given  as  X = 
max  3 max 

— 2 — 

x /pu  , x is  the  maximum  shear  stress,  and  p is  the  mean  density 
max  K max  max 

of  the  fluids.  Finally,  he  found  that  3 is  a function  of  the 

fTlaX 

Reynolds  number.  Re,  and  of  the  densimetric  Froude  number,  Fr‘.  The 
scattering  of  the  data  points,  however,  was  substantial.  In  the  case 
of  laminar  flow  at  the  interface  3 becomes  equal  to  unity. 

In  1962  Macagno  and  Rouse  [27]  studied  the  case  of  counterflow 
in  a closed  rectangular  duct.  They  found  that  the  ratio  of  the  total 
stress  to  the  viscous  stress  at  the  interface  is  a function  of  the 
Reynolds  number  and  of  the  densimetric  Froude  number.  These  numbers 
were  based  on  the  maximum  algebraic  velocity  difference  and  on  the 
distance  between  the  points  of  maximum  velocities  in  the  moving  layers. 
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Their  analysis  was  based  on  the  assumption  of  laminar  flow  at  the  fixed 
boundaries. 

In  1964  Iwasaki  [18]  studied  stratified  flows  with  a fresh  water 

upper  layer  moving  over  a quasi -stagnant  salt  water  lower  layer.  He 

2 

suggested  a close  dependence  of  f.  on  Keulegan's  number  K = Re(Fr')  . 

A summary  of  data  from  several  sources,  including  Iwasaki 's  results,  is 

given  by  Sherenkov  et  al . [39,  Figure  6]  and  by  Plate  and  Friedrich  [32, 

Figure  4].  It  is  observed  that  Iwasaki 's  data  points  [18]  fall  on  a 

straight  line  with  little  scattering;  however,  all  other  points  fall  on 

a distinctly  different  line  whose  mean  slope  is  considerably  smaller 

than  that  given  by  Plate  and  Friedrich  [32]. 

Averkiev  and  Kind  [4]  investigated  the  dynamics  of  stratified 

flows  in  cooling  ponds.  They  defined  f.  as  a function  of  the  Reynolds 

2 

stresses  only  by  the  relation  f.  = -8u  v/V  , where  V is  the  mean  cross 
sectional  velocity  and  where  u ' v 1 was  measured  directly.  They  found 
that  f.j  is  a function  of  Re  for  small  values  of  the  Re  while  for  large 
Re(Re  = Vh/v  > 10,000)  f.  was  found  to  be  constant  and  equal  to  about 
0.016. 

Shi-Igai  [40]  analyzed  the  flow  of  an  upper  fresh  water  layer 
over  a quasi-stagnant  salt  water  layer.  He  defined  the  interfacial 
shear  stresses  by  the  relation  t..  = e/V-j , where  e is  the  viscous  energy 
dissipation  per  unit  surface  of  fluid  caused  by  interfacial  harmonic 
waves  and  V-j  is  the  uniform  velocity  of  the  moving  layer.  He  derived 
the  following  equation: 
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= 32tt 
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Re(Fr')' 


(2.4.1.11) 


where  Re  = V-jh^/v,  Fr'  = J g h-j  , h-j  is  the  thickness  of  the  upper 

layer  and  a and  X are  the  amplitude  and  the  length  of  the  interfacial 

2 2 

waves,  respectively.  The  factor  A-j  = 32rr  (a/X)  was  found  to  be  about 
48.  This  value  was  based  on  his  limited  experimental  data. 

In  1966  Shi-Igai  and  Kono  [41]  studied  the  arrested  saline  wedge 
and  they  expressed  f.  by  the  relation 


A 


Re(Fr ' )‘ 


(2.4.1.12) 


where  Re  and  Fr'  are  the  Reynolds  and  densimetric  Froude  number  upstream 
of  the  wedge,  respectively,  and  A is  an  experimental  constant  equal  to 
about  80. 

In  1971  Abraham  and  Eysink  [1]  investigated  the  case  of  exchange 

flows  with  constant  densimetric  Froude  number,  Fr'  = 0.8,  on  the  basis 

of  Schijf  and  Schonfeld's  equation  [37].  Moreover,  they  reanalyzed 

Keulegan's  [23]  and  Barr's  [5,  6]  experimental  data.  They  came  to  the 

conclusion  that  f ..  is  only  a function  of  Re  for  small  values  of  Re  and 

that  f.  is  constant  for  large  values  of  Re. 

In  1971,  Polk  et  al . [33]  studied  the  dynamics  of  warm  wedges  in 

the  field  and  computed  f.  by  equations  similar  to  those  derived  by  Schijf 

and  Schonfeld  [37].  The  coefficient  f.  ranged  from  0.08  to  0.013  for 
5 7 

8x10  ^ Re  ^ 2.3x10  and  for  various  densimetric  Froude  numbers. 
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In  1971  Riddell  [36]  conducted  an  extensive  investigation  on  the 
arrested  saline  wedges.  His  experimental  results  were  reanalyzed  at 
the  Delft  Hydraulic  Laboratory  using  Schijf  and  Schonfeld's  mathematical 
model  [37]  to  obtain  values  of  f. . 

In  1972  Georgiev  [13]  studied  uniform  underflows  using  photo- 
measurements to  determine  instantaneous  velocities.  A number  of  70  to 
90  instantaneous  velocities  were  used  to  calculate  the  time-average 
velocity  and  Reynolds  stresses  in  each  point  of  the  cross  section.  The 
total  interfacial  shear  stress  was  given  by  the  relation 

x.  = - puV  + y ^ (2.4.1.13) 


where  u'  and  w1  are  the  turbulent  velocity  fluctuations  in  the  direction 
of  flow,  x,  and  the  direction  of  gravity,  z,  respectively,  while  bars 

_ 2 

indicate  average  values.  The  coefficient  f.  was  defined  as  f.  = 8r../pV 

where  V is  the  mean  velocity  of  the  moving  layer.  Georgiev' s [13]  plot 

of  experimental  data  points,  based  on  the  hydraulic  radius  of  the  lower 

layer,  indicate  that  f.  is  a function  of  the  Reynolds  number,  Re,  only. 

3 4 

However,  in  his  experiments  that  number  ranged  from  10  to  10  only. 

In  1972  Harleman  and  Stolzenbach  [16]  in  their  study  of  the  re- 
circulation in  cooling  pools,  where  heated  water  is  discharged  from  power 
plants,  reviewed  data  of  previous  investigators.  For  the  turbulent  zone 
of  their  graph  and  with  significant  scattering  of  data  points  they  sug- 
gested the  following  functional  relation: 


1 


2.0  loaRe/fT  - 0.8 


(2.4.1.14) 
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In  1974  the  Delft  Hydraulics  Laboratory  [10]  selected  and 
reanalyzed  the  best  known,  at  that  time,  data  on  stratified  flows.  This 
summary  was  extensively  discussed  in  the  author's  earlier  work  [11]  and 
some  conclusions  are  herein  presented. 

a.  There  is  no  indication  of  any  consistent  relationship 
which  correlates  f.  either  to  the  Reynolds  number,  Re, 
or  to  the  densimetric  Froude  number,  Fr1. 

b.  In  turbulent  flows,  f.  decreases  as  Re  is  increased  for 
low  values  of  the  latter. 

c.  For  large  values  of  Re,  there  is  no  clear  relationship 
between  f^  and  Re. 

Figures  2.12  and  2.13  provide  additional  bases  for  the  previous 
conclusions  as  seen  from  the  scatter  of  data. 

In  1979  Abraham,  Karelse,  and  van  Os  [2]  assumed  that  the  inter- 
facial friction  coefficient,  f. , is  a function  of  the  Reynolds  number, 

Re,  the  Richardson  number,  Ri , and  that  it  depends  also  on  the  type  of 
flow.  Moreover,  they  stated  that  the  interfacial  shear  stresses,  t. , are 
affected  by  turbulence  which  is  generated  at  the  interface  and  at  the  bed. 
Particularly,  for  arrested  saline  wedges,  they  considered  that  t.  depends 
only  on  the  interfacial  turbulence.  This  is  not  the  real  case,  however, 
since  substantial  shear  flow  takes  place  near  the  bed,  which  affects  the 
velocity  profile  within  the  interfacial  zone  and  through  it  the  inter- 
facial shear  stresses.  They  claimed  that,  for  Reynolds  number  higher 

g 

than  10  , f.  is  constant  and  equal  to  0.0032. 

Among  the  latest  advances  in  stratified  flow  is  Powell's  [34] 
study  of  flow  of  an  upper  fresh  water  layer  over  a quasi-stagnant  salt 


Figure 


Figure 


2,12  Variation  of  the  average  f.  for  underflows  according  to 
Lofquist's  experimental  data 

Source:  [10,  Figure  2.6] 
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2.13  Variation  of  the  Average  f-j  for  underflows  according  to 
Lofquist's  experimental  data 


Source:  [10,  Figure  2.7] 
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water  layer  of  appreciable  depth.  Powell  [34]  calculated  the  inter- 
facial shear  stresses,  x. , by  direct  measurements  of  viscous  and  Reynolds 
stresses.  On  the  basis  of  his  experimental  results  the  interfacial 
friction  coefficient,  f. , appeared  to  be  a function  of  the  Reynolds 
number,  Re,  only.  However,  his  data  were  based  on  a limited  range  of 
variation  of  the  relative  density  difference  and  of  some  other  flow 
parameters . 

In  their  study  of  saline  wedges  in  1980,  Hino  et  al . [17]  con- 
sidered a linear- logarithmic  velocity  distribution  for  the  upper  layer 
according  to  Monin-Obukov  similarity  theory.  They  defined  the  inter- 
facial friction  coefficient  by  the  equation 


= 2 


(1  +.B)V1 


(2.4.1.15) 


where  u*^  is  the  interfacial  friction  velocity,  k is  the  von  Karman 
constant,  6 is  the  Monin-Obukov  constant,  is  the  average  velocity  of 
the  upper  layer,  and  is  the  slope  of  the  logarithmic  plot,  i.e., 

u*  • 

s2  ■ -r  <'  + 

p 

They  also  gave  f.  as  a function  of  Re(Fr')  with  Fr'  as  an  independent 
parameter.. 

The  discrepancy  among  various  formulas  and  methods  in  estimating 
values  of  f..  can  be  demonstrated  with  an  application  to  the  South  Pass 
of  the  Mississippi  River  for  the  following  data  [31]:  discharge, 
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Q = 2,837  m'Vsec.;  width,  b = 500  m;  depth,  h = 15  nr,  fresh  water 

velocity,  V g=  0.38  m/sec.;  relative  density  difference,  Ap/p  = 0.02; 

-6  2 

kinematic  viscosity,  v = 10  m /sec.;  densimetric  Froude  number 
upstream  of  the  wedge,  FrQ  = 0.22;  and  measured  salinity  intrusion, 

L = 22.5  km.  It  is  noted  that  the  critical  velocity  for  initiation  of 
entrainment  is  of  the  order  of  10  percent  of  the  actual  velocity;  there- 
fore, entrainment  does  take  place  through  the  theoretical  interface  and 
the  flow  within  the  interfacial  zone  is  expected  to  be  turbulent.  The 
following  results  are  obtained: 

a.  From  Schijf  and  Schonfeld's  [37]  model,  f^  = 0.00185 

b.  From  Keulegan's  [22,  24]  equation  for  laminar  flow  around 
the  interface,  f.  = 2.5xl0~^ 

c.  From  Bata  and  Knezevich's  [8]  equation  for  Re  = 

244.4xl05,  f.  = 0.74xl0"5 

d.  From  Shi-Igai's  [40,  41]  equation,  after  properly  con- 
verting Re  and  Fr1  to  mean  flow  conditions  over  the 

_5 

wedge,  f.  = 0.96x10  , and 

e.  From  Harleman  and  Stol zenbach ' s [16]  equation, 

f.  = 0.0071 

This  limited  review  led  to  the  following  conclusions:  (1)  The 

actual  value  of  f . , based  on  field  measurements,  may  exceed  the  values 
given  by  equations  based  on  laminar  flow  conditions  around  the  inter- 
face by  two  to  three  orders  of  magnitude,  (2)  f.  cannot  be  expressed  as 
a function  of  the  Reynolds  number.  Re,  only,  and  (3)  the  functional 
dependence  of  f.  on  the  Keulegan  number  seems  to  be  rather  weak. 
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These  conclusions  motivated  the  author  to  investigate  the 
interfacial  friction  coefficient,  for  the  first  time,  in  his  earlier 
work  [11].  His  analysis  was  based  on  a one-dimensional  model  similar 
to  that  by  Schijf  and  Schonfeld  [37]  for  open  channels.  The  main  results 
are  given  in  Figures  2.9,  2.10,  and  2.11.  The  graphs  in  the  figures  are 
the  result  of  many  trials  for  the  correlation  of  f.  to  the  gross  flow 
parameters  and  fluid  properties  and  can  be  used  as  a first  estimate  of  f-. 

2.4.2  Bed  Shear  Stresses  in  Arrested  Saline  Wedges 

The  only  experimental  work  on  the  bed  shear  stresses,  Tq,  was 
conducted  by  Keulegan  [22,  24]  (Figures  2.5,  2.6).  In  Figure  2.6,  the 
bed  shear  stresses,  Tq,  are  larger  than  the  interfacial  shear  stresses, 

T.j , at  least  by  one  order  of  magnitude  for  large  Reynolds  numbers. 

2.4.3  Velocity  Profiles  at  the  Interfacial  Zone 

Lofquist  [26]  was  one  of  the  first  investigators  who  studied 
systematically  the  interfacial  velocity  distribution.  He  gave  the  non- 
dimensional  velocity  distribution  by  the  functional  relationship 


u 


Ry  - yt)(au/3y)raax 


max 


(2.4.3.1 ) 


where  V. , y.  are  the  velocity  and  the  depth,  respectively,  at  the  inter- 
face, u*  is  the  friction  velocity  at  the  interface,  given  by  the  equation 
u*  = J Timax//p’  Timax  1S  maxinium  interfacial  shear  stress.  Equation 
(2. 4. 3.1)  is  given  graphically  by  Lofquist  [26]  in  Figure  2.14. 

In  1978  Csanady  [9]  reanalyzed  Lofquist's  experimental  data  [26] 
to  study  the  turbulence  structure  around  the  interface.  His  nondimensional 
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Figure  2.14  Velocity  profile 

Source:  [26,  Figure  14] 
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velocity  distribution  law,  defined  as  "law  of  the  interface,"  cor- 
responds to  the  "wall  law"  for  fixed  boundaries  and  is  given  by  the 
equation 


where  k is  the  von  Karman  constant,  vg  is  the  effective  viscosity, 
defined  as 


(du/dy)^  is  the  velocity  gradient  at  the  interface  and  1C  is  the 
Keulegan  number  at  the  interface  defined  as 


The  law  of  the  interface  raises  some  questions.  Specifically, 
the  inverse  of  the  coefficient  of  the  logarithmic  law  is  equal  to  von 
Karman's  constant  for  fixed  boundaries,  where  no  transport  of  mass  and 
momentum  takes  place  through  the  boundary.  However,  such  transport  does 
take  place  through  the  interface  of  the  two-layered  stratified  system. 

Powel  [34],  in  his  recent  and  aforementioned  doctoral  disserta- 
tion, studied  the  velocity  and  density  distribution  near  the  interface. 
The  nondimensional  velocity  distribution  is  given  in  Figure  2.15  [34, 
Figure  27].  For  a region  over  the  interface  the  graph  has  a logarithmic 


(2. 4. 3. 2) 


Ve  (du/dy) . 


(2. 4. 3. 3) 


(2. 4. 3. 4) 
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Figure  2,15  Normalized  mean  velocity  profile  from  Powell's  data 


Source:  [34,  Figure  27] 
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form.  Powell  [34]  postulated  that  the  transport  of  mass  and  momentum 
through  the  interface  may  modify  the  turbulence  structure  in  the  vicin- 
ity of  the  latter  and  he  gave  the  following  general  form  of  the 
velocity  distribution  law: 


FA> tn 


u*(y  - yfr 


V 


+F2<Ki> 


(2.4.3.5) 


where  F^(K. ) and  F2 ( ) are  functions  of  the  Keulegan  number  defined  at 
the  interface,  Powell's  data  [34]  were  not  sufficient  to  permit  the  deriva- 
tion of  a law  for  either  F^(K. ) or  F2(K-). 

Hino  et  al . [17]  also  made  use  in  their  aforementioned  recent 
work  of  the  Monin-Obukov  similarity  theory  and  claimed  that  velocity 
profiles  obey  the  "log-linear"  law  or  the  "bi-log"  law. 

In  the  case  of  arrested  saline  wedges  Keulegan  [22,  24]  gave  a 
nondimensional  velocity  distribution,  shown  in  Figure  2.16,  without  any 
information  on  the  flow  parameters  and  salinity. 

The  dimensional  or  non-dimensional  velocity  distribution  laws 
given  by  various  investigators  differ  not  only  for  different  types  of 
stratified  flows,  but  also  for  the  same  type  of  flow.  The  graphs  given 
by  Georgiev  [13],  Polk  et  al . [33],  and  Sherenkov  et  al . [39]  have  a 
logarithmic  form  at  the  interfacial  region.  The  graphs  presented  by 
Averkiev  and  Kind  [4]  appear  to  be  parabolic  while  those  by  Bata  and 
Knezevich  [8]  have  a distribution  ranging  from  parabolic  to  logarithmic. 

This  limited  review  leads  to  the  following  conclusions: 

(1)  The  shape  of  the  profile  varies  from  parabolic  to  logarithmic. 
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Figure  2.16  Velocity  distribution  in  a cross  section  of  arrested 
saline  wedge,  after  Keulegan 

Source:  [24,  Figure  11.7] 
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(2)  Csanady's  law  [9]  gives  the  inverse  of  the  coefficient  of  the 
logarithmic  term  equal  to  von  Karman's  constant,  although  Powell's 
[34]  recent  research  seems  to  suggest  that  this  coefficient  must  be  a 
function  of  the  Keulegan  number.  (3)  There  is  no  law  of  variation  for 
the  equivalent  von  Karman  coefficient  nor  for  the  constant  factor  in  the 
logarithmic  law.  The  objective  of  the  work  herein  presented  has  been  to 
provide  some  answers  to  the  specific  problems  and  questions  isolated  in 
this  review  and  to  fill  a number  of  gaps  in  our  knowledge  on  stratified 


flows. 


CHAPTER  III 

THEORETICAL  BACKGROUND 
3.1  Fundamental  Equations 

For  Newtonian  and  incompressible  fluids  with  constant  viscosity 
the  equations  of  conservation  of  momentum,  known  as  Navier-Stokes 
equations,  are 


P 


$ = 

dt 


— TT- 

- VP  + pb 


+ yV2\T 


(3.1.1) 


where  t = u^T  + U2j"  + u3l<  is  the  velocity  vector  with  f,  j",  and  t the 
unit  vectors  in  the  directions  of  the  positive  Cartesian  coordinate  axes 
x-| , x2,  x^,  respectively. 


3x, 


t + 


3x, 


J + 


3x. 


is  the  vector  operator  del , & = b^T  + b2j’  + b^t  is  the  body  force  per 
unit  volume,  p is  the  local  fluid  density,  P is  the  pressure,  and  y is 
the  molecular  viscosity.  Moreover,  if  the  fluid  is  nondiffusive  and 
nonconducti ve , the  equations  of  incompressibility  and  continuity  have 
the  forms 


2f>  = + 

Dt  3t 


(3.1.2) 


V • t = 0 


(3.1.3) 
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In  a gravitational  field  Equations  (3.1.1)  and  (3.1.3)  can  be  written 
in  the  form 


3U.  3U. 

J 


3P  _ 

3x7  " Pgi  + y 3x73x7 
J J 


(3.1.4) 


and 


where  i,  j = 1,2,3,  and  g.  are  the  instantaneous  velocity  and  gravi- 
tational component  in  the  x^  direction. 

Introducing  next  for  turbulent  flows  in  Equations  (3.1.4)  and 
(3.1 .5)  the  notations 

U.  = u.  + u! 
l l l 

P = p + p' 
p = + p' 

where  u^ , p,  are  the  local  time-averaged  values  and  u! , p',  and  p' 
are  the  fluctuating  components  of  U. , P,  and  p,  respectively,  and  taking 
into  consideration  that 

uT  = p1-  = ^ = 0 (3.1.6) 


the  following  equations,  for  steady  state  conditions,  are  obtained: 
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3^ 

3Xj 


+ p'u 


3u! 

' 1 + 

J 3xd 

3u! 

3Xj 

3u . 

+ P'U-  "T “ 

J 3Xj 

UjP' 

3Ui  _ 

. 3]L 

y32u- 

n 4-  1 

J 3x. 


9xi 


3x  -3x . 

J J 


(3.1.7) 


(3.1.8) 


0 


(3.1.9) 


Equation  (3.1.7)  combined  with  Equations  (3.1.8)  and  (3.1.9)  after  sim- 
plification takes  the  form 


P'u: 
u.  + i 

J P 0 


3u 


i + 1 3 


3Xj  Pjt  SXj 


Uj  (p'U,:>  + 


* - 9; 


1 3 


p*  3XJ 


p4ij  + 2p2 


3^ 
3x . 


3u, 


3X-1 


- (uiuj> 


53  5*  + J.  55E3I  + jl 

pt  3Xj  pt  Ui  3Xj  p,  tUlUjJ  3Xj 


3p ' 


The  sum  of  the  three  last  terms  of  the  previous  equation  is  zero.  This 
follows  from  Equations  (3.1.2),  (3.1.8),  and  (3.1.9).  The  final  equa- 


tion then  becomes 
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(3.1.10) 


where 


(3.1.11) 


Here  t. • is  the  stress  tensor.  The  first  term  of  the  right  hand  member 
i J 3 


represents  normal  stresses  due  to  the  thermodynamic  pressure.  The  terms 
in  the  parentheses  represent  the  viscous  stresses,  the  last  term  repre- 


Equations (3.1.8),  (3.1.9),  and  (3.1.10)  are  the  governing  equa- 
tions for  flow  of  any  incompressible  fluid  with  variable  density.  There 
are,  however,  more  unknowns  than  equations.  For  this  reason  we  have  to 
resort  to  assumptions  in  order  to  obtain  appropriate  mathematical  models. 

In  stratified  flows  with  small  rates  of  entrainment  and  distinct 
interface,  the  two  layers  are  normally  considered  as  homogenous  with 
the  exception  of  a zone  of  small  width  around  the  interface.  Moreover, 
if  the  measurements  take  place  at  the  centerline  of  the  channel  and  if 
the  flow  is  considered  symmetrical  without  secondary  currents,  the 
assumption  of  two-dimensional  flow  will  not  deviate  from  reality.  There- 
fore, in  the  outer  region  of  the  interfacial  zone,  the  equation  of  motion 
(3.1.10)  can  be  written  in  the  form 


sents  the  Reynolds  stresses  and  6..  is  the  Kronecker  delta. 

• J 
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3u.  1 3x.. 

U.  t — — ■ = q.  + - — 

J 3x.  yi  P 9x. 


(3.1.12) 


and 


Tij  = ' p5ij  + 2v  2 


3ui 


9u  . 


9x  • 9x . 

J ij 


pu  !u '. 
i J 


(3.1.13) 


In  the  vertical  direction,  the  equation  of  motion  coincides  with 
the  equation  of  hydrostatic  pressure  distribution.  The  equations  of 
motion  and  continuity  in  the  direction  of  flow  have  the  form 


lj!y.+  v— =-g  — ^ _ 1 + 1 + Hi* 

9x  9y  y dx  p 9x  p 9z  9y 


(3.1.14) 


9u 

9x 


+ 


(3.1.15) 


where  u and  v are  the  velocity  components  in  the  flow  direction  x and  in 

the  vertical  direction  y,  respectively,  z is  the  elevation  of  the  bottom 

of  the  channel,  and  x , x are  the  shear  stresses  in  the  direction  x, 

zx  y x 

acting  on  planes  normal  to  z and  y directions,  respectively. 

3.2  The  One-Dimensional  Model 

On  the  basis  of  the  Schijf  and  Schonfeld's  assumption  [37]  and 
referring  to  Figure  3.1,  Equations  (3.1.14)  and  (3.1.15)  for  the  one- 
dimensional flow  obtain  the  forms 


V 


dV 

dx 


dz0  1 dp  1 9x 

dx  ~ p dx  p 9y 


9 


(3.2.1) 
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Figure  3.1  Definition  sketch  of  an  arrested  saline  wedge 


V-jh^  = const. 


(3.2.2) 


where  V is  the  local  mean  cross  sectional  velocity,  in  general,  V-j  is 
its  value  in  the  upper  layer  and  h-j  is  the  thickness  of  the  upper  layer. 
Equations  (3.2.1)  and  (3.2.2)  are  next  applied  to  each  layer  separately 
(Figure  3.1). 

a.  Fresh  water  layer 

The  integration  of  Equation  (3.2.1)  over  any  cross  section  of 
the  upper  layer  gives 


, ^ 

1 dx 


= - g 


dx 


_ 1 d£  _ 
p dx 


ph- 


(3.2.3) 


where  t.  is  the  interfacial  shear  stress.  In  Equation  (3.2.3)  the 
following  relations  are  next  introduced  together  with  Equation  (3.2.2) 


(3.2.4) 


P = P0  + pg(h  - y) 


T 


i 


(3.2.5) 


where  Sq  is  the  bottom  slope  of  the  channel,  y is  the  vertical  distance 
from  the  bed  for  y Z h , Pq  is  a reference  pressure,  h is  the  depth  of 


46 


the  system,  h = h-j  + h , hg  is  the  depth  of  the  lower  layer,  and  p is 
the  density  of  the  fresh  water  layer.  After  some  simplifications. 
Equation  (3.2.3)  obtains  the  form 


dh,  dh  V?  dh,  f.V? 

I + s 1_  _J_  + i 1 _ s 

dx  dx  gh^  dx  8gh-|  " °0 


(3.2.7) 


b.  Salt  water  layer 

Applying  the  same  procedure  as  in  the  upper  layer  and  consider- 
ing that  the  density  difference  between  the  two  layers  is  small,  the 
following  equation  is  obtained 


1 


dh,  dh 

_L  + _ 1 

dx  dx 


V 

8gh 


2 

1 

s 


(3.2.8) 


Here  the  term  of  the  kinetic  energy  and  the  term  of  bottom  shear  stress 
were  neglected,  because,  according  to  the  one-dimensional  model,  Vg— the 
mean  cross  section  velocity  of  the  lower  layer— has  to  be  taken  equal  to 
zero.  Indeed,  the  shear  stresses  at  the  bottom  are  defined  by  the 
relation 


, pV' 

T0  f0  8 


(3.2.9) 


therefore,  if  Vs  is  near  zero,  then,  by  definition,  Tq  has  to  be  taken 
as  zero.  Assuming  next  constant  total  depth,  i.e., 


hQ  = h = hi  + hs  = const. 


(3.2.10) 
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Equations  (3.2.7)  and  (3.2.8)  can  be  combined  in  the  form 


-1  + 


V1 

— Qh. 

0 - I 


dh  -j 
dx 


vl 

8^g 

P 


(3.2.11) 


The  second  term  in  the  parentheses  of  the  left  hand  member  is  the  square 
of  the  densimetric  Froude  number  for  the  upper  layer.  That  is, 


(Fr')2 


( Ap/ p ) gh-j 


Fr-j  can  be  written  in  the  form 


(Frp2 


(Fri)‘ 


(3.2.12) 


(3.2.13) 


Equation  (3.2.11),  after  introducing  Equations  (3.2.10),  (3.2.12),  and 
(3.2.13),  takes  the  form 


8h, 


dx 


(Fpi)' 


31 


(3.2.14) 


Schijf  and  Schonfeld  [37]  considered  the  saline  wedge  as  a long  internal 
wave.  Therefore,  its  length  depends  on  the  condition  at  the  ocean  end, 
that  is  at  x = 0.  The  wedge  becomes  stationary  when  the  conditions 
there  become  critical.  The  critical  depth,  h-j  = h-j  is  given  by  Equa- 
tion (3.2.13)  for  Frj  = 1,  that  is 
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(prJ) 


2/3 


(3.2.15) 


Integration  of  Equation  (3.2.14)  from  -x  to  0 with  the  introduction  of 
Equation  (3.2.15)  at  x = gives 


Ii2< 

8 h. 


3(Fr^)2/3 


1 

4 


10  ^ 


2/31 

J 


hl 

1 

[hl 

5 i 

M 

ho 

StF^)2 

h0j 

4(Fr’)2 

> 

hi 

2thl 

N 

ho 

(3.2.16) 


where  f..  is  the  average  value  of  f.  over  the  length  of  the  wedge,  L. 

That  length  is  the  value  of  x in  Equation  (3.2.16)  for  h^/hg  = L.  That  is 


fjL  = 2h0 


5<Fl0> 


2 - 2 + 3(Fr^)2/3  - | (Fr^)4/3 


(3.2.17) 


The  same  model  and  methodology  applies  to  warm  wedges  [3,  7,  16, 

33]. 

It  should  be  noted  that,  although  the  described  analysis  leads 
to  a handy  practical  model  for  the  intrusion  length  and  the  wedge  shape, 
it  does  not  give  the  detailed  velocity  distribution  nor  the  actual  sig- 
nificant shear  stresses.  Moreover,  there  is  no  way  of  estimating 
analytically  f..  Use,  therefore,  of  Equations  (3.2.16)  and  (3.2.17) 
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would  require  previous  calibration  of  f.  on  the  basis  of  field  and/or 
laboratory  measurements.  Even  so,  it  can  readily  be  shown  that  f.  is 
very  sensitive  to  changes  of  the  flow  parameters.  As  a result,  the  use 
of  these  equations  for  cases  beyond  the  calibration  conditions  becomes 
unreliable. 

3.3  Integration  of  the  Equations  of  Motion: 

Complete  Intrusion 

Equation  (3.1.15)  is  first  multiplied  by  u and  then  added  to 
Equation  (3.1.14)  to  yield 


9u^  + 5(uv)  _ 


3x  9y  ^ d> 


a.  Fresh  water  layer 

Equation  (3.3.1)  integrated  over  the  cross-section  A-,  = bh 


gives 


dA 


(3.3.2) 


The  integrals  are  given  by  the  following  relations: 
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3u 

3x 


dA  = 


fbh,  dvi 


dx 


, „2  dh  , ...2  dhs 
dx  i dx 


(3.3.3) 


where 


aV-j'A.j  = 


u2dA 


A, 


(3.3.4) 


and  a is  the  momentum  correction  coefficient. 


ay 


= bU 


2 dh 
dx 


2 dhs 
bV.  — r— 
l dx 


bvivs 


(3.3.5) 


For  a pressure  given  by 


P = P0  + pg(h  - y)  for  y > hc 


(3.3.6) 


the  pressure  integral  becomes 


JAn 


H dA  ■ *»>i  & 


(3.3.7) 


The  shear  stress  terms  finally  take  the  forms 


1 

P j 


Stzx  2tVw 


3z 


dA  = - 


(3.3.8) 


1 

P J 


dA  = - 1 bT. 
a 8y  p 1 
A1 


(3.3.9) 
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Here  b is  the  width  of  the  channel,  U is  the  velocity  at  the  free  sur- 
face, V.j  is  the  velocity  at  the  interface,  vg  is  the  net  rate  of  salt 

water  entrainment  per  unit  area,  and  x is  the  sidewall  shear  stress. 

r w 

Replacing  the  integrals  by  the  given  functions  in  Equation 
(3.3.2),  the  following  equation  is  obtained: 


jy  dV?  , 

dH_  + _a 1_  + 1 

dx  2g  dx  pgh-j 


Ti  + 


2th, 

w I 


Vs_ 

ghi 


(3.3.10) 


where  dH/dx  is  the  slope  of  the  piezometric  line. 

The  entrainment  velocity  vg  is  equal  to  qs/bl_Q,  where  q$  is  the 
rate  of  salt  water  entrainment  over  the  saline  wedge.  The  latter  has 
been  given  by  Keulegan  [22,  24],  In  the  analysis  of  data,  the  term 
Vi vs/ghl  was  f°uncl  to  be  of  minor  importance  with  a maximum  effect  of 
5 percent.  Because  of  this  and  because  entrainment  does  not  take  place 
continuously  over  the  interface,  but  at  isolated  spots  in  the  form  of 
breaking  of  interfacial  waves,  it  was  neglected  in  the  subsequent  analy- 
sis. 

Next,  expressing  xw  by 


pV‘ 


T — f 

w w 8 


(3.3.11) 


and  solving  Equation  (3.3.10)  with  respect  to  x./p  the  following  equation 


is  obtained: 
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-J-  “ - g(h  - hs)  jH-  - | (h 


dvl  fJh  - Uv? 

hs>  -dF  - - ' 4b  - <3-3'12> 


which  integrated  from  -L  to  0,  gives 


T. 


PV0 


h0  dH  T 
^ v2  dx  1 
0 


ah0  T fw0h0  T 
2L  *2  “ 4b  *3 


(3.3.13) 


Here  "T  and  dH/dx  are  the  average  values  of  t and  dH/dx,  respectively, 
over  the  wedge,  fwg  is  the  wall  friction  coefficient  upstream  of  the 
wedge  and  I-j,  I ^ , Ig  are  dimensionless  functions  given  by  Keulegan  [24], 
after  a numerical  integration  of  his  experimental  data  [Equations 
(2. 4. 1.4),  (2. 4. 1.5),  and  (2. 4. 1.6)  for  I , I.,  I ] . 

a D C 

Partheniades  et  al . [31]  developed  the  following  equation  for 
the  interface: 


x 

L 


Kl 

h 

T 

CO 

4 - 2 

h 

(1  -Fr’)2/Sl 

scj 

l scj 

V 

— 

f h 

4 + 

s 

,hSCj 

1 - (Fr^)Z/3 


r 


6(Fr')4/3  + 3(Fr^)2/3  + 1 


-10 

fh  ] 

s 

h 

2 

(Fr^)4/3  + (Frg)2/3  + 1 

sc 

— — 

6(Fr’)4/3  + 3(Frg)2/3  + 1 


(3.3.14) 
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where  hsc  is  the  critical  depth  of  the  lower  layer  at  the  ocean  end  of 
the  channel,  that  is  at  x = 0. 

Using  the  last  equation,  the  same  investigators  derived  the 
following  functional  relationships  for  the  integrals  I-j , I?,  Ig. 


[5 

(Fr^)2  + 6(Frg)4/3  + 3(Fr^)2/3  + l] 

X1 

3 

6 ( Fr^ ) 4/3  + 3(Frg)2/3  + 1 

Io  = 2 


(FrC)) 


-2/3 


(3.3.16) 


I 


3 


0.3 

5 + 255  (Frg)13/6 

- 221  (Fr^)2  - 39  (Fr')17/6] 

1 - (Fr^)2^ 

3 

"^(Frg)4/3  + 3(Frg)2/3  + f| 

(3.3.17) 


These  relationships  differ  from  those  given  by  Keulegan  [22,  24]  and  the 
difference  increases  as  the  value  of  Frg  decreases,  particularly  below 
0.4  (Figure  3.2).  In  Figure  3.2  Keulegan's  [22,  24]  nondimensional 
curve  for  the  arrested  saline  wedge  is  compared  with  the  theoretical  set 
of  curves  given  by  Equation  (3.3.14).  Keulegan  [22,  24]  used  a limited 
range  of  variation  for  Frg,  something  that  misled  him  to  the  conclusion 
that  the  wedge  shape  is  independent  of  Frg  and  function  of  only  the 
relative  distance,  x/L,  i.e.. 


(3.3.18) 
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Figure  3.2  Dimensionless  shape  of  arrested  saline  wedges 


Source:  [31] 
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For  values  of  x/L  between  0.2  and  0.6,  there  is  satisfactory  agreement 
between  the  theoretical  and  experimental  data,  when  FrQ  takes  the  values 
used  in  Keulegan's  experiments  [22]. 

b.  Salt  water  layer 

The  integral  form  of  Equation  (3.3.1)  over  a cross  section  of 
the  lower  layer,  A$  = bhg , is 


s 


dA 


] lE  dA  + — - — 

p + Ap  jA  3x  p + Ap  JA 


s 


— dA 


(3.3.19) 


where  the  integrals  are  given  by  the  following  relations: 


(3.3.20) 


where 


(3.3.21) 


s 


and  g is  again  the  momentum  correction  coefficient. 
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A 


3^1  dA  = bV2  ^ + bV.vs 


(3.3.22) 


For  a pressure  distribution  given  by 


p = pn  + (h  - hjpg  + g(p  + Ap)  (he  - y)  for  y < h 


the  pressure  distribution  integral  becomes 


(3.3.23) 


1 


p + Ap  , 


A, 


IE.  dA 
3x 


pgbh  dh-.  dh 

~dx  + gbhs  "dT 


(3.3.24) 


Finally,  the  shear  stress  integral  has  the  form 


p + Ap 


3x 


zx 


3z 


dA  + 


^ dA 
3y 


L-  s 


bfTi  + T„) 


P + Ap 


(3.3.25) 


Equation  (3.3.19)  thus  takes  the  form 


d(h_  + zn)  Q d(vfh  ) V . v 


p dH_  + Ap  v s Cr  + J3 


p + Ap  dx  p + Ap 


dx 


ghs  dx 


i s'  + is 


9h< 


Ti  + T0 
g(p  + Ap)h 


= 0 


(3.3.26) 


which,  upon  integration  from  -L  to  0 neglecting  V^Vg/gh  , reduces  to 


+ To 


gh0  dH 


' PV0 


Id  + 


Ap 

S-  h2  + 


Vq  dx  4 ' 2VqL  sc 
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(3.3.27) 


Here  Tg  is  the  average  bed  shear  stress  over  the  length  of  the 
wedge,  (dZg)/(dx)  = - Sg  is  the  average  bottom  slope,  £ = ( ) / ( Vg ) , 

V-jc  is  the  fresh  water  critical  velocity  at  x = 0 and  1^  is  an  integral 
given  by 


The  analysis  to  this  point  is  valid  for  complete  intrusion,  that 
is  for  cases  where  the  salt  water  depth  at  the  ocean  end  of  the  channel 
is  critical  and  equal  to  hsc<  It  is  similar  in  principle  to  the  analy- 
sis used  by  Keulegan  [22,  24]  with  the  following  differences:  (1)  a 

more  detailed  analytical  equation  for  the  wedge  profile  has  been  used 
instead  of  one  single  experimental  curve  and  (2)  the  actual  measured 
values  of  a and  3 for  each  experiment  were  used,  instead  of  the 
constant  values  of  a = 1.028  and  6 = 0.146  used  by  Keulegan  [22,  24]. 
Figure  3.3  gives  graphically  the  integrals  I-j,  Ig , I^.  Moreover,  it 
gives  the  integrals  1^  [Equation  (6.4.1)]  and  Ig,  Ig  [Equations  (A. 3) 
and  (A. 4)  of  the  Appendix]. 


This  is  the  case  where  only  part  Lg  of  the  arrested  saline  wedge 
L is  in  the  channel  (Figure  3.4).  This  case  wi 11  be  di scussed  in  Chapter  V. 


I 


4 


(3.3.28) 


3.4  Integration  of  the  Equations  of  Motion: 
Partial  Intrusion 
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Figure  3.3  Variation  of  the  integrals  I with  respect  to  Frg 
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a.  Fresh  water  layer 

Replacing  in  Equation  (3.3,12)  by  h-h$  and  dividing  both 
2 

members  by  Vg  the  following  equation  is  obtained: 


Ti 


PV 


gh 


1 dH 


ah1  dV^ 


,2  dx 


2Vg  dx 


f h,V? 
w 1 1 

4bV? 


which  can  also  be  written  in  the  form 


gh,  dH  ch0  dV1  fwh, V, 


,2 


\/t-  ..2  dx  i»  dx 

pvo  vo  vo 


4bV, 


In  Equation  (3.4.2),  the  Blasius  equation 


(3.4.1) 


(3.4.2) 


fw  = 


0.316 


Re 


0.25 


(3.4.3) 


is  next  introduced,  where  the  Reynolds  number.  Re,  is  given  by 


Re 


4v]Rhi 

v 


(3.4.4) 


and  where  the  hydraulic  radius  R^  is 


bh- 


hl  (b  + 2h1) 


(3.4.5) 


Thus,  Equation  (3.4.2)  takes  the  form 
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^idH  _ ^0 

.,2  " .,2  dx  Vn  dx 

»vo  vo  0 


0.0559  k2 
,1.25  h0 


V0h0 


0.25 


(b  + 2h]) 


0.25 


(3.4.6) 


The  last  equation,  integrated  from  -Lg  to  0 (Figure  3.4),  takes 


the  form 


Ti  . 9hl  dH  oh0  hs0 


2 .,2  dx  hnn  L r 


pV0  V0 


‘10  '-0 


where 


0.0559  u 

V 

0.25 

(b  + 2h1)0,25r 

,1.25  n0 
b 

V0h0j 

hl 

(3.4.7) 


o J 


h^dx 


(3.4.8) 


and 


(b  + 2h1)0’25' 

1 

r 0 (b  + 2h1)0,25 

hl 

L0 

dx 


(3.4.9) 


Here  T.  is  the  average  value  of  the  interfacial  shear  stresses 
over  the  actual  length  of  intrusion,  Lg.  The  integrals  of  Equations 
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OPEN  CHANNEL 


TANK 


Figure  3.4  Definition  sketch  for  partial  intrusion  of  an  arrested 
saline  wedge 
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(3.4.8)  and  (3.4.9)  have  been  computed  by  numerical  integration  on  the 
basis  of  the  trapezoidal  rule. 

b.  Salt  water  layer 

Equation  (3.3.26),  after  neglecting  the  bottom  slope,  is  first 
written  in  the  form 
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(3.4.10) 


Integration  of  the  last  equation  from  -LQ  to  0 gives 
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where  the  integral 


s " L, 


h dx 
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has  been  computed  by  numerical  integration. 
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c.  Equations  for  T-j  and  Tq 


For  a direct  evaluation  of  "T  and  TQ  by  Equations  (3.4.7)  and 
(3.4.11),  a reliable  measurement  of  the  average  slope  of  the  piezo- 
metric line,  dH/dx,  is  required.  This  was,  in  fact,  the  procedure 
followed  by  Keulegan  in  his  early  studies  [22,  24].  The  experimental 
measurement  of  dH/dx,  however,  has  been  found  to  be  very  difficult, 
since  the  head  loss  over  the  effective  length  of  the  flume  is  normally 
of  the  order  of  a few  tenths  of  a millimeter.  Attempts  to  measure  that 
head  loss  optically  through  a laser  beam  required  an  excessively  long 
time  due  primarily  to  the  slowness  of  achieveing  equilibrium  in  the  two 
static  tanks,  as  well  as  to  a number  of  factors  affecting  the  measure- 
ments. Moreover,  the  above  equations  are  quite  sensitive  to  dH/dx  so 
that  a measurement  error  of  the  order  of  10  percent  may  cause  an  error 
in  T.  of  the  order  of  50  percent. 


distribution  in  the  lower  layer.  Indeed,  the  shear  stress  distribution 
within  the  wedge  is  given  by  the  equation 


The  velocity  head  is  very  small  and  can  be  neglected;  then,  introducing 
p from  Equation  (3.3.23),  we  get 


This  difficulty  was  bypassed  by  assuming  a linear  shear  stress 
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(3.4.13a) 
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The  free  surface  slope  dH/dx  is  very  small,  as  has  been  dis- 
cussed, and  can  be  considered  as  constant.  Also,  apart  from  two 
relatively  small  segments  around  the  wedge  toe  and  the  critical  section, 
the  interface  is  very  nearly  a plane  surface  and  [d(h$  + ZQ)]/(dx)  is 
also  very  nearly  constant.  Therefore,  it  is  reasonable  to  assume  linear 
variation  of  t with  y within  the  wedge.  From  Figure  3.4  it  follows  that 


Ti  + T0 


(3.4.14) 


where  h^  is  the  distance  from  the  interface  to  the  point  of  maximum 
velocity  in  the  wedge.  Assuming  next  that  A is  independent  of  x,  the 
following  simple  relation  for  the  average  shear  stresses  T..  and  TQ  is 
obtai ned 
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(3.4.15) 


Elimination  of  the  mean  slope  of  the  free  surface,  dH/dx  from  Equations 
(3.4.7)  and  (3.4.11)  and  replacement  of  TQ  by  T.. (X  - 1)  in  Equation 
(3.4.15)  leads  to  the  equation 
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(b  + 2h1) 


0.251 


(Ah1  + h$) 


(3.4.16) 


All  the  terms  of  the  second  member  of  Equation  (3.4.16)  are  either  known 
or  can  be  directly  measured.  With  T.  known  the  average  bed  shear  stress 
can  be  evaluated  by  means  of  Equation  (3.4.15). 

For  infinite  channel  width  and  for  Tq  = 0,  Equation  (3.4.16) 

becomes 
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(3.4.17) 


The  corresponding  friction  coefficient,  f . , in  this  case  should 
be  similar  to  the  one  given  by  Schijf  and  Schonfeld's  model  [37]. 

The  analysis  presented  in  section  3.3  constitutes  an  improvement 
as  compared  to  Keulegan's  original  approacn  [22,  24],  particularly  since 
it  eliminates  the  need  of  measurement  of  dH/dx  and  since  it  takes  into 
consideration  partial  intrusion. 

The  expected  relative  error  in  the  estimate  of  X is  consider- 
ably lower  than  the  error  for  dH/dx,  while  the  error  in  the  estimate  of 
shear  stresses  is  smaller  than  that  for  \ [Equation  (3.4.16)]. 


CHAPTER  IV 

EXPERIMENTAL  APPARATUS  AND  INSTRUMENTATION 


4.1  Introduction 

The  experimental  salinity  intrusion  flume  was  designed  at  the 
University  of  Thessaloniki,  in  Greece,  and  at  the  University  of  Florida, 
in  Gainesville,  by  Dr.  Partheniades  and  the  author.  This  apparatus  was 
constructed  at  the  Coastal  and  Oceanographic  Laboratory  of  the  latter 
university. 

The  construction  of  the  flume  started  in  the  summer  of  1977 
and  was  completed  in  the  fall  of  1978,  with  the  exception  of  some  minor 
details  which  were  completed  in  the  beginning  of  1979. 

The  apparatus  basically  consists  of  a variable  slope  plexiglass 
flume,  the  associated  fresh  water  and  salt  water  tanks,  and  the  fresh 
water  circulation  system. 

The  main  measuring  instruments  were  the  hot  film  anemometers,  a 
specially  designed  electronic  conductivity  meter  and  the  orifice  meters. 

4.2  The  Main  Experimental  Apparatus 

In  laboratory  studies  on  salt  water  intrusion  a constant  salin- 
ity water  basin  has  to  be  maintained  at  the  downstream  end  of  the 
channel.  This  is  rather  difficult,  requiring  excessively  large  dimen- 
sions. For  example,  Keulegan  [24]  recommends  a basin  width  equal  to  ten 
times  the  channel  width  and  a basin  length  equal  to  twenty  times  the 
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channel  length.  Space  limitations,  however,  impose  severe  restric- 
tions . 

Therefore,  the  following  apparatus  has  been  designed  and  con- 
structed as  a best  compromise  between  the  ideal  requirements  and  the 
physical  and  technical  limitations.  That  apparatus  is  outlined  in 
Figure  4.1  and  consists  of  a variable  slope  open  channel  with  the 
necessary  tanks  and  other  accessories.  The  channel  has  a width  of  45.7 
cm  and  a length  of  about  20.0  m,  so  that  an  adequate  length  of  salt 
water  intrusion  could  form  inside  the  channel.  Moreover,  the  toe  of 
the  wedge  should  not  exceed  a certain  distance  from  the  upstream  end  of 
the  channel,  otherwise  significant  mixing  could  take  place.  The  results 
of  such  a mixing,  which  has  been  observed  in  some  experiments,  is  the 
development  of  three  or  more  layers  and  the  destruction  of  the  sharp 
interface. 

For  large  Reynolds  numbers,  large  mean  fresh  water  velocities 
and  large  cross  sectional  dimensions  are  required;  however,  physical 
limitations  impose  restrictions  on  both  variables.  Considering  the 
mean  fresh  water  velocity,  there  is  the  additional  restriction  on  the 
interfacial  shear  stresses.  These  stresses  must  be  lower  than  an  upper 
limit,  to  avoid  high  rates  of  entrainment.  In  the  case  of  small 
entrainments  with  still  distinct  interface,  Keulegan  [22,  24]  defined  a 
critical  velocity,  V , for  the  moving  layer.  For  a mean  velocity  larger 
than  the  critical  one,  mixing  is  expected  to  take  place  at  significant 
rates.  In  many  experiments  in  the  present  study  the  mean  fresh  water 
velocity  was  larger  than  Vc;  however,  the  interface  remained  distinct. 
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Finally,  a rectangular  cross  section  45.7  cm  x 45.7  cm  (18  in 
x 18  in)  and  a length  of  20.24  m (66.5  ft)  have  been  selected.  The 
bottom  and  the  walls  of  the  channel  have  been  constructed  of  trans- 
parent plexiglass  1.27  cm  (0.5  in)  thick. 

The  channel  is  supported  on  a steel  frame  (Figures  4.1,  4.2, 

4.3,  and  4.4)  through  a plywood  board.  To  allow  free  expansion  and 
contraction  due  to  temperature  changes,  the  upstream  end  of  the  channel 
was  free  while  the  downstream  end  was  connected  with  the  fixed  tank,  T-j , 
through  a 5 cm  long  thick  rubber  forming  a flexible  and  waterproof 
hinge.  In  addition  to  the  free  slippage  of  the  channel,  this  connection 
also  permitted  a rotation  of  the  channel  around  the  joint,  B (Figure 
4.1),  when  the  slope  was  varied. 

The  channel  has  been  constructed  in  parts  of  2.5  m length, 
placing  a rubber  tape  2 mm  thick  between  them.  For  cross  sectional 
rigidity  vertical  sidewall  supports  and  cross  bars  over  the  channel  have 
been  provided. 

The  steel  frame  is  supported  on  the  joint  B and  on  two  hydraulic 
jacks,  one  at  the  center  and  the  other  at  the  upstream  end  of  the 
channel.  Through  a gear-chain,  the  central  jack  could  be  moved  up  and 
down  over  a height  equal  to  one-half  that  of  the  upstream  jack.  In  this 
way,  a variation  of  the  slope  of  the  channel  up  to  4 °/oo  is  possible 
while  a special  gauge  installed  near  the  upstream  jack  measures  the 
change  of  the  elevation  of  the  bottom  of  the  channel  with  an  accuracy 
of  one  hundredth  of  an  inch. 

Tank  T-j  (Figures  4.1,  4.2,  4.3,  and  4.4)  was  used  as  a salt 
water  tank.  Its  bottom  and  its  three  sidewalls  are  made  of  steel  and 
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Figure  4.1  Main  experimental  apparatus 
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Figure  4.2  Downstream  side  of  the  apparatus 


Figure  4.3  Front  view  of  the  apparatus 
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Figure  4.4  Rear  view  of  the  apparatus 


Figure  4.5  Rear  view  of  the  apparatus  with  gate 
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the  fourth  sidewall  is  made  of  transparent  plexiglass  1.27  cm  (0.5 
in)  thick.  This  permits  a macroscopic  control  of  mixing  of  the  two 
layers  in  the  tank  as  well  as  a control  of  the  salt  water  elevation. 

The  continuous  observation  of  the  interface  level  is  necessary  for  good 
operation.  This  level  should  not  be  lower  than  the  channel  bottom. 

The  tank  was  operated  like  a two-layered  wide  channel  with  the  upper 
layer  moving  over  a quasi-stagnant  lower  layer. 

The  tank  length  is  3 m and  its  depth  is  1.3  m.  The  steel  parts 
of  the  tank  were  covered  with  insulating  coating  and  were  subsequently 
painted  by  plastic  dye. 

To  eliminate  the  disturbance  of  stratified  flow  at  the  entrance 
of  the  tank  due  to  the  abrupt  change  of  the  depth,  a steel  plate  E was 
used  (Figures  4.1  and  4.3).  The  slope  of  the  plate  to  the  horizontal 
direction  is  32.5  degrees.  The  tank  has  two  plexiglass  gates  G-j  and 
(Figures  4.1,  4.2,  and  4.3)  at  its  downstream  and  the  upstream  ends. 

Gate  G^  regulates  the  two-layer  system  in  the  channel.  Pipe  P-|  trans- 
ports fresh  water  to  the  recirculation  tank  T^.  The  upper  end  of  P-j  is 
at  the  level  of  the  bottom  of  the  channel;  therefore,  the  elevation  of 
the  interface  must  be  higher  than  the  upper  end  of  P-j. 

To  avoid  overflow  of  salt  water,  gate  G-j  has  a waterproof  fit 
with  the  tank  wall.  A similar  construction  was  used  for  gate  G^  at  the 
upstream  wall  of  the  tank.  Gate  G^  is  waterproof  and  it  separates  the 
fresh  water  in  the  channel  from  the  salt  water  tank  before  the  beginning 
of  the  experiments.  During  the  experiments  gate  is  raised  to  permit 
the  downstream  flow  of  fresh  water  and  the  upstream  motion  of  the  salt 
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water.  The  raising  and  lowering  of  the  gates  is  carried  out  by  the 
rotation  of  disks  D-j  and  D^.  The  endless  screws  S-|  and  pass  through 
the  disks  D-j  and  respectively,  and  the  screws  have  a fixed  connec- 
tion with  the  gates.  Tank  T-j  is  connected  to  the  emptying  pipe  P^g. 

Tank  (Figures  4.1,  4.2,  and  4.3),  with  a length  of  3.0  m and 
a width  of  1.52  m,  is  used  as  a recirculation  tank  for  the  fresh  water. 
The  tank  is  connected  to  the  axial  flow  pump  A through  the  suction  pipe 
P-|2*  An  appropriate  depth  of  water  must  be  maintained  over  the  crest 
of  P-|2  for  good  operation  without  air  suction.  This  depth  and  the 
diameter  of  the  suction  Dipe  determined  the  depth  of  water  in  T^  to 
0.75  m.  The  steadiness  of  this  depth  is  controlled  by  an  overflow 
pipe  P^.  Tank  T^  has  been  placed  on  in  a way  that  leaves  a part  of 
T^  free.  This  is  necessary  for  macroscopic  control  of  the  steadiness 
of  the  fresh  water  flow.  The  control  of  the  fluctuation  of  the  free 
water  surface  in  T ^ in  every  experiment  was  particularly  important, 
since  this  variation  was  significant  for  short  time  intervals.  In  par- 
ticular, when  the  fresh  water  flow  rate  was  increasing,  the  free  surface 
of  the  system  in  the  channel  also  increased;  however,  the  free  surface 
in  was  lowered.  Therefore,  there  was  a danger  of  a drop  of  the  free 
surface  in  T^  below  the  crest  of  pipe  P^.  To  avoid  this  danger,  the 
volume  of  fresh  water  in  T^  was  replenished  by  tap  water  from  pipe  P^ 
connected  to  the  tap  water  faucet.  It  is  recalled  that,  after  estab- 
lishing a steady  flow,  valve  of  pipe  P^  was  closed.  Tank  T^  is  also 
connected  to  the  emptying  pipe  P^  and  to  pipe  P^  through  which  fresh 
water  is  being  recirculated.  This  tank  is  made  of  steel  and  it  is 
coated  and  painted  in  the  same  way  as  tank  T-j . 
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The  fresh  water  is  pumped  into  the  upstream  tank  Tg  (Figure  4.1) 
from  tank  Tg.  There  the  turbulence  is  dampened  by  flow  through  a honey- 
combed section.  The  recircul ation  of  water  takes  place  through  pipes 
P-J4  of  7.6  cm  diameter  and  P-jg  of  30.4  cm  diameter,  whose  downstream 
ends  are  not  connected  to  the  channel.  Therefore,  the  upstream  end  of 
the  channel  can  be  moved  in  the  vertical  direction  freely.  The  pipes 
are  connected  at  the  location  L and  the  operation  of  one  or  the  other 
depends  on  the  required  discharge.  The  latter  is  regulated  by  the  axial 
flow  pump  A.  Each  one  of  the  pipes  has  an  orifice  meter.  Orifice  0-| 
with  a diameter  D = 2.5  cm  is  in  pipe  P14  and  orifice  0g  with  a diameter 
D = 7.6  cm  is  in  the  pipe  P-jg. 

Pipes  P-j , Pg,  Pg,  and  Pg  and  valves  V-j,  Vg,  Vg,  and  have  an 
appropriate  connection  to  drive  the  recirculated  fresh  water  either  to 
tank  Tg  or  to  the  sewer,  through  pipe  P-j-j,  when  the  mixing  of  the  two 
layers  is  significant. 

Pipe  Py  is  branched  into  pipes  Pg  and  Pg  by  an  appropriate 
operation  of  valves  Vg  and  Vg.  Pipe  P -j  either  drives  salt  water  from 
the  tank  of  preparation  of  salt  water  to  tank  T-j  through  pipe  Pg,  or 
it  empties  the  tank  of  preparation  of  saltwater  through  pipe  Pg.  In 
the  last  case  the  salt  water  is  driven  through  the  pipes  Py,  Pg,  P^g, 
and  P-j^.  The  length  of  pipe  Pg  was  so  chosen  that  the  system  of  pipes 
Pg  and  Pg  operates  as  a reversed  siphon.  Therefore,  to  avoid  the  danger 
of  salt  water  intrusion  reaching  the  upstream  end  of  the  channel,  the 
level  of  salt  water  in  Tank  T^  was  lowered  using  the  reversed  siphon. 

In  this  way,  salt  water,  instead  of  being  transferred  to  the  sewer,  was 
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stored  in  the  mixing  tank  for  a new  usage  during  the  experiment. 

The  salt  water  from  the  mixing  tank  was  transferred  to  tank  T-j  by  a 
pump. 

The  measuring  instruments,  that  is  the  hot  film  anemometer 
probe,  the  conductivity  meter  probe  and  the  point  gauges,  were  supported 
on  a carriage  moving  on  rails  (Figure  4.1). 

4.3  Instrumentation 

The  electronic  system  of  instruments  used  in  this  study  is  about 
the  same  as  the  one  used  by  Powell  [34].  For  this  reason  much  of  this 
analysis  on  the  instrumentation  is  basically  a repetition  of  that  given 
by  Powell  [34].  For  more  information  see  reference  [34],  This  system 
contained  transducers,  signal  conditioners,  and  an  analog-to-digi tal 
(A/D)  tape  recorder.  The  output  signal  from  each  transducer  was  passed 
through  one  or  more  signal  conditioners  and  was  sampled  and  recorded 
into  magnetic  tape  by  the  kinemetrics  Model  DDS-1103  A/D  tape  recorder. 
Figure  4.7  shows  the  general  arrangement  of  the  instrument  package. 

The  input  voltage,  being  in  A/D  recorder  between  -5.0  V and 
+5.0  V,  is  resolved  into  4,096  discrete  values.  Therefore,  the  input 
resolution  for  this  unit  was  limited  to  t 0.0012  V.  The  probes  measur- 
ing temperature,  conductivity,  and  velocity  were  mounted  at  the  end  of 
the  transverse  rod  of  the  point  gauges.  All  this  system  of  probes  and 
gauges  were  supported  on  the  carriage.  Figure  4.8  gives  a general  view 
of  electronic  instrumentation  used  in  this  study. 

The  measurement  of  temperature  in  the  vicinity  of  the  conductiv- 
ity and  hot  film  probes  was  obtained  with  the  thermistor  probe.  That 
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Figure  4.6  Rear  view  of  the  apparatus  showing  the  laser  instrument 


TEMPERATURE  TRANSDUCER: 
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Figure  4.8  Outline  of  electronic  instrumentation 
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probe,  Model  SF2-42A68,  manufactured  by  Victory  Engineering  Company, 
has  a sensor  of  1.5  mm  in  diameter.  The  output  from  the  thermistor  was 
passed  through  a low-pass  filter  with  a cutoff  frequency  of  20  Hz.  The 
output  from  the  filter  was  recorded  by  the  A/D  recorder. 

The  density  at  any  point  of  the  cross  section  of  the  stratified 
system  was  determined  by  the  conductivity  and  temperature.  The  conduc- 
tivity meter  used  in  this  study  was  designed  by  Doddington  [38]  and 
constructed  in  the  Engineering  Sciences  Laboratory  of  the  University  of 
Florida  in  Gainesville  (Figure  4.9).  It  is  made  of  acrylic  plastic  and 
has  one  stainless  steel  and  one  platinum  electrode.  During  its  opera- 
tion, the  surrounding  fluid  is  continuously  drawn  into  the  interior 
chamber  through  a 0.5  mm  slit  and  a 0.3  mm  diameter  vertical  hole  at  the 
bottom  of  the  probe.  The  fluid  exits  through  the  plastic  tube  connect- 
ing the  probe  and  the  vacuum  pump.  The  electrical  current  density 
between  the  electrodes  is  small  except  in  the  vertical  hole,  where  it 
becomes  large.  Thus,  the  resistance  between  the  electrodes  at  any 
instant  is  governed  by  the  conductivity  of  the  fluid  in  the  vertical 
hole,  not  by  the  resistance  at  the  fluid-electrode  interface.  This 
reduces  the  effects  caused  by  electrode  contamination  and  pla-ting.  For 
more  details  about  the  conductivity  meter  see  Sheppard  and  Doddington 
[38].  The  vacuum  pump  should  work  on  low  suction  pressures  to  have 
small  flow  rates  through  the  probe.  An  indicative  maximum  flow  rate  is 

3 

about  30  cm  /min,  which  does  not  affect  the  conductivity. 

The  output  voltage  from  this  transducer  becomes  again  the  input 
to  a low-pass  filter  with  a 20-Hz  cutoff  frequency  and  the  output  from 
the  filter  is  also  recorded  digitally  on  a magnetic  tape. 
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PUMP 


Figure  4.9  Cross  section  of  conductivity  probe 
Source:  [38] 
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A hot  film  anemometer  with  constant  temperature,  Thermo-system 
Model  1050  general  purpose  anemometer,  was  used  to  measure  both  the  mean 
velocity  and  the  fluctuating  components  of  velocity  in  the  plane  of  the 
sensors.  The  two  sensors  are  part  of  an  x-type  hot  film  probe.  Thermo- 
system Model  1243-60w  NaCl . 

The  output  from  each  anemometer  is  first  passed  to  a thermosys- 
tems Model  1057  signal  conditioner.  This  conditioner  suppresses  the 
excessive  DC  voltage  to  make  the  signal  compatible  with  the  input 
requirements  of  the  A/D  recorder.  From  the  conditioner  the  signal  is 
sent  to  a low-pass  filter  with  a cutoff  frequency  of  17  Hz  and  from 
there  the  signal  goes  to  the  A/D  recorder  where  it  is  ditigized  and 
recorded  onto  a magnetic  tape. 

The  overheating  ratio  is  given  as 

r ' 1 + “c(ts  - ‘e> 


where  is  the  operational  resistance  of  the  sensor,  Rg  is  the  sensor 
resistance  in  environmental  temperature,  aQ  is  the  temperature  coef- 
ficient, tg  is  the  environmental  temperature,  and  t is  the  operational 
temperature  of  the  sensor.  For  the  hot  film  used  in  this  study  the 
value  of  ac  is  0.0026/°C,  To  avoid  creation  of  bubbles  on  the  sensors 

by  the  high  temperature,  an  overheat  ratio  of  Ru/R  = 1.05  was  chosen. 

h e 

To  transform  the  output  voltage  into  flow  velocity,  King's 
law  has  been  used.  That  law  reads 


E 


2 


A + B(pu)1/n 


(4.3.2) 
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where  A and  B are  constants  depending  on  the  temperature,  u is  the  fluid 
velocity,  1/n  was  taken  equal  to  0.45,  and  E is  the  output  voltage. 

The  low-pass  filters  were  used  to  prevent  contamination  of  the 
signal  by  high  frequency  noise,  particularly  60-cycles  noise,  which  is 
picked  up  from  the  large  pumps,  electric  motors,  and  other  support 
equipment  in  the  laboratory.  Spectral  analysis  of  the  velocity  data, 
done  by  Powell  [34],  has  indicated  that  for  these  experiments  negligible 
energy  exists  above  7 Hz;  therefore,  these  low  pass  filters  have  little 
or  no  effect  on  the  true  signal  from  their  respective  transducers. 

In  addition  to  the  aforementioned  electronic  instrumentation, 
some  other  instruments  were  used,  such  as  orifices,  thermometers , and 
pycnometers. 


CHAPTER  V 

EXPERIMENTAL  PROCEDURE  AND 
CALIBRATION  OF  THE  INSTRUMENTS 

5.1  Introduction 

The  hydrodynami c behavior  of  arrested  saline  wedges  has  been 
studied  by  using  fresh  and  salt  water,  with  the  latter  dyed  by  a red 
food  color  dye.  The  relative  density  difference,  Ap/(p  + Ap),  ranged 
from  8 °/oo  to  93  °/oo. 

The  following  primary  parameters  were  independently  varied: 

(1)  the  fresh  water  flow  rate,  Q,  (2)  the  ratio  Ap/(p  + Ap),  and  (3)  the 
fresh  water  depth,  hg,  upstream  of  the  wedge. 

Each  group  of  experiments  was  conducted  during  a day  and  was 
preceeded  by  the  following  stages  of  preparatory  work: 

1.  Preparation  of  the  salt  water  in  the  mixing  tank.  The 
salt  water  was  subsequently  transferred  to  tank  T-j  up  to 
the  desirable  depth  (Figure  4.1).  Gate  was  closed 
during  this  phase 

2.  Calibration  of  the  electronic  instruments.  When  the 
measurement  duration  was  more  than  10  to  12  hours,  a 
checking  of  the  instruments  calibration  was  necessary, 
particularly  for  the  hot  film  sensors  and  the  thermistor 

3.  Placing  and  fixing  of  the  probes  on  the  carriage  K 
(Figure  4.1) 
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4.  Gate  G^  was  raised  and  the  operation  of  axial  pump 
A was  started  at  the  same  time  (Figure  4.1).  The 
measurements  started  a few  hours  after  the  onset  of 
flow  to  allow  the  establishment  of  the  maximum  degree  of 
steadiness . 

The  preliminary  work  usually  started  at  seven  o'clock  in  the  morning  and 
ended  early  in  the  afternoon.  The  main  experiments  followed  right 
afterwards.  The  fresh  water  flow  rate  was  measured  through  orifices 
0-]  and  O2 ; next,  the  total  depth  of  the  two  layers,  the  length  of  salt 
water  intrusion,  and  the  depth  of  the  saline  wedge  at  every  meter  of 
length  were  determined.  Finally,  the  instantaneous  velocities,  conduc- 
tivities, and  temperatures  at  several  points  in  a selected  cross-section 
were  measured. 


5.2  Preparation  of  the  Salt  Water  in  Tank  T] 

The  fresh  water,  the  salt,  and  the  dye  were  mixed  in  the  mixing 
tank  until  a uniform  salt  water  mixture  of  the  desirable  density  was 
obtained.  The  mixing  tank  volume  was  about  one  cubic  meter.  For  this 
reason  and  before  starting  the  experiments,  some  time  was  necessary  to 
fill  tank  T-|  up  to  the  required  depth.  A volume  equal  to  seven  times  the 
volume  of  the  mixing  tank  was  required. 

According  to  Equations  (3.4.16)  and  (3.4.17),  the  length  of  salt 
water  intrusion  in  the  channel,  Lg,  is  an  increasing  function  of  the 
depth  at  the  mouth  of  the  channel,  h-jg,  and  a decreasing  function  of  the 
average  interfacial  shear  stress  T. . Therefore,  an  approximate  estimate 
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of  hsQ  was  necessary  so  that  LQ  would  be  less  than  the  length  of  the 
channel.  This  was  made  possible  by  a rough  calculation  of  an  arrested 
saline  wedge  with  a length  of  intrusion  equal  to  LQ  (Figure  3.4). 

For  values  of  hg  equal  to  or  less  than  10  cm,  the  levels  of 
free  surfaces  of  the  fresh  and  salt  water  on  the  respective  side  of  gate 
G2  were  approximately  the  same.  For  values  of  hQ  higher  than  15  cm  the 

fresh  water  level  was  higher  than  the  salt  water  level  by  a depth  h . 

9 

This  was  done  in  order  to  avoid  any  disturbance  of  the  flow  at  the  posi- 
tion of  gate  at  the  start  of  the  experiment  due  to  an  unbalanced 
pressure  force.  The  operation  in  this  case  was  the  following:  The 

depth  hg  was  approximately  estimated  and  tank  T^  was  filled  with  salt 
water  up  to  the  bottom  of  the  channel.  Then  the  axial  pump  was  started 
at  very  small  discharge  to  fill  the  channel  and  the  tank  T-j  to  a depth 
hg  above  the  bottom  of  the  channel.  During  this  operation  valves  V-j  and 
V2  remained  closed,  gate  G2  remained  open,  and  gate  G-j  was  lifted  to  a 
height  determined  on  the  basis  of  the  desirable  depth  hg  upstream  of  the 
wedge.  Fresh  water  was  replaced  in  tank  T2  during  pumping  through  pipe 
P-|2.  When  the  fresh  water  depth  in  tank  T-|  was  equal  to  hg,  gate  G2 
was  closed  while  fresh  water  was  supplied  to  the  channel  and  salt  water 
was  supplied  to  the  tank  T-j  from  the  mixing  tank  through  pipe  Pg.  The 
discharge  of  pipe  Pg  was  small  in  the  beginning  to  avoid  any  interfacial 
mixing  and  it  was  increased  later  on.  When  the  depth  of  fresh  water  was 
approximately  equal  to  the  desirable  depth,  any  supply  of  fresh  or  salt 
water  was  ceased.  To  avoid  any  mixing  of  fresh  water  in  the  channel  and 
the  salt  water  in  tank  T-j  due  to  a possible  small  leakage  through  gate  G2, 
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another  gate  was  installed,  as  shown  in  Figure  4.1.  The  thermistor 
and  the  hot  film  sensors  were  next  calibrated. 

5.3  Calibrations 

The  thermistor  was  calibrated  through  a mercury  thermometer  with 
a resolution  of  ± 0.05  degrees  Celcius  (°C),  assumed  to  be  the  same  as 
the  resolution  of  the  thermistor. 

Each  calibration  curve  corresponds  to  a daily  series  of  exper- 
ments.  In  the  curve  of  Figure  5.1  a small  degree  of  nonlinearity  can 
be  observed.  However,  the  linearization  of  the  curve  on  the  basis  of 
the  method  of  least  squares  error  gives  an  error  less  than  + 0.1°C. 

The  hot  film  anemometer  sensors  were  calibrated  through  a Thermosys- 
tems, Inc.  calibrator  for  water,  Model  1127. 

Figure  5.2  gives  typical  calibration  curves  for  the  two  sensors 
of  an  experimental  day. 

The  method  of  least  squares  error  was  applied  again  for  the  cor- 
relation curves  and  the  square  of  correlation  coefficient  was  of  the 
order  of  99.9  percent.  Otherwise,  the  relative  error  between  the  real 
velocities  and  the  corresponding  values  on  the  calibration  curve  would 
be,  at  least,  of  the  order  of  10  percent.  It  follows  from  Figure  5.2 
that  for  small  velocities,  the  correlation  is  nonlinear;  however,  we  can 
accept  the  linearity  for  velocities  as  low  as  3.0  mm/sec. 

Figure  5.3  gives  a partial  view  of  the  instruments  used  for  the 
calibration  of  the  thermistors  and  the  sensors. 

The  conductivity  probes  were  calibrated  just  prior  to  the  exper- 
imental work  and  the  calibration  curves  were  frequently  checked  during 
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Figure  5.1  Calibration  curve  for  the  thermistor 
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the  conduct  of  each  experiment.  The  calibration  was  repeated  after 
five  months  in  order  to  detect  the  time  effect  on  the  probe.  That 
effect  was  found  to  be  negligible. 

The  calibration  started  with  the  correlation  of  the  output  of 
the  probe  with  the  relative  density  difference  Ap/(p  + Ap)  and  with 
temperature  (Figure  5.4).  The  ratio  of  Ap/(p  + Ap)  in  the  first  cali- 
bration ranged  from  5 °/oo  to  40  °/oo  and  in  the  second  one  from  5 °/oo 
to  85  °/oo.  The  temperature  ranged  from  14.5°C  to  29°C.  The  variation 
of  Ap/ (p  + Ap)  with  the  temperature  was  approximately  linear.  Through 
the  application  of  the  method  of  least  squares  error,  the  following 
relation  was  derived  first  for  a temperature  of  20.4°C  (Figure  5.5) 
and  later  on  extended  to  any  temperature. 


. AP 

(p  + Ap) 


= - 0.0031  + 0.087E  - 0.00067  (T  - 20.4) 


(5.3.1) 


where  E is  the  probe  output  in  volts,  T is  the  temperature  in  degrees 
centigrade,  and  Ap/ (p  + Ap)  is  given  per  thousand  (°/oo).  Equation 
(5.3.1)  is  applicable  only  for  Ap/(p  + Ap)  = 40  °/oo.  The  calibration 
curve  for  Ap/(p  + Ap)  larger  than  40  °/oo  is  given  in  Figure  5.6.  In 
this  figure,  the  calibration  curve  for  Ap/(p  + Ap)  > 40  °/oo  has  the 
form 


..  AP 

(p  + Ap) 


- 0.0021  + 0.125E  - 0.00067  (T  - 20.4) 


(5.3.2) 


and  the  calibration  curve  for  Ap/(p  + Ap)  ^ 40  °/oo  has  the  form 


= - 0.0032  + 0.0875E  - 0.00067  (T  - 20.4)  (5.3.3) 
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Figure  5.2  Calibration  of  hot  film  sensors 
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Figure  5.3  Instruments  used  for  calibration  of  thermistor  and  sensor 
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Figure  5.4  Calibration  curves  for  the  conductivity  probe  for  various 
temperatures 
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Figure  5.5  Calibration  of  the  conductivity  probe  for  T = 20.4°  C 
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Figure  5.6  Calibration  of  the  conductivity  probe  for  Ap/(p  + Ap)  up  to 
85  °/oo 
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The  difference  between  Equations  (5.3.1)  and  (5.3.3)  appears  to  be 
neglibible.  Therefore,  the  time  had  no  serious  effect  on  the  conduc- 
tivity probe.  The  curves  of  Figure  5.6  have  different  slopes  because 
the  conductivity  meter  had  been  designed  for  Ap/(p  + Ap)  ^ 40  °/oo. 

The  calibration  curves  for  the  orifice  meters  were  given  by  the 
company;  however,  they  were  checked  in  place.  The  continuous  curve  in 
Figure  5.7  is  the  company's  calibration,  while  the  points  correspond  to 
the  experimental  check.  After  linearizing  the  given  curves,  application 
of  the  method  of  the  least  squares  error  led  to  the  calibration  formulas 

Q = 0.0048(Ah)0-488  (5.3.4) 

for  the  orifice  with  a diameter  D = 2.5  cm  and 

Q = 0.0711 (Ah)0,499  (5.3.5) 

for  the  orifice  with  a diameter  D = 7.6  cm.  Here  Ah  is  expressed  in  meters 

3 

of  mercury  column  and  Q in  m /sec. 

A figure  similar  to  Figure  5.7  was  obtained  for  the  orifice  with 
diameter  D = 7.6  cm. 

5.4  The  Placing  of  the  Instruments  on  the  Carriage 
After  their  calibration,  the  thermistor  and  the  hot  film  probe 
were  attached  to  the  shaft  of  one  point  gauge  while  the  conductivity 
probe  was  attached  to  the  shaft  of  a second  point  gauge.  These  gauges 
were  permanently  supported  on  the  carriage  K (Figure  4.1).  The  horizon- 
tal distance  between  the  hot  film  probe  and  the  thermistor  was  about  2 cm 
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Figure  5.7  Orifice  calibration  curve  for  D = 2.5  cm 
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and  that  between  the  hot  film  and  conductivity  probe  was  about  8 cm. 

Thus,  the  high  frequency  fluctuations  of  temperature  at  the  position  of 
the  thermistor  and  of  conductivity  at  the  position  of  the  conductivity 
probe  were  not  affected  by  the  high  frequency  fluctuations  of  tempera- 
ture and  density  in  the  neighborhood  of  the  hot  film  sensors.  The 
existence  of  such  fluctuations  introduces  nonexisting  fluctuations  in  the 
velocity  measurements. 

Because  temperature,  density,  and  velocity  varied  with  the  depth, 
the  measurements  were  taken  at  the  same  horizontal  level.  First,  these 
instruments  were  leveled  with  respect  to  the  horizontal  surface  of  stag- 
nant water.  Next,  the  level  of  the  bottom  of  the  channel  was  determined. 
After  all  these  preparatory  steps,  the  experiments  were  ready  to  start. 

Figure  5.8  shows  the  establishment  of  the  level  of  the  bottom  of 
the  channel  with  only  fresh  water  in  it  before  starting  the  experiments. 

5.5  Experiments 

The  outlined  preparatory  steps  usually  took  about  six  hours, 
provided  that  no  problems  were  encountered  in  the  establishment  of  the 
calibration  curves  of  the  instruments.  One  of  the  usual  problems  was 
the  formation  of  very  small  air  bubbles  invisible  to  the  eye  on  the  x- 
type  hot  film  probe.  The  existence  of  the  bubbles  was  detected  after 
the  definition  of  calibration  curve,  when  the  square  of  correlation 
coefficient  of  the  curve  was  less  than  99.8  percent.  In  such  a case,  the 
calibration  of  the  hot  film  sensors  was  repeated.  For  this  reason,  the 
preliminary  work  sometimes  continued  for  more  than  eight  hours.  When 
the  air  bubbles  were  visible,  the  required  time  was  not  so  long. 
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Figure  5.8  Establishment  of  bottom  elevation  of  the  channel  before 
starting  the  experiments 
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When  everything  was  ready,  the  movable  gate  G^  was  raised  and 
the  axial  pump  A was  turned  on  with  a nearly  simultaneous  raising  of 
gate  G2  and  an  opening  of  valves  V-j  and  V2.  Valve  V4  was  closed.  In 
this  way,  all  the  fresh  water  discharge  was  driven  to  the  sewer  through 
the  pipes  , P2,  Pg,  and  P-p.  This  was  necessary  because  by  raising 
gate  G2  a salt  water  front  was  formed  while  moving  upstream  (Figure 
5.9).  The  interfacial  mixing  of  the  front  was  significant  because  the 
fresh  water  conditions  had  not  yet  reached  their  final  steady  state. 

Thus,  the  fresh  water  contained  large  quantities  of  salt  and  that  vio- 
lated the  requirement  of  small  entrainment  and  sharp  interface.  The 
mixing  was  normally  controlled  through  visual  observation  and  in  certain 
instances  through  sampling  from  the  recirculation  tank  T2 . The  volume  of 
fresh  water  in  T2  was  replaced  through  the  discharge  pipe  P-j 2 , which  was 
connected  to  the  tap  water  faucet. 

In  order  to  maintain  the  toe  of  the  saline  wedge  downstream  of 
the  head  water  tank,  the  fresh  water  flow  rate  and  the  depth  of  salt 
water  at  the  mouth  of  the  channel  had  to  be  properly  controlled.  The 
flow  rate  was  controlled  through  the  axial  variable  speed  pump  and  the 

l 

salt  water  depth  at  the  mouth  of  channel  was  controlled  through  siphon 
P8"P7* 

As  it  was  pointed  out  earlier,  the  length  of  salt  water  intru- 
sion was  very  sensitive  to  flow  and  salinity  fluctuations;  for  this 
reason,  the  pump  and  siphon  operations  had  to  be  done  very  carefully. 

When  the  upper  layer  was  clear,  valves  V2  (emptying  the  water  to  the 
sewer),  Vg  (operating  the  siphon),  and  (discharging  tap  water)  were 
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Figure  5.9  Upstream  motion  of  salt  water  front 


Figure  5.10  Example  of  smooth  interface 
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closed  while  valve  (through  which  the  upper  layer  was  transferred  to 
tank  T2  for  circulation)  was  opened.  Next,  the  flow  rate  for  each 
experiment  was  established  through  the  axial  pump  (Figure  4.1).  This 
work  was  difficult  because  the  motion  of  salt  water  front  had  to  be 
continuously  observed  and  the  discharge  of  pump  A had  to  be  properly 
controlled.  The  establishment  of  steady  conditions  required  normally 
one  hour  and  in  certain  special  and  extreme  cases  more  than  four  hours. 

The  entrainment  of  salt  water,  although  small,  was  eroding  the 
depth  of  salt  water  and  diminishing  the  length,  Lg,  of  the  saline  wedge. 
For  small  interfacial  shear  stresses,  the  decrease  of  Lg  was  negligible; 
however,  for  larger  stresses,  the  erosion  of  the  salt  water  layer  was 
large  and  the  assumption  of  stationary  saline  wedge  was  no  longer  valid. 
To  avoid  this,  Tank  T-j  was  continuously  supplied  with  salt  water  in 
appropriate  discharges  through  pipe  Pg.  For  smooth  interfaces  (Figure 
5.10)  and  low  mean  velocity  in  the  upper  layer,  the  interfacial  shear 
stresses  were  low  and  mixing  was  small;  therefore,  the  depth  of  salt 
water  in  any  cross  section  was  approximately  constant. 

When  the  duration  of  the  experiment  was  of  the  order  of  half  an 
hour,  the  decrease  of  intrusion  length,  Lg,  was  negligible  compared  to 
the  original  length.  Most  of  the  experiments  of  series  I,  L,  N,  ,0,  P, 
and  many  of  the  remaining  series  belong  to  this  category. 

When  the  mean  velocity  of  the  upper  layer  was  large,  interfacial 
waves  appeared  (Figure  5.11).  For  even  larger  velocities,  these  waves 
started  breaking.  The  upper  limit  of  the  mean  velocity  of  the  upper 
layer  was  established  from  the  form  of  wave  breaking.  Figure  5.12  gives 
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Figure  5.11  Example  of  wavy  interface 
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V, 


Figure  5.12 


(a)  Eddy  formation  due  to  the  interfacial  wave  breaking; 

(b)  beginning  of  gradual  mixing  of  the  two  layers 


Source:  [14] 
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two  characteristic  forms  of  breaking.  In  the  first,  for  moderate 
interfacial  shear  stresses,  a small  volume  of  injected  fluid  is  diffused 
in  the  upper  layer,  but  the  largest  part  of  it  is  returned  back  to  the 
saline  layer.  In  this  case  mixing  is  a factor  of  stability.  In  the 
second  case  of  high  interfacial  shear  stresses,  the  total  injected 
fluid  in  the  upper  or  in  the  lower  layer  is  almost  totally  diffused. and 
mixed  with  the  ambient  water.  This  is  the  beginning  of  gradual  mixing 
of  the  two  layers.  In  the  experiments  of  the  present  study  the  upper 
limit  of  mean  velocity  of  the  fresh  water  layer  was  determined  by  the 
instability  of  the  first  form  (Figure  5.12). 

The  following  flow  and  salinity  parameters  were  measured: 

(1)  the  depths,  h , of  salt  water  intrusion  in  every  meter  of  its  length; 

(2)  the  salt  water  intrusion  length,  Lq;  and  (3)  the  manometer  readings. 
These  preliminary  measurements  were  followed  by  the  detailed  measurements 
of  velocities,  temperatures,  and  salinities  at  several  points  of  the 
cross  section  from  the  bottom  to  the  free  surface  by  the  three  instru- 
ments attached  to  the  point  gauges.  One  of  the  point  gauges  was  used  to 
establish  the  levels  of  the  free  surface  and  of  the  bottom,  from  which 
the  depths  of  each  layer  of  the  system  were  calculated.  Next,  the  posi- 
tion of  the  thermistor,  the  hot  film  sensors,  and  the  conductivity  probe 
were  defined  by  the  point  gauges.  The  A/D  recorder  was  turned  on  to 
record  the  data  of  velocities  and  conductivities  on  the  tape.  The 
measurements  took  place  in  a cross  section  within  the  middle  third  of  the 
saline  wedge.  The  recording  was  automatically  stopped  before  continuing 
the  measurements  at  another  point  of  the  cross  section.  Most  of  the 
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points  were  located  around  the  interface  and  in  the  lower  layer.  In 
the  lower  part  of  the  salt  water  layer  the  sensors  were  directed  down- 
stream, because  in  this  part  the  salt  water  motion  was  in  the  upstream 
direction.  Above  this  part  of  the  lower  layer  the  sensors  were  directed 
upstream. 

For  every  point  of  the  cross  section,  four  series  of  data  were 
recorded  on  the  tape,  that  is,  the  data  of  the  thermistor,  of  the  con- 
ductivity meter,  and  of  the  two  hot  film  sensors.  Using  the  calibration 
curves,  the  I.B.M.  electronic  computer  was  transforming  the  data  from 
volts  to  data  in  temperature,  density,  and  velocities,  respectively.  In 
particular,  to  calculate  the  instantaneous  velocities  at  a point,  the 
temperature  and  the  density  at  the  point  had  to  be  known  [Equation 
(4.3.2)]. 

As  mentioned  in  the  previous  chapter,  the  high  frequency  fluctua- 
tions in  the  temperature  and  in  the  conductivity  may  introduce  fluctu- 
ations in  velocity  measurements  which  do  not  actually  exist.  To  reduce 
the  effect  of  this  source  of  errors,  the  time  average  values  of  tempera- 
ture and  density  were  introduced  in  Equation  (4.3.2)  instead  of  their 
instantaneous  values. 

In  order  to  estimate  the  magnitude  of  errors  from  such  a source, 
Powell  [34]  performed  the  following  simple  experiment.  He  determined  the 
instantaneous  velocity  and  all  the  correlations,  first  for  constant 
temperature  and  density  equal  to  their  corresponding  averaged  values  and 
next  for  instantaneous  values  of  temperature  and  density.  The  comparison 
of  the  results  showed  that  the  relative  difference  between  these  two  sets 
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of  values  was  always  less  than  0.1  percent.  Therefore,  the  temperature 
and  density  fluctuations  are  not  considered  to  be  a significant  source 
of  error. 


CHAPTER  VI 

ANALYSIS  OF  EXPERIMENTAL  DATA 
6.1  Introduction 

There  were  two  objectives  that  directed  the  analysis  of  the 
experiments:  firstly,  the  derivation  of  universal  laws  for  the  inter- 

facial and  bed  friction  coefficients,  f ..  and  fQ,  respectively,  and, 
secondly,  the  derivation  of  similar  laws  for  the  velocity  and  density 
distribution,  particularly  in  the  neighborhood  of  the  interfacial  zone. 

The  main  flow  and  salinity  parameters  ranged  within  the  follow- 
ing limits:  (a)  as  stated  in  the  previous  chapter,  the  relative  density 

difference  ratio  Ap/(p  + Ap)  ranged  between  8 °/oo  and  93  °/oo, 

(b)  the  mean  velocity,  Vg,  upstream  of  the  toe  of  the  wedge  varied 
from  1.86  cm/sec  to  10.01  cm/sec  and  the  velocity  of  the  upper  fresh 
water  layer,  V-j,  in  the  cross  section  of  measurements  ranged  from  2.21 
cm/sec  to  10.91  cm/sec,  (c)  the  depth  of  flow  upstream  of  the  toe  of 
the  wedge,  hg,  ranged  from  6.55  cm  to  34.57  cm,  and  (d)  the  Reynolds 
number,  upstream  of  the  toe  of  the  wedge,  based  on  the  equivalent 
pipe  diameter,  ranged  from  6,700  to  62,000. 

The  cross  section  of  velocity  and  density  measurements  was 
located  at  about  4.5  meters  from  the  downstream  end  of  the  channel. 

6.2  Density  Distribution  Law 

In  each  position  of  the  conductivity  probe,  more  than  two 
hundred  instantaneous  data  of  density  were  taken.  These  data  were 
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averaged  to  obtain  the  mean  value  of  density  at  this  position.  The 
variation  of  density  with  the  depth  was  significant  in  the  vicinity  of 
the  interface,  becoming  negligible  at  some  distance  away  from  it,  of  the 
order  of  a few  millimeters.  For  this  reason,  most  of  the  measurements 
were  restricted  to  the  immediate  neighborhood  of  the  interface.  Knowing 
the  position  and  the  corresponding  density,  the  density  profiles  have 
been  drawn  for  each  experiment.  Figures  A.l  to  A. 5 of  the  Appendix  give 
examples  of  density  profiles  for  some  experiments  belonging  to  the  one 
particular  group.  Figure  6.1  gives  the  parameters  of  a density  inter- 
facial  zone.  The  width  6 of  that  zone  is  defined  as  the  vertical 
distance  between  the  intersection  points  of  the  vertical  lines  at 

AV(p  + Ap)  ’ 0 and  AV(p  + Ap)  " Ap/(p  + Ap)  and  the  tangent  to  the 
curve  at  Ap/2(p  + Ap),  where  the  subscript  z indicates  mean  local  condi- 
tions with  respect  to  time.  The  assumption  of  small  entrainment  means 
that  the  mixing  takes  place  within  the  interfacial  zone  whose  width  is 
very  small  compared  to  the  depths  of  the  two  layers.  Outside  the  inter- 
facial zone  the  layers  are  considered  homogeneous.  The  mean  value  of 
6/2  has  been  found  to  be  3 mm  with  a divergence  of  ±1  mm.  Therefore, 
after  comparing  the  values  of  6/2  to  the  depths  of  the  layers  we  can 
consider  that  the  previous  assumption  is  valid  (Table  2). 

The  theoretical  interface  is  defined  as  the  locus  of  points  at 
which  the  density  is  p ♦ Ap/2(p  + Ap).  The  distance  between  the  inter- 
face and  the  bottom,  h$,  is  defined  as  the  depth  of  the  lower  layer  and 
the  distance  between  the  interface  and  the  free  surface,  h,,  is  defined 
as  the  depth  of  the  upper  layer. 
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Figure  6.1  Definition  sketch  for  the  interfacial  zone 
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The  density  distribution  curve  is  approximately  linear  in  its 
central  portion  and  it  curves  at  its  two  ends  near  the  boundaries. 
Powell  [34]  gave  the  following  density  distribution  law. 


p£  - p 

Ap 


1 

2 


1 - tanh 


6/2 


(6.2.1) 


The  symbols  of  Equation  (6.2.1)  have  been  defined  in  Figure  6.1. 

Figure  6.2  gives  Powell's  [34]  experimental  data  and  Figure  6.3 
represents  some  experiments  of  the  present  study,  belonging  to  group  5, 
with  Ap/ (p  + Ap)  = 0.0225  approximately  equal  to  the  value  used  by 
Powell  [34]  (Table  1 and  Appendix).  In  both  figures  the  curves  given 
by  Equation  (6.2.1)  have  been  drawn  and  they  appear  to  agree  satisfac- 
torily with  the  experimental  data  for  0.1  ^ Ap£/(p  + Ap)  ^ 0.9. 

When  the  velocity,  , of  the  upper  layer  was  less  than  the 
critical  entrainment  velocity,  V , given  by  Keulegan  [22,  24],  the 
interface  was  smooth  (Figure  5.10).  For  V-j  larger  than  V , random 
interfacial  wavelets  were  formed  moving  downstream  with  a slowly 
increasing  wavelength  and  amplitude.  In  many  cases  the  wavelength  was 
of  the  order  of  some  centimeters  and  the  amplitude  of  some  millimeters. 
For  larger  velocities  V-j,  the  wave  crests  become  sharper  and,  after 
breaking,  the  two  fluids  were  mixed  at  the  interface  [Figure  5.12(a)]. 
The  latter  in  all  the  experiments  was  distinct  and  easily  identifiable. 

6.3  Velocity  Distribution 

The  time-averaged  horizontal  velocity  component,  u,  was  evalu- 
ated from  the  instantaneous  velocities  computed  from  the  hot  film 
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Figure  6.2  Nondimensional  density  distribution,  Powell's  data 


Source:  [34,  Figure  18] 
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Figure  6.3  Nondimensional  density  distribution  (data  from  present 
study) 
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Table  1.  Grouping  of  experiments 


Group 


Experiments 


Ap 

p + Ap_ 


7oo 


1 

2 

3 

4 

5 

6 

7 

8 


P1 0 ’ P9  ’ P7  ’ P6  ’ P5 ,06  ’ P4 5 N4  ’ P2 

i2,i3,i4,n1 ,n3,o3,o5,p1 

I-|  » 0 1 5 O2 5 Q-j 

L2,L6’L7’L9,V4 

,V-j  >Vg,Vg,S-]  ,S2,S4,U4,Ug , AA-j , AA2 

u1  ,u2,u3,v2,w4,w5,w6,y] ,y2,y4,y5,z1 , 

Z2’Z3’Y3 

wrw2,w3 

DA-j , DA2 , DA3 , DA4 , DAg , EA-j , EA2 , EA3 , EA4 , 

G A -|  ,GA2 


8.8 

10.4 
12.0 
17.8 

22.5 

35.5 

48.5 

82.0 
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anemometer  recordings.  Knowing  the  velocity  and  the  corresponding 
position,  the  dimensional  velocity  distribution  was  determined  for  each 
experiment.  Examples  are  given  in  Figures  A.l  to  A. 5 of  the  Appendix. 
Figure  6.4  is  a definition  sketch  with  the  following  parameters: 

a.  The  position  of  the  density  interface,  as  defined  in 
section  6.2  from  the  density  profiles 

b.  The  position  of  the  velocity  interface— This  interface 
is  defined  at  the  point  where  u = U/2  where  U is  the 
free  surface  velocity.  This  point  approximately  coin- 
cides with  the  inflection  point  at  the  interfacial 
zone  of  the  velocity  distribution  curve  (see  section 
2.3). 


c.  The  distance,  6u,  between  the  velocity  interface  and 
the  upper  limit  of  the  zone  of  shear  flow — This  limit 
is  considered  to  be  in  a position  where  the  velocity 

is  98  percent  of  U.  Above  this  limit,  the  flow  is  con- 
sidered frictionless  (Figures  A.l,  A. 2,  A. 4,  and  A. 5 
of  the  Appendix).  This  limit  corresponds  to  the  limit 
of  the  boundary  layer  for  homogeneous  fluids. 

d.  The  distance,  6u-j,  between  the  velocity  interface  and 
the  lower  inflection  point  of  the  velocity  distribution 
curve— A short  distance  below  that  point  the  velocity 
becomes  negative  (see  Figure  2.16  and  Figures  A.l  to 

A. 5 of  the  Appendix). 

e.  The  distance,  6u2>  given  as  6u^  = 5u  + 6u -j 
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Figure  6.4  Definition  sketch  for  velocity  profiles 
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f.  The  distance,  n,  between  the  two  interfaces 

g.  The  parameter,  e = 56u-j , defined  for  a better  repre- 
sentation of  the  experimental  data 

h.  The  distance,  yg,  between  the  bottom  and  the  velocity 
interface,  defined  as 

yQ  = hs  + n (6.3.1 ) 

For  a dimensionless  presentation,  the  velocity  has  been  non- 
dimensional  i zed  by  U and  the  vertical  distance  by  Su.  The  latter  has 
been  taken  with  reference  to  the  velocity  interface,  i.e.,  data  have 
plotted  versus  (y-yg)/6u.  After  some  trials,  the  dimensionless 
velocity  distributions  have  been  separated  into  two  distinct  groups. 

The  first  group  contains  the  experiments  of  Figure  6.5  and  the  second 
the  experiments  of  Figures  6.6  and  6.7.  The  separation  of  the  second 
group  into  two  subgroups  was  based  on  the  ratio  Ap/(p  + Ap).  Figure 

6.6  contains  the  experiments  with  Ap/(p  + Ap)  A 20  °/oo  and  Figure  6.7 
the  experiments  with  Ap/(p  + Ap)  > 20  °/oo.  The  curve  in  Figure  6.5 
differs  from  the  corresponding  curves  of  Figures  6.6  and  6.7;  Figure 
6.5  contains  mainly  experiments  with  small  velocities,  Vg,  upstream  of 
the  wedge.  The  dimensionless  velocity  distribution  of  Figures  6.6  and 

6.7  is  similar  to  the  one  given  by  Powell  [34]  (Figure  2.15).  Powell 
[34]  found  out  that  6u  has  no  significant  correlation  with  the  mean 
velocity  of  the  upper  layer  and  the  downstream  distance  or  with  the 
Reynolds  number  based  on  these  two  quantities.  Therefore,  the  flow  is 
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Figure  6.5  Dimensionless  velocity  distribution 
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Figure  6.6  Dimensionless  velocity  distribution  [Ap/(p  + Ap)  = 20  °/oo] 
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Figure  6.7  Dimensionless  velocity  distribution  [Ap/(p  + Ap)  > 20  °/oo] 
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believed  to  be  fully  developed.  This  is  not  valid  for  the  experiments 
of  Figure  6.5.  In  the  previous  figures  the  mean  curves  were  drawn  and 
transferred  onto  a semi  1 ogari thmi c paper.  Figure  6.5  corresponds  to 
Figure  6.8  and  Figures  6.6  and  6.7  correspond  to  Figure  6.9.  The  mean 
curves  of  Figures  6.6  and  6.7  coincide.  The  comparison  of  Figures  6.8 
and  6.9  suggests  that  in  fully  developed  flow, the  dimensionless  velocity 
distribution  is  logarithmic  and  it  consists  of  two  branches  intersecting 
each  other  at  approximately  the  point  where  (y  - yQ)/Su  = 0.15.  This 
value  is  very  near  the  limit  of  validity  of  the  wall  law  for  a solid 
boundary  for  homogeneous  fluids.  This  significant  observation  will  be 
discussed  later. 

The  subsequent  analysis  has  been  applied  to  the  experiments  of 
the  second  group,  where  the  flow  is  considered  fully  developed  (experi- 
ments of  Figures  6.6  and  6.7). 

Within  the  region  of  shear  flow  there  are  three  distinct  zones 
of  the  dimensional  velocity  distribution  curves:  (a)  the  zone  6u, 

(b)  the  zone  6u-j,  and  (c)  the  zone  in  the  immediate  neighborhood  of  the 
velocity  interface  (Figure  6.4).  To  study  the  flow  structure,  the  ver- 
tical distances  were  represented  in  each  zone  by  the  nondimensional 
relations  (y  - y0)/SuT,  (y  - yQ)/c  = (y  - yQ)/56u-| , and  (y  - yQ)/6u2, 
respectively.  The  velocities  were  nondimensionalized  by  the  free  sur- 
face velocity,  U. 

By  correlating  the  dimensionless  lengths  to  dimensionless 
velocities  a mathematical  expression  for  the  velocity  distribution  was 
sought  by  which  viscous  stresses  could  be  reliably  computed  at  any  point 
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Figure  6.8  Dimensionless  velocity  distribution  (data  from  Figure  6.5) 
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Figure  6.9  Dimensionless  velocity  distribution  (data  from  Figures  6.6  and  6.7) 
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of  the  cross  section.  The  general  laws  of  velocity  distribution  were 
derived  in  paragraph  6.5  in  terms  of  the  friction  velocity,  u*. 

First,  the  velocity  distribution  in  the  immediate  neighborhood 
of  the  velocity  interface  was  studied.  The  plots  in  Figures  6.10  and 
6.11  clearly  suggest  that  this  distribution  is  linear;  the  correlation 
coefficient,  R,  is  equal  to  98  percent  for  the  data  of  Figure  6.10  and 
98.6  percent  for  the  data  of  Figure  6.11.  The  linear  distribution  zone 
has  a width  between  2 mm  and  6 mm  with  an  average  of  4 mm.  The  linear 
velocity  distribution  can  be  satisfactorily  approximated  by  the  equa- 
tion 


+ 0.5 


(6.3.2) 


where  the  symbols  are  given  in  Figure  6.4. 

Next,  the  velocity  distribution  in  the  zones  of  width  6u  and  6 u-j 
have  been  simultaneously  studied.  A wide  scattering  of  data  points 
was  first  observed  in  each  one  of  the  two  original  groups  for 
Ap/(p  + Ap)  = 20  °/oo  and  Ap/(p  + Ap)  > 20  °/oo,  respectively. 

This  scattering  is  indicative  of  some  dependence  of  dimension- 
less velocities  on  Ap/ (p  + Ap)  in  the  aforementioned  zones,  in  spite  of 
the  fact  that  in  the  zone  of  linear  velocity  distribution  the  dimension- 
less velocity  profile  appears  to  be  practically  independent  of 
Ap/(p  + Ap).  To  further  elucidate  that  dependence,  it  has  been  decided 
to  separate  the  experiments  into  eight  groups  according  to  Ap/(p  + Ap). 

Table  1 shows  the  grouping  of  experiments  with  each  group  iden- 
tified by  its  mean  value  of  Ap/(p  + Ap). 
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Figure  6.10  Dimensionless  velocity  distribution  in  the  immediate  neighborhood  of  velocity  inter- 
face for  Ap/(p  + Ap)  = 20  °/oo 
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Figure  6.11  Dimensionless  velocity  distribution  in  the  immediate  neighborhood  of  velocity  inter- 
face, for  Ap/(p  + Ap)  > 20  °/oo 
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The  graphs  of  Figures  6.12  to  6.19  give  the  dimensionless 
velocity  distribution  outside  of  the  linear  zone  for  the  eight  groups 
of  experiments  (Table  1).  These  graphs  suggest  a logarithmic  velocity 
distribution  in  the  zones  5u  and  Su-j  with  values  of  the  correlation 
coefficient  R ranging  from  61  percent  to  98.1  percent.  In  the  zone  6u 
the  semi  logarithmic  plot  consists  of  two  linear  segments  which  inter- 
sect approximately  at  the  point  (y  - y0)/6u  = 0.15.  It  is  noteworthy 
that  in  boundary  layers  with  a width  6,  and  with  fixed  boundaries,  the 
point  for  which  y/S  = 0.15  is  the  limit  of  validity  of  the  wall  law. 

We  can  define  tentatively,  therefore,  the  value  of  (y  - yg)/6u  = 0.15 
as  the  limit  of  validity  of  the  "law  of  the  interface"  which  cor- 
responds to  the  wall  law  for  solid  boundaries.  In  the  zone  5u-|  pre- 
liminary analysis  also  leads  to  a similar  logarithmic  law.  These 
previous  laws  were  subsequently  used  for  the  calculation  of  viscous 
stresses . 


6.4  Laws  of  Variation  of  Interfacial 
and  Bed  Friction  Coefficients 

The  local  friction  coefficients  f ^ and  fg  at  the  interface  and 
at  the  bed,  respectively , are  defined  by  the  relations 


pv? 


and 


8t, 


PV 


(6.4.1) 


where  t.  and  Xg  are  the  local  shear  stresses  at  the  interface  and  at 
the  bottom,  respecti vely . The  shear  stresses  in  a turbulent  flow  are 
given  by  the  equation 
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Figure  6.12  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (group  1) 
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Figure  6.13  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (group  2) 
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Figure  6.14  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (group  3) 
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Figure  6.15  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (group  4) 
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Figure  6.16  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (group  5) 
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Figure  6.17  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (group  6) 
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Figure  6.18  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (group  7) 
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(6.4.2) 


where  x^  are  the  viscous  stresses  and  x^.  are  the  Reynolds  stresses. 

The  viscous  stresses  were  obtained  from  the  velocity  distribution  curve 
while  the  Reynolds  stresses  are  given  by  the  correlation  of  the  turbu- 
lent velocity  fluctuations  u'  and  v1. 

For  the  evaluation  of  viscous  stresses,  the  mean  curves  of 
each  group  of  experiments  in  Figures  6.12  to  6.19  were  used.  However, 
a more  careful  study  of  the  graphs,  particularly  in  the  region 
0 < (y  - yg)/6u  ^ 0.15,  led  to  the  observation  that  the  slope  of  the 
straight  line  differs  from  experiment  to  experiment  as  well  as  from 
the  mean  straight  line  of  the  experimental  group.  Therefore,  the 
straight  line  of  each  particular  experimental  velocity  distribution  was 
used  for  the  computation  of  the  viscous  stress  instead  of  the  average 
line.  These  curves  are  presented  in  Figures  A. 6 to  A. 12  of  the  Appen- 
dix, as  an  example  for  a part  of  the  experiments  of  group  5 (Table  1). 
The  velocity  distribution  law  in  zone  6u-j , in  general,  is  expressed  by 
two  straight  lines  intersected  at  point  (y  - Yg)/e  3 0.1.  That  point 
has  been  defined  as  the  limit  of  the  law  of  the  interface  for  zone  6u-| . 


The  general  form  of  the  laws  of  the  interface  are 


U 


u 


(6.4.3) 


for  the  6u  zone  and 


(6.4.4) 
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for  the  6u-|  zone,  where  a2»  33*  b^,  are  functions  of  the  flow  and 
salinity  parameters  and  e is  taken  equal  to  56u^ . In  the  immediate 
neighborhood  of  the  velocity  interface  the  1 inear  law  [Equation  (6.3.2)] 
has  been  used. 

Figure  6.20  constitutes  a definition  sketch  for  the  velocity 
distribution  on  a semi  logarithmic  paper. 

In  the  zone  6u-j,  points  1 and  2 define  the  limits  of  the  region 
of  the  law  of  the  interface.  In  the  same  way,  points  3 and  4 define 
similar  limits  in  the  zone  6u.  Between  points  2 and  3 the  linear  law 
is  valid.  The  viscous  stresses  at  any  position  y of  the  cross  section 
can  be  expressed  by  the  following  relationships: 


a3U 

TV  * 11  (y  - y0) 
for  the  6u  zone 

a2U 

Tv  = y (yQ  - y) 

for  the  6u-|  zone  and 

a-,11 


(6.4.5) 


(6.4.6) 


(6.4.7) 


in  the  immediate  neighborhood  of  the  velocity  interface,  where  a-j  is 
the  slope  of  straight  lines  in  Figures  6.10  and  6.11  [in  Equation 
(6.3.2)  a-j  = 5],  From  Equations  (6.4.5)  to  (6.4.7)  and  from  Figure  6.20, 
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it  follows  that  the  above  stresses  reach  their  maximum  value  at  one 
of  the  positions  0,  1,  2,  3,  or  4.  In  particular,  at  position  4 the 
viscous  stresses  are  lower  than  at  any  other  position  of  line  lg 
[Equation  (6.4.5)],  because  the  slope  of  1^  is  less  than  the  slope  of 
lg.  Moreover,  the  value  of  y - yg  in  1^  is  larger  than  the  correspond- 
ing value  of  y - yg  at  any  point  of  lg. 

For  the  experiments  of  groups  1,  2,  3,  and  4 with  Ap/(p  + Ap) 

= 20  °/oo,  the  corresponding  positions  of  points  0,  1,  2,  and  3 
(Figure  6.20)  are  (yg  - y)/e  = 0.0,  (yg  - y)/e  = 0.12  [and  in  some 
experiments  (y0-y)/e  = 0.077],  (yQ  - y)/e  = 0.047  and  (y  - yQ)6u 
= 0.037,  respectively.  For  the  experiments  of  groups  5,  6,  7,  and  8 
with  Ap/(p  + Ap)  > 20  °/oo  the  corresponding  positions  of  points  0,  1, 

2,  and  3 are  (yg-  y)/e  = 0.0,  (yg  - y)/e  = 0.12  [and  in  some  experiments 
(y0  - y)/e  = 0.08],  (yQ  - y)/e  = 0.048  [and  in  some  experiments 
(yQ  - y)/e  = 0.028],  and  (y  - y0)/6u  = 0.023  [and  in  some  experiments 
(y  ~ Yg ) /Su  = 0.04].  In  some  of  the  points  0,  1,  2,  3 correspond  two 
positions.  For  example,  at  point  1 of  the  groups  with  Ap/(p  + Ap) 

5 20  °/oo,  correspond  the  positions  (yQ  - y)/e  = 0.12  and  0.077;  that 
means  that  in  some  experiments  point  1 corresponds  to  (yg  - y)/e  = 0.12 
and  in  the  remaining  experiments  point  1 corresponds  to  (yg  - y)/e 
= 0.077.  The  maximum  Reynolds  stresses,  x.^.,  have  been  defined  in  the 
immediate  vicinity  of  points  0,  1,  2,  and  3.  Next,  the  total  stresses 

Ti0  = ^Tiv^0  + ^Tit^0  P°int  0 have  been  calculated.  Here  (x..  )g 
are  the  viscous  and  (x^g  are  the  Reynolds  stresses  at  (yg  - y)/e  = 0.0. 
Next,  the  Reynolds  stresses  and  the  viscous  stresses  x..  at  points 
0,  1,  2,  and  3 have  been  nondimensional ised  by  x.g. 


The  results  are 
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given  in  Figures  6.21  and  6.22.  In  these  figures  the  left  vertical 

axis  represents  the  relative  distance  of  the  position  from  the  velocity 

interface  with  respect  to  6u  or  with  respect  to  e = 56u^ . The  right 

vertical  axis  represents  the  absolute  value  of  the  above  distances  in 

millimeters  and  the  horizontal  axis  represents  the  relative  Reynolds 

stresses  x-*  or  viscous  stresses  x-  with  respect  to  the  total  stresses 
it  i v 

x.j  at  (yQ  - y)/e  = 0.0  (that  is,  x.g).  From  Figures  6.21  and  6.22  and 
from  the  governing  equations  of  the  velocity  distribution,  the  total 
stress  x j = x.jv  + for  each  experiment  has  been  calculated  at  points 
0,  1,  2,  and  3;  the  highest  of  these  values  was  taken  as  the  local 
interfacial  stress  x-.  The  local  interfacial  coefficient  f^  was  calcu- 
lated from  the  relationship 


fi  ■ 


^Tiv  + Tit^max 


(6.4.8) 


where  (x.  + T-j-t-)rTiax  is  the  maximum  value  of  the  total  shear  stress  in 
the  cross  section  of  measurements,  except  in  the  vicinity  of  the  bottom. 

Next,  on  the  basis  of  the  theoretical  analysis  developed  in 
section  3.4,  the  mean  interfacial  friction  coefficient  over  the  length 
of  salt  water  intrusion  has  been  calculated  by  means  of  Equation  (3.4.16) 
The  mean  value  of  the  depth  h$  has  been  numerically  calculated  and  that 
of  h-|  has  been  calculated  from  the  relation  h-j  = hg  - h . The  depth  of 
the  lower  layer  at  the  mouth  of  the  channel,  hsg,  is  the  mean  value  of 
the  depths  in  three  successive  cross  sections  in  a distance  between 


them  of  about  2 m. 
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Figure  6.21 


Dimensionless  values  of  Reynolds  and  viscous  stresses 
for  Ap/(p  + Ap)  ^ 20  °/oo 


140 


Figure  6.22  Dimensionless  values  of  Reynolds  and  viscous  stresses 
for  Ap/(p  + Ap)  > 20  o/oo 
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Table  2 gives  the  values  of  the  parameters  entering  in  Equation 
(3.4.16)  as  well  as  the  dimensionless  values  of  Su  and  6u-j . In  column 
26,  the  values  of  f..  based  on  Equation  (6.4.8)  are  listed,  while  those 
values  estimated  from  Equation  (3.4.16)  are  listed  in  column  27. 

A comparison  of  f.  suggests  that  these  two  values  of  f^  based 
on  direct  measurements  and  on  the  integration  of  the  equations  of  motion, 
respectively,  are  very  close  to  each  other;  therefore,  f..  can  be 
reasonably  approximated  from  Equation  (3.4.16). 

The  values  of  f.  calculated  on  the  basis  of  Schijf  and 
Schonfeld's  model  [37]  [Equation  (3.2.17)]  are  listed  in  column  29. 

Since  for  the  reasons  explained  in  Chapter  IV,  the  length  of  salt  water 
intrusion,  Lq,  in  the  channel  was  smaller  than  the  full  length,  L,  of 
the  arrested  saline  wedge  (Figure  3.4),  the  calculation  of  L had  to  be 
first  computed  in  the  following  way. 

For  each  experiment  of  given  densimetric  Froude  number,  Fr^, 
the  critical  depth  of  the  fresh  water  layer,  h.|  , was  computed  from 
Equation  (3.2.15).  The  critical  salt  water  depth  h$c  = hg  - h-jc  has 
next  been  used  to  nondimensionalize  the  wedge  profile,  as  hs/h  . 

From  Equation  (3.3.14)  or  from  the  graph  of  Figure  3.2,  the  total 
intrusion  length,  L,  was  estimated  for  several  values  of  hs/h$c  and 
distances  x^  from  the  channel  entrance.  The  mean  value  of  L,  thus 
determined,  with  the  exception  of  certain  data  that  appeared  to  be  in 
gross  error,  was  taken  to  represent  the  true  total  intrusion  length,  L, 
of  the  wedge. 

The  number  of  the  rejected  values  was  small  in  comparison  to 
their  total  number.  Finally,  f.  has  been  estimated  by  Schif  and 


Table  2.  Values  of  parameters  introduced  in  Equation  (3.4.16) 


No 

Exp 

Syr: 

Ad 

h 

hsl 

A 

v! 

L 

0 

L 

p+Ac 

0 

nl 

n 

V1 

1 

h10 

o/oo 

cm 

cm 

cm 

cm 

sec 

cm 

cm 

cm 

cm 

sec 

m 

m 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

ii 

12 

13 

14 

15 

16 

17 

1 

P10 

O 

7.9 

14,01 

11,76 

0.77 

2.21 

1.60 

1.41 

12.35 

11.34 

2.11 

14.16 

42.6 

1,033 

0.325 

2 

P9 

• 

8.2 

14,06 

11.81 

0.86 

2.47 

1,60 

1.41 

12.26 

11,33 

2,73 

14.02 

37.5 

1.041 

0,266 

3 

P7 

3 

8.7 

34,17 

11.74 

0.81 

3.01 

1,73 

1.41 

12.08 

11,07 

2.92 

14,66 

32,4 

1.052 

0.218 

4 

P6 

C 

9.0 

14.22 

11.56 

0.65 

3.35 

1.91 

1.39 

11.96 

10.95 

3.23 

15.27 

29.3 

1.027 

0,193 

5 

P5 

e 

9.3 

14.28 

11.51 

0.75 

3.61 

1.97 

1.41 

11.94 

10.75 

3.48 

15.46 

27.6 

1.045 

C.242 

6 

°6 

e 

9.4 

11.84 

9.62 

0.45 

2.76 

1.57 

1 .41 

10.00 

9.09 

2.65 

12.21 

23.0 

1.030 

0.240 

7 

P4 

6 

9.7 

14.34 

11.46 

0.76 

3,76 

2.08 

1.39 

11.75 

10.47 

3.66 

17.04 

27.6 

1.047 

0.190 

8 

N4 

o- 

9.5 

16.85 

14.35 

0.26 

2.97 

1,80 

1.39 

14.93 

13.82 

2.88 

15,63 

52.3 

1.014 

C.260 

9 

P2 

9 

9.8 

14.33 

12.43 

0.56 

3.61 

1.35 

1.41 

12.65 

11.86 

3.55 

9 19 

22.6 

1.017 

0.173 

10 

l2 

A 

10,8 

10.82 

9.25 

0.55 

3.49 

1.12 

1.40 

9.15 

8.44 

3,52 

8,46 

16. C 

1 056 

0.350 

11 

h 

A 

10,6 

14.13 

11,48 

0.87 

3.82 

1.90 

1.40 

11.  S3 

10.61 

3.67 

15.82 

27.7 

1.050 

0.333 

12 

U 

4 

10.5 

14.39 

12.34 

0.02 

3.88 

1.40 

1.46 

12.62 

11,27 

8.80 

11.19 

23.5 

1.015 

C.I76 

13 

Ni 

L 

10.3 

17.05 

14.15 

0,35 

3.48 

2,05 

1.42 

14,94 

13.62 

3.30 

18.30 

52.2 

1.020 

0.219 

14 

N3 

4 

10,1 

16.94 

14.44 

0.51 

3.11 

1.80 

1.39 

15.10 

13.84 

2.97 

15,66 

49.7 

1.017 

0.246 

15 

°3 

4 

10.7 

11.97 

9.70 

0.44 

3.19 

1.62 

1.40 

10.09 

9.30 

3.06 

10.83 

20.7 

1.036 

0.181 

16 

°5 

A 

10.0 

11.91 

9.53 

0.46 

2.97 

1.68 

1.42 

9,92 

8.83 

2.86 

13.95 

24.3 

1.023 

0.198 

17 

P1 

t 

10,0 

14.36 

12.11 

0.40 

4.29 

1,60 

1.41 

12.36 

11.34 

4.24 

11.10 

2 i . 5 

1.026 

0.165 

18 

h 

E 

12.0 

10.60 

8.20 

0.44 

3.68 

1.7C 

1.41 

8.84 

7.75 

3.42 

12.96 

19.1 

1.013 

0.261 

19 

°1 

□ 

11.1 

12.12 

10.36 

0.33 

3,67 

1,26 

1.40 

10.06 

9.44 

3.76 

8.50 

14.9 

1.021 

0.139 

20 

°2 

■ 

11.1 

12.01 

9.52 

0,46 

3,58 

1.79 

1.39 

9.79 

8.98 

3.48 

12,62 

19.1 

1.041 

0.134 

21 

Q1 

a 

12.5 

11.93 

9.16 

0.16 

4.31 

1.97 

1.41 

9.77 

8.73 

4.04 

18.92 

15,8 

1.018 

0.192 

22 

<>3 

b 

12.6 

11.81 

9.21 

0.35 

3,70 

1.85 

1.41 

9.58 

8.79 

3.14 

12.80 

22.6 

1.032 

0.213 

23 

^4 

9 

11.6 

11.76 

9.66 

0.80 

3,41 

1.50 

1.40 

10.18 

9.19 

3.23 

11.47 

23.3 

1 .055 

0.323 

24 

% 

» 

11.3 

11.74 

9.48 

0.57 

3.30 

1.61 

1.40 

10.13 

9.14 

3.12 

11.86 

24.8 

1.036 

0.201 

25 

% 

11.3 

11.65 

9.05 

0.18 

2.99 

1,85 

1.41 

9.58 

8.72 

2. S3 

12.97 

24.5 

1 .006 

0.222 

26 

h 

c 

12.1 

33.73 

30,75 

0.29 

4.37 

2.12 

1.41 

31.52 

30.36 

4.26 

17.85 

158.3 

1.007 

0.164 

27 

ca3 

m 

12.8 

29.50 

26.43 

0.79 

6,54 

2,07 

1.46 

26.34 

24.97 

6.56 

12.44 

35,2 

1.020 

0.221 

28 

L2 

h 

16.6 

6.54 

5.24 

0.33 

3,96 

0.95 

1.37 

5,26 

4.71 

3.94 

7.62 

9.2 

1.056 

0.250 

29 

L6 

— 1 

17.3 

11.22 

8.94 

0.31 

4.57 

1.53 

1.40 

9.28 

8.35 

4.40 

12.91 

19.8 

1.022 

0.190 

30 

L7 

+ 

16.6 

11.32 

9.37 

0.44 

4.41 

1.40 

1.39 

9.49 

8,62 

4,35 

9. "72 

15.7 

1 .035 

0.213 

31 

L9 

X 

16.8 

11,39 

9.12 

0.32 

4.54 

1.57 

1.45 

9.46 

8.22 

4.42 

1-1 . 28 

16.2 

1.045 

0,198 

32 

V4 

•* 

18.9 

29.64 

27.81 

0.40 

3.68 

1,28 

1.43 

28.23 

27.62 

3.62 

11.15 

246,7 

1.011 

0,203 

33 

R 

o 

20.9 

11.66 

9.45 

0.37 

4,26 

1.60 

1.41 

10,61 

9.54 

4.04 

10.63 

27.4 

1.021 

0,237 

34 

S5 

• 

20.9 

19.01 

16.53 

0.39 

3.76 

1.78 

1.39 

17.08 

16.04 

3.64 

16.23 

71.1 

1,022 

0.261 

35 

T1 

3 

20,7 

29.76 

27.14 

0.31 

4.86 

1.82 

1.44 

27.91 

26.81 

4.69 

16.78 

167.6 

1.011 

0.119 

36 

V1 

C 

22.5 

20,42 

18.59 

0.50 

3.69 

1.33 

1.38 

18.96 

18.29 

3.55 

13.58 

113.7 

1.019 

0.339 

37 

V5 

w 

19.8 

34.57 

32.45 

0.12 

4.17 

1.52 

1.40 

32.49 

31.70 

4.16 

15.34 

203.8 

1,007 

2.123 

38 

V6  J 

© 

20.9 

34.18 

31.61 

0.18 

3.60 

1 .82 

1.41 

32.19 

31.10 

3.53 

18.73 

128.0 

1.005 

0.100 
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Table  2— Continued 


£ 

T 

A 

Fr 

0 

5u/h 

ns  . 

!r 

Rerr 

2 f. 

1 

f . 

T 

f . 
1 

f . 

l 

7. 

2 

r 

£ 

' 0 

Ho. 

°C 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

0.24? 

27. 9( 

0,211 

0.17? 

0.58/ 

0.765 

3.909 

12.2 

0,0254 

0,0239 

0.0340 

0,0334 

1.40 

0.025C 

0,0100 

1 

0.22  * 

27.9: 

0.221 

0,195 

0.682 

0.539 

3.721 

16.0 

0.0200 

0.0201 

0.0275 

0,0310 

1.40 

0.020C 

0.0080 

2 

0.233 

27.8C 

0.229 

0.226 

0.675 

0.488 

3.402 

25.4 

0.0184 

0.0155 

0,0231 

0,0251 

1 .40 

0,0183 

0.0072 

3 

0.343 

27 ,8C 

0.236 

0.243 

0.653 

0.473 

3.248 

31.6 

0.0147 

0,0135 

0.0199 

0,0231 

1.40 

0.015C 

0.0050 

L 

0.318 

27. 7C 

0.236 

0.255 

0.692 

0.513 

3.148 

37.4 

0.0122 

0.0129 

0.0199 

0.0216 

1 .40 

0,01 30 

0.0052 

c 

0,333 

27.05 

0.281 

0.214 

0,634 

0,543 

3.519 

18.8 

0.0245 

0,0213 

0.0298 

0.0312 

1.40 

0,0240 

0.0095 

6 

0.239 

27.60 

0.238 

0.257 

0,621 

0,409 

3.132 

40.7 

0.0130 

0,0134 

0,0196 

0.0212 

1 .40 

0,0125 

0.0050 

7 

0.301 

27.90 

0.189 

0.203 

0.519 

0.516 

3.634 

28.6 

0,0167 

0.0172 

0,0248 

0.0241 

1.40 

0,01 65 

0,0060 

8 

0.277 

27.50 

0.225 

0,267 

0.742 

0.542 

3.053 

45.0 

0.0125 

0.0133 

0.0198 

0,0224 

1 .40 

0,0117 

0.0047 

0 

0,317 

26.70 

0.277 

0.278 

0,673 

0.455 

2.970 

37,1 

0,0179 

0.0137 

0.0209 

0.0219 

1 ,40 

0.0155 

0.0062 

10 

0.246 

26.80 

0.230 

0.256 

0,616 

0,568 

3.140 

44.3 

0.0132 

0.0132 

0.0203 

0,021 1 

1 .40 

0,0135 

0.0055 

i i 

0.501 

27.10 

0.219 

0.259 

0,744 

0.322 

3.116 

46.4 

0,0138 

0,0160 

0.0238 

0,0226 

1 .40 

0.0135 

0.0054 

12 

0.360 

26,80 

0.187 

0.220 

0,654 

0.379 

3.460 

39.5 

0.0154 

0.0143 

0,0217 

0.0200 

1 .40 

0.0148 

0.0059 

13 

0,294 

26.90 

0.188 

0.205 

0.519 

0.658 

3,612 

31.9 

0.0171 

0.0178 

0.0262 

0.0249 

1 .40 

0,0170 

0.0068 

14 

0.245 

26,80 

0.263 

0.230 

0.659 

0,578 

3.365 

27,3 

0.0208 

0.0176 

0.0269 

0,0316 

1 .40 

0.01 88 

0.0075 

15 

0.439 

27,05 

0.270 

0.220 

0.747 

0.462 

3.460 

22.2 

0,0197 

0.0195 

0.0285 

0.0301 

1 ,40 

0,0220 

0.0088 

16 

0.349 

27.30 

0.225 

0,308 

0,693 

0.446 

2.765 

64.3 

0,0115 

0.0101 

0,0148 

0.0163 

1 ,40 

0.0117 

0,0047 

17 

0,353 

26.60 

0.301 

0.255 

0,591 

0,409 

3.148 

33.4 

0.0170 

0.0152 

0,0215 

0.0232 

1.40 

0,0172 

0,0069 

1 8 

0,303 

26.10 

0,267 

0.272 

0,610 

0.364 

3,015 

42.7 

0,0181 

0.0163 

0,0232 

0.0285 

1 .40 

0,0158 

0,0063 

19 

0.322 

26.50 

0,277 

0,248 

0,608 

0,369 

3,206 

34.5 

0,0169 

0.0149 

0,0214 

0.0282 

1 ,40 

0,0172 

C , 0069 

20 

0,469 

25.20 

0.270 

0.274 

0,660 

0,352 

3.000 

49.8 

0.0141 

0.0168 

0,0234 

0,0259 

1 .40 

0,0150 

0,0060 

21 

0.423 

27.10 

0.272 

0.211 

0,669 

0,552 

3.550 

28.0 

0.0178 

0.0218 

0.0305 

0,0355 

1,40 

0.0195 

0,0078 

22 

0.245 

27.25 

0,271 

0.242 

0.667 

0.475 

3.257 

34.2 

0.0172 

0.0163 

0,0233 

0.0242 

1.40 

0.0172 

0,0069 

23 

0,301 

27.30 

0.264 

0.236 

0.670 

0,484 

3.310 

30.8 

0,0178 

0,0168 

0,0241 

0.0243 

1 .40 

0,0187 

0,0075 

24 

0.469 

27.25 

0.277 

0.205 

0.637 

0,454 

3,612 

21,8 

0.0233 

0.0232 

0.0322 

0.0348 

1 .40 

0.0235 

0.0094 

25 

0,389 

27.35 

0,098 

0.199 

0.185 

0.407 

3.677 

111.5 

0.0124 

3.0106 

0,0165 

0,0168 

1.41 

0.0117 

0,0047 

26 

0.283 

27.10 

0.116 

0.304 

0.198 

0.475 

2.791 

267.5 

0.0091 

3.0103 

0,0169 

3,0180 

1 ,48 

0.0088 

0,0041 

27 

0.459 

24.80 

0,490 

0.307 

0.676 

0.461 

2,765 

26.4 

0,0233 

3,0192 

3.0287 

3,0289 

1 ,40 

0,0230 

0,0092 

28 

0.368 

25.80 

0.287 

0,264 

0.629 

0,402 

3,076 

68.7 

0.0151 

3.0129 

3.0185 

3.0215 

1 .40 

0,0156 

3 ,0062 

29 

0.447 

26.00 

0.280 

0.267 

0.670 

0,678 

3,052 

69.1 

0.0150 

3,0170 

3.0233 

3.0286 

1 .40 

0,0153 

3.0061 

30 

0.452 

26.50 

0,271 

0,268 

0.462 

3.477 

3,045 

70.8 

0.0158 

,0172 

3,0234 

3.0235 

1 ,40 

0,0152 

3,0061 

31 

0,320 

27.35 

0,101 

0,147 

0.318 

3.496 

.,346 

85,3 

0.0150 

.0156 

3.0233 

3,0192 

1 .40 

0,0140 

3,0056 

32 

0,309 

27.15 

0.265 

0.238 

3.685 

3,923 

5,293 

78,3 

0.0184  C 

.0156 

’,0220 

,0213 

1 .40 

0,0174 

3,0070 

33 

0.314 

27.50 

0.163 

3.166 

3.481 

3,445 

,076 

68.8 

0.0173  C 

.0231 

3,0328 

.0324 

1 .40 

0,0183 

3,0073 

34 

0.363 

27.30 

0,109 

3.179 

3.378 

3.368 

,909 

59.9 

0.0143  C 

,0126 

.0194  C 

.0180 

1 .42 

0.0138 

3,0058 

35 

0,238 

26.90 

0,153 

3,155 

3.425 

3,616  4 

.228 

72.3 

3,0182  0 

,0172  C 

.0244  C 

.0255 

1.40 

0.0180 

. 007  2 

36 

0.456 

27.30 

0.093 

3.151 

3.223 

.332  4 

.286  1 

23.0 

3,0164  0 

.0168  C 

.0250  C 

,0255 

1 .40 

0.0146 

.0059 

37 

0,390 

’7.10 

0,096 

3,125 

'.314  C 

.417  4 

.706 

82,3 

3.0188  0 

,0225  C 

,0343  0 

.0343 

1,40 

0,0170 

.0068 

38 
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2 

3 

4 

C 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

39 

S1 

A 

24.  C 

9.82 

7.25 

0.09 

5. 2( 

1.82 

1.4 

7.96 

7.25 

4.74 

11.76 

18.5 

1.021 

0.167 

40 

S2 

i 

24.1 

9.62 

7.59 

0,46 

4.21 

1.45 

1,42 

7.84 

7.19 

4.11 

13.11 

23.1 

1,03c 

0.220 

41 

S4 

A 

23.8 

19,29 

16.69 

0.49 

4.81 

1.85 

1.4] 

17,12 

16.16 

4.75 

16.60 

59,8 

1.01? 

0.192 

42 

U4 

L 

25.1 

10.13 

7.78 

0.14 

3.71 

1.7C 

1.38 

8.42 

7.41 

3.49 

16.70 

29.3 

1.017 

0.22S 

43 

U5 

4 

23.5 

9,95 

7.93 

0.42 

3.01 

1.42 

1,42 

8.43 

7.60 

2.89 

15.30 

33,7 

1 .027 

0,246 

44 

AAj 

+ 

23,9 

30,49 

27.32 

0.28 

7.96 

2.02 

1.57 

27.20 

26.26 

7.99 

13.37 

53.5 

1.005 

0.170 

45 

aa2 

A 

23.4 

30.73 

27.71 

0.24 

7.95 

1.92 

1.57 

27.74 

26.63 

7.94 

12,12 

53.5 

1.005 

0.169 

4G 

U1 

□ 

33,2 

7.06 

5,54 

0.14 

4.25 

1.12 

1.36 

5.82 

5,29 

4.04 

9.69 

16.5 

1.012 

0.156 

47 

U2 

■ 

33.8 

6.91 

5.14 

0,29 

3.68 

1.27 

1.39 

5.56 

4.84 

3.41 

13,47 

21.5 

1.039 

0.216 

43 

U3 

a 

30,7 

10.19 

8.47 

0.32 

4.10 

1.27 

1.35 

8.83 

8.21 

3.94 

11.27 

30.4 

1.024 

0,210 

49 

V2 

c 

32.0 

21.45 

19,35 

0.21 

5.80 

1.45 

1.45 

18,88 

18,35 

5.94 

13.00 

55.7 

1.010 

0.162 

50 

W4 

E 

40.2 

10.43 

8.47 

0.30 

4.08 

1.41 

1.39 

8.83 

8.21 

3.91 

16,26 

39.1 

1.020 

0,201 

5! 

W5 

a 

37.1 

10.30 

8.45 

0,26 

3.43 

1.30 

1.42 

8.89 

8,18 

3.26 

15.95 

44.4 

1.028 

0.209 

52 

U6 

CD 

30.7 

10.11 

8.26 

0.33 

2.86 

1,30 

1.42 

8.81 

8.79 

2.69 

16.82 

48.6 

1.021 

C , 330 

53 

Y1 

B 

32.5 

20.31 

17.76 

0.52 

7.25 

1,70 

1.50 

13,16 

17.18 

7.09 

13.7]“ 

47.1 

1.020 

0,224 

54 

Y2 

6 

34.0 

21.12 

18,57 

0.02 

8.50 

1.65 

1.55 

18.57 

16.74 

8.35 

15.13 

37.2 

1.017 

0,179 

55 

Y4 

9 

45.7 

30,90 

27.90 

0.48 

7.17 

1.95 

1.54 

28.  OS 

26.90 

7.12 

19.08 

107.6 

1.016 

0.196 

56 

Y5 

0 

31,9 

30.24 

27,49 

0.14 

6.25 

1.82 

1.51 

27.68 

26.87 

6,21 

17.19 

118,3 

l.OOS 

0.147 

57 

Z1 

-a 

32.8 

6,83 

5.29 

0.28 

5.18 

1.09 

1.41 

5.43 

4.65 

5.04 

9.26 

10.6 

1.030 

0.204 

58 

Z2 

■ 

32,6 

6.99 

6.57 

0,10 

5,51 

1.02 

1.39 

5.59 

4.96 

5,49 

10,46 

8.7 

1.019 

0.149 

59 

Z3 

■ 

35,5 

9.84 

7.62 

0.11 

5.17 

1.57 

1.41 

7.90 

7,04 

4.93 

14.98 

22.5 

1.022 

0.165 

60 

Y3 

m- 

34.4 

31.85 

28,30 

0.65 

8.97 

1,90 

1.58 

28.89 

27,32 

8.79 

16.93 

80.6 

1.012 

0.229 

61 

W1 

0 

49.5 

6,73 

5.15 

0,07 

5.44 

1.13 

1.40 

5.30 

4.75 

5.28 

8.86 

12.0 

1,070 

0.149 

62 

W2 

♦ 

49.0 

6.68 

5.83 

0,00 

5.00 

1,30 

1.42 

5.22 

4.65 

4.63 

10,96 

15.6 

1,014 

0,196 

63 

W3 

<t> 

47.5 

6.55 

4.97 

0.20 

4.08 

1.13 

1.40 

5.24 

4.65 

3,88 

11.64 

17.3 

1.024 

0.263 

64 

DAj 

1 

93.0 

7.24 

5,44 

0.25 

6.11 

1.25 

1.44 

6.06 

5.17 

5.48 

16.29 

29.7 

1.029 

0.184 

65 

da2 

93.5 

7.10 

5,39 

0.16 

6.24 

1.21 

1.41 

5. 68 

5,00 

5.93 

14.00 

20.6 

1.028 

0.140 

66 

da3 

L 

90.0 

7.26 

5.43 

0.06 

6.49 

1.28 

1.43 

5,98 

5,26 

5.89 

13.71 

23.5 

1.016 

0.114 

67 

da4 

J 

85.8 

9.08 

7.20 

0.15 

5.73 

1,33 

1.42 

7.80 

6. 98 

5.29 

15.93 

37.7 

1.023 

0,137 

68 

DAs 

-L 

82.4 

9.01 

6.94 

0.20 

6.24 

1.42 

,.46 

7.64 

6.78 

5.67 

15.86 

33.2 

1.021 

0.169 

69 

£A1 

H 

80.1 

20.95 

18.68 

0.23  1 

0.9! 

1.27 

.79 

18.63 

7.52 

0,94 

14,06 

47.7 

1.015 

0.166 

70 

EA2 

H 

78.3 

20.42 

18.62 

0.35 

9,30 

1,05 

.71 

18.79 

8.09 

9.22 

10.84 

64.9 

1.015 

0.211 

71 

E/U 

0 

T 

88.2 

30.16 

27.79 

0,38  1 

0.86 

1.32 

.80 

17,70  2 

6.66 

0,90 

17,03 

120.5 

1.011 

0.157 

72 

EA4 

+ 

82.3 

29.14 

27,02 

1.26 

7.73 

1.32 

.61 

’7 ,33  2 

6,61 

7.65 

16.90 

220.8 

1 .013 

0,133 

i 6 

OA] 

X 

70,4 

20,86 

8.68 

5.44 

8.41 

.33 

.64 

8.90  1 

7.99 

8,31 

15.64 

78.8 

1.016 

0.185 

74 

GA2 

-» 

72.6 

21.55 

8.93 

1.24 

8.47 

.62  1 

.62 

9.69  1 

8.49 

8.14 

15.40- 

80.4 

1.015 

0.164 

/ l r)t  \ 0 , 25 

( b ^ 2h,  ) 

A = ( r — 
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.18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

0.477 

25.90 

0.323 

0,253 

0.467 

0.392 

3.164 

83.2 

0,0189 

0.0138 

0.0193 

0.0226 

1 .40 

0.0163 

0 ,0067 

3 37 

0.259 

27.80 

0.329 

0.222 

0.614 

0.551 

3.441 

82.6 

0.0178 

0.0147 

0.0213 

0.0250 

1 .40 

0.0190 

0,0076 

40 

0,246 

27.60 

0.169 

0,199 

0.411 

0.461 

3.677 

133.5 

0.0148 

0.0154 

0.0228 

0.0254 

1 .41 

0.0143 

0,0059 

41 

0.427 

26,80 

0.325 

0.184 

0.684 

0.491 

3,848 

44.4 

0.0212 

0.0229 

0.0314 

0.0328 

1 .40 

0.0220 

0.0088 

42 

0.327 

26.80 

0.313 

0.162 

0.594 

0,605 

4.130 

28,7 

0.0301 

0,0260 

0.0368 

0.0378 

1 .40 

0.0275 

0.0110 

43 

0.415 

26.20 

0.103 

0.267 

0.271 

0.338 

3.053 

527.1 

0,0089 

0,0105 

0.0177 

0.0204 

1 .55 

0.0052 

0.0051 

44 

0,346 

26.50 

0,113 

0,270 

0.243 

0.42S 

3.030 

518.6 

0,0101 

0.0110 

0.01S8 

0.0205 

1 .55 

0.0093 

0,0051 

45 

0,450 

25.80 

0.448 

0,220 

0,661 

0.452 

3.460 

59,2 

0,0238 

0.0185 

0.0238 

0,0261 

1.40 

0.0230 

0,0092 

46 

0,358 

26.30 

0,468 

0,181 

0.651 

0.652 

3.885 

37.4 

0.0265 

0.0226 

0.C285 

0.0317 

1 .40 

0.0285 

0.01 14 

47 

0.380 

26.90 

0.292 

0.195 

0,703 

0.551 

3.721 

70.0 

0.0209 

0.0195 

0.0264 

0.0276 

1.40 

0,0218 

0.0087 

46 

0.380 

26.80 

0.146 

0.202 

0.373 

0.304 

3.644 

247.5 

0.0131' 

0.0165 

0,0247 

0.0295 

1.45 

0.0130 

0.0059 

49 

0.355 

27,35 

0.298 

0.163 

0.676 

0.425 

4.116 

71 .4 

0,0228 

0,0203 

0.0285 

C.0337_ 

1 .40 

0,0212 

0.0083 

50 

0,396 

27.45 

0,301 

0,145 

0,684 

0.390 

4,376 

46 . 6 

0.0241 

0.0258 

0.0365 

0.0383 

1.40 

0,0268 

0,0107 

51 

0.378 

27.25 

0.307 

0.134 

0.694 

0.446 

4.552 

27.8 

0,0315 

0.0307 

0.0430 

0.0410 

1 .40 

0.0330 

0,0132 

52 

0.246 

26.30 

0.158 

0.249 

0.514 

0,465 

3.198 

382.8 

0,0108 

0.0105 

0,0167 

0,0192 

1.51 

0.0113 

0,0057 

53 

0.480 

26,30 

0.153 

0.276 

0.399 

0,337 

2.985 

552.3 

0.0100 

0.0128 

0.0193 

0,0195 

1.55 

0.0102 

0,0055 

54 

0.335 

26.30 

0.105 

0,197 

0,303 

0.372 

3.699 

470,6 

0.0094 

0,0133 

0.0219 

0.0232 

1.53 

0.0108 

C.0057 

55 

0.432 

26,50 

0.106 

0.185 

0.291 

0,305 

3.837 

331.8 

0.0108 

0.0124 

0.0209 

0.0240 

1 .49 

0.0117 

0,0057 

56 

0,368 

27.10 

0.463 

0.270 

0.641 

0.414 

3.030 

93.1 

0,0208 

0.0159 

0,0216 

0.0230 

1.40 

0,0190 

0.0076 

57 

0,481 

27,50 

0.467 

0.294 

0.607 

0.279 

2.858 

116.1 

0.0197 

0.0151 

0.0201 

0.0226 

1,40 

0,0175 

0.0071 

58 

0.492 

27.30 

0.336 

0.216 

0,723 

0.242 

3.499 

107.1 

0,0168 

0,0172 

0.0241 

0,0281 

1.40 

0.0181 

0.0072 

59 

0.335 

26.50 

0.243 

0.146 

0,472 

3.248 

775,1 

0.0102 

0.0105 

0.0182 

0.0187 

1.62 

0.0094 

0.0058 

6C 

0.603 

25.60 

0,500 

0,230 

0.639 

0,175 

3.365 

111.7 

0.0205 

0.0208 

0.0271 

0.0293 

1 .40 

0,0205 

0,0082 

61 

0.535 

26.30 

0,492 

0,202 

0.599 

0.281 

3.644 

80.9 

0,0207 

0.0212 

0.0286 

0.0312 

1.40 

0,0232 

0,0093 

62 

0.397 

26.50 

0,501 

0.177 

0,644 

0.412 

3.934 

55,1 

0.0261 

0.0254 

0,0342 

0.0394 

1.40 

0.0277 

0,0111 

63 

0.281 

27,40 

0.429 

0.179 

0.642 

0,542 

3.909 

181.1 

0.018S 

0.0201 

0.0285 

0,0279 

1.43 

0,0195 

0,0034 

64 

0.401 

28,15 

0.455 

0.186 

0.596 

0,286 

3.825 

198.7 

0.0182 

0,0203 

0,0280 

0.0317 

1 .44 

0.01S7 

C. 0082 

65 

0.421 

28,15 

0,422 

0.192 

0,456 

0,281 

3,755 

208.9 

0.0202 

0.0190 

0.0273 

0.0265 

1.44 

0.0182 

O.OOSO 

66 

0.421 

28.15 

0.364 

0,174 

0.652 

0.359 

3.972 

181.3 

0.0183 

0,0220 

0.0312 

0.0261 

1,43 

0.0195 

0,0084 

67 

0.387 

28,00 

0.371 

0,178 

0,686 

0,409 

3,921 

212.1 

0.0157 

0,0199 

0.0286 

0.0279 

1.44 

0,0180 

0.0030 

68 

0.385 

27.50 

0.163 

0.240 

C.  292 

0,284 

3,275 

1452.0 

0,0092 

0.0106 

0,0184 

0.0215 

1.76 

0.0093 

0,0071 

69 

0.398 

27.20 

0.163 

0.214 

0,644 

0,393 

3.519 

1033.0 

0,0094 

0,0109 

0,0173 

0.0207 

1 . 63 

0.0100 

0,0068 

70 

0.290 

27.00 

0,116 

0.196 

0.274 

0,309 

3.710 

1796.0 

0.0090 

0.0108 

0,0185 

0.0204 

1.81 

0.0038 

0.0071 

71 

0.323 

27.20 

0,117 

0.148 

0,274 

0.331 

4.331 

770.0 

0,0116 

0,0119 

0.0206 

0,0208 

1.62 

c.oiio 

0,0067 

72 

0,315 

27.20 

0.164 

0.198 

0,593 

0.408 

3.688 

762.0 

0,0102 

0,0118 

0,0187 

0,021 1 

1.61 

0.0110 

C.0067 

73 

0,373 

26.10 

0.163 

0.190 

0.262 

0.435 

3.778 

733,0 

0.0110 

0.0150 

0 .0241 

0.0236 

1.61 

0,0112 

0,0058 

74 
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Sc honf e Id's  equation  [37]  (Table  2,  column  29)  and  also  on  the  basis  of 
Equation  (3.4.17)  which  is  derived  from  the  aforementioned  Equation 
(3.4.16)  for  Tg  = 0,  that  is,  for  conditions  similar  to  those  assumed 
by  Schijf  and  Schonfeld  [37]  (Table  2,  column  28).  It  is  observed  that 
there  is  close  agreement  between  the  values  of  these  last  columns. 

The  values  of  f.  in  columns  28  and  29  exceed  the  corresponding 
values  in  columns  26  and  27,  based  on  direct  measurements  and  on  the 
rigorous  integration  of  the  equation  of  motion,  respectively , by  a 
margin  of  the  order  of  50  percent.  This  gives  an  idea  of  the  order  of 
magnitude  of  error  introduced  by  the  Schi jf-Schonfeld ' s [37]  simplified 
model . 

The  correlation  of  f. , as  obtained  by  each  one  of  the  above 

methods,  to  the  pertinent  flow  and  salinity  parameters  followed  the 

general  line  of  the  author's  early  work  [11].  After  several  trials  it 

2 

was  found  that  f..  can  be  best  correlated  to  the  number  ReFr  with  the 
relative  density  difference  Ap/(p  + Ap)  taken  as  a parameter.  Here  Fr 
is  the  nondensimetric  Froude  number. 


2 

The  dimensionless  variable  ReFr  is  given  by  the  equation 


(6.4.9) 


where 


I 


10 


(6.4.10) 


[6(Frg)4/3  + 3(Fr^)2/3  + 1] 


Ip  is  graphically  given  in  Figure  3.3. 
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The  correlations  for  the  f.  computed  from  Equations  (6.4.8) 
and  (3.2.17)  are  presented  in  Figures  6.23  to  6.30.  In  the  same 
figures  values  of  f.  computed  on  the  basis  of  Equations  (3.4.16)  and 
(3.4.17)  have  been  plotted.  Figures  6.31  and  6.32  summarize  the  curves 
of  Figures  6.23  to  6.30  corresponding  to  the  cases  TQ  f 0 and  TQ  = 0, 
respectively. 

In  Figures  6.33  (a),  (b),  and  (c),  data  from  other  investiga- 
tions are  compared  with  the  corresponding  data  of  the  present  study. 

In  Figures  6.33  (a)  and  6.33  (b),  Powell's  [34]  experimental  results 
for  Ap/(p  + Ap)  = 0.022  have  been  replotted  and  compared  with  the 
results  of  this  study.  When  f.  is  plotted  against  the  Reynolds  number 
only,  as  done  by  Powell  [34],  the  two  sets  of  data  deviate  strongly 

[Figure  6.33  (a)];  however,  a close  agreement  is  reached  when  both  sets 

2 

are  plotted  in  Figure  6.33  (b)  against  ReFr  . In  particular,  the  drop 

2 

of  the  curve  for  low  values  of  ReFr  agrees  with  the  trend  observed  in 
the  author's  experiments.  Powell  [34]  gave  as  a law  of  variation  of  f. 
the  straight  line  given  in  Figure  6.33  (a)  and  he  claimed  that  the  point 
for  Re  = 22,000  seemed  to  be  in  error.  This  point,  however,  as  we  see 
in  Figure  6.33  (b),  is  in  agreement  with  the  general  law  derived  from 
this  study  when  plotted  against  the  appropriate  parameter  (Figure  6.31). 

The  curve  corresponding  to  Ap/(p  + Ap)  = 0.020  has  been 
replotted  in  Figure  6.33  (c)  together  with  the  reanalyzed  laboratory 
data  by  Abraham  and  Eysink  [1]  and  one  field  data  point  from  the  Mis- 
sissippi estuary  [31].  The  agreement  is  quite  good.  This  suggests  that 
until  the  curves  of  Figure  6.31  are  implemented  with  field  data,  they  can 
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Figure  6.23  Correlation  of  f.  with  ReFr  and  [(Ap)/(p  + Ap)]  (group  1) 
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Figure  6.24  Correlation  of  f.  with  ReFr  and  [ ( Ap ) / ( p + Ap)]  (group  2) 
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Figure  6.25  Correlation  of  f.  with  ReFr  and  [(Ap)/(p  + Ap)]  (group  3) 
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Figure  6.26  Correlation  of  f.  with  ReFr  and  [(Ap)/(p  + Ap)]  (group  4) 
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Figure  6.27  Correlation  of  f.  with  ReFr  and  [(Ap)/(p  + Ap)]  (group  5) 
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igure  6.28  Correlation  of  f.  with  ReFr  and  [(Ap)/(p  + Ap)J  (group  6 
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Figure  6.29  Correlation  of  f.  with  ReFr  and  [(Ap)/(p  + Ap)]  (group  7) 
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Figure  6.30  Correlation  of  f.  with  ReFr  and  [ (Ap)/(p  + Ap)]  (group  8) 
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Figure  6.31.  Variation  of  f.  with  ReFr  and  [ (Ap)/(p  + Ap)]  for 
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Figure  6.32  Variation  of  f.  with  ReFr  and  [(Ap)/(p  + Ap)]  for 
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(b) 

Figure  6.33  Comparison  of  data  from  other  investigations  with  the 
corresponding  data  of  the  present  study 


0.08h  • ABRAHAM/EYSINK  (1970) 
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Figure  6.33— Continued 
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be  tentatively  extrapolated  following  the  extrapolation  of  the 
Ap/(p  + Ap)  = 0.020  curve  of  Figure  6.33  (c). 

The  bed  friction  coefficient  fg  [Equation  (3.4.15)]  was 
evaluated  in  the  following  way.  For  each  experiment  the  ratio 
X = hs/hs]  [Figure  3.4  and  Equation  (3.4.14)]  has  been  computed.  In 
Figure  6.34,  the  quantity  X - 1 = Tg/T.  is  plotted  against  ReFr2.  It 
is  observed  that  this  ratio  is  independent  of  the  salinity  parameter 

Ap/ (p  + Ap)  and  that  it  reaches  a constant  value  equal  to  0.4  for  ReFr2 

2 ? 

= 100.  For  ReFr  > 100  the  ratio  Tg/T^  increases  with  ReFr  . In  particu- 

2 

lar,  for  ReFr  = 2,000  the  bottom  shear  stresses,  Tg,  are  equal  to 

about  85  percent  the  value  of  the  corresponding  interfacial  shear 

stresses,  . Knowing  X and  f . , the  bed  friction  factor,  fg,  can 

readily  be  computed.  In  the  computation  of  fQ,  the  values  of  f.  have 

been  taken  from  the  graph  of  Figure  6.31  (f..  = f..  ) and  are  tabulated 

in  column  31  of  Table  2.  Likewise,  the  values  of  X have  been  taken 

equal  to  X , that  is  equal  to  the  value  correspond!' ng  to  the  curve  of 

Figure  6.31.  These  values  are  tabulated  in  column  30  of  Table  2. 

2 

The  correlation  of  fg  with  ReFr  for  different  values  of 
Ap/ (p  + Ap)  is  given  in  Figure  6.35. 

The  graphs  of  Figures  6.31  and  6.34  or  6.35  can  be  applied  to 
practical  problems,  for  the  calculation  of  f.  or  fg,  because  the  intro- 
duced parameters  can  be  readily  defined.  A comparison  of  Figure  6.35 
with  Figure  2.5  leads  to  the  justification  of  the  earlier  conclusion, 
based  on  the  review  of  earlier  work,  that  both  the  interfacial  and  bed 
friction  factors  cannot  be  universally  expressed  as  a function  of  one 
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Figure  6.34  Variation  of  Tn/T.  with  ReFr 
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of  the  commonly  used  dimensionless  parameters  only  (Re,  Fr^,  or 
Keulegan's  number). 


According  to  section  6.3,  the  experiments  of  this  study  indi- 
cated that  the  velocity  distribution  follows  a linear  law  within  a 
relatively  thin  zone  around  the  velocity  interface,  while  above  and 
below  that  zone  the  velocity  distribution  becomes  logarithmic. 

The  objective  of  this  section  is  the  detailed  analytical 
formulation  of  these  laws.  This  was  accomplished  in  the  following  way. 
First,  from  the  values  of  f..  (Table  2,  column  31)  and  from  the  mean 
velocity  V-j,  of  the  upper  layer  at  the  cross  section  of  measurements, 
the  friction  velocity,  u*,  was  computed  for  each  experiment  on  the 
basis  of  the  equation 


The  velocity  U/2,  corresponding  approximately  to  the  velocity  at 
the  interface,  has  been  taken  as  a reference  velocity.  The  velocity 
difference  u - U/2  was  next  nondimensional ized  by  u*;  the  corresponding 


The  dimensionless  velocities  have  been  plotted  in  Figures  6.36  to  6.43 
against  the  dimensionless  vertical  distance  with  Ap/(p  + Ap)  as  a 
parameter.  Here  bars  indicate  average  values.  These  plots  lead  to  the 
following  law  for  the  interface: 


6.5  Velocity  Distribution  Laws 


(6.5.1) 
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Figure  6.36  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (group  1) 
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Figure  6.37  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (group  2) 
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Figure  6.38  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (group  3) 
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Figure  6.39  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (group  4) 
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Figure  6.40  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (group  5) 
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Figure  6.41  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (group  6) 
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Figure  6.42  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (group  7) 
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Figure  6.43  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (group  8) 
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-u  - u/2  = X £n 
u*  k3 


u*(y  - ynT 


v 


+ C- 


for  the  upper  zone,  Su,  and 


u - U/2  _ 1 


£n 


u*(yQ  - yT 


v 


+ c. 


(6.5.2) 


(6.5.3) 


for  the  lower  zone  6u-| , where  k3,  k^,  c^,  and  c3  are  coefficients  to  be 
defined  later.  As  a first  approximation,  the  mean  values  of  these 
coefficients  from  the  graphs  of  Figures  6.36  to  6.43  are  taken  and 
introduced  into  Equations  (6.5.2)  and  (6.5.3)  which  then  take  the 
form 


u - U/2 


0.37 


In 


My 


yor 


- 1.61 


(6.5.4) 


for  the  upper  zone  and 


U/2  _ 


0.24 


£n 


u*(yn  - yT 


v 


+ 3.70 


(6.5.5) 


for  the  lower  zone. 

A more  careful  study  of  the  graphs,  however,  leads  to  the  obser- 
vation that  the  dimensionless  velocity  distribution  curve  differs  for 
each  experiment  (examples  in  Figures  A. 13  to  A. 19  of  the  Appendix). 

More  specifically,  the  slope  of  the  straight  line  portion  of  each 
experiment,  sometimes,  deviates  from  the  slope  of  the  mean  straight  line 
of  the  correspond!- ng  group  of  experiments.  The  analysis  of  the 
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experimental  data,  on  the  basis  of  straight  lines  of  each  experiment 
instead  of  the  mean  line  of  the  group  is  presented  in  Table  3,  where 
a2,  a3,  b2,  and  b3  are  coefficients  of  Equations  (6.4.3)  and  (6.4.4). 
Then  k2,  k^,  c2,  and  c3  of  Equations  (6.5.2)  and!  (6.5.3)  are  correlated 
with  the  flow  and  salinity  parameters.  The  correlations  for  k2  and  k^ 
are  given  in  Figures  6.44  and  6.45,  respectively,  and  the  correlations 
for  c2  and  c3  are  presented  in  Figure  6.46.  The  straight  lines  of 
Figure  6.44  can  be  expressed  by  the  relations 


and 


The  equations  for  the  curves  of  Figure  6.46  are 

c3  = 9.50  in  |k3|  + 7.97  (6.5.8) 

and 

c2  = - 6.37  in  |k2|  - 4.71  (6.5.9) 

The  correlation  coefficient  of  Equations  (6.5.8)  and  (6.5.9) 
appears  to  be  satisfactory.  From  the  graphs  of  Figures  6.44  and  6.45 
and  from  Equations  (6.5.6)  and  (6.5.7)  it  can  be  concluded  that  k2 


0.77 


v/6u1/hs  V/ReFr‘' 


1/2 


(6.5.7) 


0.77 


(6u/hQ)(ReFr2)0-4 


1/3 


(6.5.6) 


Table  3.  Correlation  of  k and  c coefficients  with  the  flow  and  salinity 
parameters 
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Figure  6.44  Correlation  of  k?  with  the  flow  and  salinity  parameters 


GR.  SYM.  GR.  SYM. 


178 


o 


Figure  6.45  Correlation  of  k,.  with  the  flow  and  salinity  parameters 
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Figure  6.46  Correlation  of  k9  and  k~  with  c9  and 
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and  k^  are  not  constants  but  they  depend  on  the  flow  and  salinity 
parameters.  Therefore,  Csanady's  conclusion  [9]  that  k3  is  constant 
and  equal  to  von  Karman's  universal  constant  is  not  correct.  In  the 
same  way  c2  and  c^  depend  on  the  same  parameters,  since  they  are 
expressed  in  terms  of  k2  and  k3-  The  graphs  of  Figures  6.36  to  6.43 
or  the  graphs  of  Figures  6.44  to  6.46  and  Equations  (6.5.2)  and  (6.5.3) 
clearly  define  the  velocity  distribution  in  the  zones  Su  and  6u-j . A 
first  approximation  of  the  velocity  distribution  in  the  zones  6u  and  Su-j 
is  given  by  Equations  (6.5.4)  and  (6.5.5),  where  k and  c are  constants. 

A better  approach  to  the  problem  is  obtained  by  the  equations  of  Figures 
6.36  to  6.43  where  k and  c are  functions  of  Ap/(p  + Ap).  These 
procedures  can  be  readily  applied  to  practical  problems.  Another 
approach  is  obtained  by  means  of  the  Equations  (6.5.2)  and  (6.5.3) 
through  Equations  (6.5.6)  to  (6.5.9). 
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CHAPTER  VII 

SUMMARY  AND  CONCLUSIONS 


The  dynamics  of  arrested  saline  wedges  and,  specifically,  the 
functional  dependence  of  the  shear  stresses  at  the  interface  and  at  the 
bed  and  the  velocity  distribution  in  the  neighborhood  of  the  interface 
have  been  experimentally  investigated.  The  experiments  were  conducted 
in  a special  open  flume  equipped  with  all  the  accessories  for  the 
necessary  measurements  and  control.  A total  of  74  experiments  have 
been  conducted  for  a wide  range  of  values  of  the  primary  independent 
parameters.  The  relative  density  difference  Ap/(p  + Ap)  ranged  between 
8 °/oo  and  93  °/oo.  The  analysis  of  the  experimental  data  of  the 
present  work  confirmed  an  earlier  derived  law  on  the  density  distribu- 
tion [34]  of  the  form 


p£  ~ p = 1 
p + Ap  2 


tanh 


(7.1) 


The  interfacial  friction  coefficients,  f. , have  been  computed  by 
four  different  approaches: 

1.  From  direct  measurements  of  the  velocity  profiles  and 
Reynolds  stresses  near  the  interface  on  the  basis  of 
the  equation  [Equation  (6.4.8)] 


T 


imax 


f du 

»y 


- pu  v 


imax 


fipVl 

8 


(7.2) 
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where  x.  is  the  maximum  value  of  the  shear  stresses 
imax 

in  the  neighborhood  of  the  interface. 

2.  Through  the  integration  of  the  equations  of  motion  in 
the  vertical  direction  and  subsequent  averaging  through 
a second  integration  in  the  horizontal  direction  over 
the  actual  intrusion  length,  Lg,  making  use  of  the 
measured  velocity  profiles,  the  shape  of  the  interface 
and  other  pertinent  quantities.  This  equation  finally 
takes  the  form  [Equation  (3.4.16)] 


Ti . 

ffApl 
.1  p 

gh$g  + 26(1  + Ap/p)e2v20 

h l oh 1 " 2aV0h0hs 

1 

C\J  o 
> 

2V0h10(Xhl  + hs> 

0.0559  u2 
b1.25  h0 


v 


V0h0 


0.25 


(b  + 2h1) 


(Xh1  + h$) 


V,)  2 


(7.3) 


where  bars  indicate  mean  values  of  the  correspondi ng 
parameters . 

3.  From  Schijf  and  Schonfeld's  equation  [37]  [Equation 


(3.2.17)] 


2h, 


f.  = 

l 


5(FrJ) 


2 - 2 + 3(Fr^)2/3-|(Fr’)4/3 


(7.4) 
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where  Fr^  is  the  densimetric  Froude  number  upstream  of 
the  wedge. 

4.  From  a simplification  of  Equation  (7.3)  for  tq  = 0, 

an  assumption  used  by  Schijf  and  Schonfeld  [37]  in  their 
mathematical  model.  Equation  (7.3)  reduces  then  to  the 
form  (3.4.17) 


(Ap/p)gh1hs0 


2V0h0 


(7.5) 


The  bed  friction  coefficient,  fQ,  has  been  computed  on 
the  basis  of  a linear  relation  between  interfacial  and 
bed  shear  stresses,  that  is, 

■f0pv? 

T0  = U - 1)T,  - -V-  (7.6) 

This  relationship  was  utilized  in  the  derivation  of  Equation 
(7,3),  The  analysis  of  the  experimental  data  on  the  coefficients  f.  and 
fg  led  to  the  following  conclusions: 

1.  The  best  correlation  of  f..  and  fg  with  the  flow  and 
salinity  parameters  has  been  obtained  when  these  coef- 
ficients were  plotted  against  the  dimensionless 
2 

parameter  ReFr  = 4V-j/vg  with  Ap/(p  + Ap)  as  an  inde- 
pendent parameter  (Figures  6.23  to  6.32  and  Figure  6.35). 
The  values  of  fg  are  substantially  lower  than  those  of 
f-.  This  conclusion  and  the  derived  functional 
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dependence  of  both  coefficients  on  these  two  distinct 
parameters  are  drastically  different  from  the  results 
of  all  known  previous  studies,  but  in  agreement  with 
the  author's  earlier  work  [11]  regarding  f.  alone. 

2.  The  graphs  of  Figures  6.23  to  6.30  give  the  variation 

2 

of  f.  with  respect  to  ReFr  for  different  values  of 
Ap/(p  + Ap).  The  lower  curve  of  each  graph  is  based  on 
the  approach  gi ven  by  Equation  (6.4.8)  or  Equation  (7.2) 

The  agreement  with  the  results  based  on  Equation  (3.4.16) 
or  Equation  (7.3)  is  very  good.  The  upper  curve  is  based 
on  the  approach  given  by  Equation  (3.4. 17)  or  Equation  (7.5) 
where  tg  is  considered  equal  to  zero.  The  agreement  with 
the  results  based  on  Schijf  and  Schonfeld's  mathematical 
model  [37]  [Equation  (3.2.17)  or  Equation  (7.4)]  is  also 
very  good.  Figures  6.31  and  6.32  contain  the  summary  of  the 
lower  curves,  corresponding  to  tq  f 0,  and  the  upper 
curves  corresponding  to  xQ  = 0 of  Figures  6.23  to  6.30, 
respectively.  From  these  graphs  the  error  in  the  values  of 
f - can  be  readily  estimated  when  the  Schijf  and  Schdnfeld 

O 

mathematical  model  [37]  is  used.  For  example,  for  ReFr 
= 200  and  Ap/(p  + Ap)  = 22.5  °/oo  the  relative  error  is 
54  percent  of  the  real  value. 

3.  The  graph  of  Figure  6.34  gives  the  law  of  variation  of 
the  ratio  of  bed  shear  stresses  to  the  interfacial  shear 
stresses,  Tg/T...  It  is  observed  that  (a)  the  ratio  Tg/T. 
is  independent  of  the  density  variation  and  (b)  the  above 
ratio  is  constant  and  equal  to  0.4  for  low  values  of 
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ReFr2,  especially  for  ReFr2  ^ 100,  while  for  ReFr2 
> 100,  Tq/T.  is  approximately  given  by  the  equation 

Jin  (ReFr2)  - 0.26  (7.7) 


with  a correlation  coefficient  of  95.2  percent. 

4.  The  law  of  variation  of  fQ  with  the  flow  and  salinity 
parameters  is  given  by  the  graph  of  Figure  6.35.  The 
values  of  fg  were  calculated  from  the  following  equation 

fQ  = f;U  - 1)  (7.8) 


which  was  derived  from  Equation  (7.6)  with  T^  substituted 
~2 

by  f -} pV-j /8 . The  val-ues  of  f.  and  A - 1 were  taken  from  the 
graphs  of  Figures  6.31  and  6.34,  respectively. 

The  laws  given  by  the  graphs  of  Figures  6.31,  6.34,  and  6.35  can 
be  easily  applied  to  practical  problems. 

Considering  the  flow  structure  within  the  interfacial  zone,  the 
following  results  and  conclusions  have  been  obtained: 

1.  There  exists  a relatively  thin  core  near  the  density 
interface  within  which  the  dimensionless  velocity  dis- 
tribution is  linear  (Figures  6.10  and  6.11).  This 
distribution  is  given  by  the  equation 


u - U/2 
U 


= 5 


6U2 


(7.9) 


The  symbols  have  been  defined  in  Figure  6.4.  The  linear 
law  appears  to  be  independent  of  salinity. 
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2.  Outside  of  the  linear  zone  the  dimensionless  velocity 
distribution  follows  a logarithmic  law  (Figures  6.12 
to  6.19  and  Figures  6.36  to  6.43)  which  can  be  expressed 
for  the  zone  6u  (Figure  6.4)  by  the  equations 


u 

U 


a^  £n 


(7.10) 


or 


u - U/2 
u* 


u*(y 


v 


(7.11) 


and  for  the  zone  6u^  (Figure  6.4)  by  the  equations 


u^ 

U 


&n 


y0  - y 


56u, 


+ b. 


(7.12) 


or 


fu*(y0  - y)'1 


V 


+ Cr 


(7.13) 


where  the  values  of  the  parameters  a^,  &2>  b^,  b^,  and 
^3’  ^2’  c3’  c2  are  91ven  figures  6.12  to  6.19  and 
Figures  6.36  to  6.43,  respectively,  with  each  figure 
corresponding  to  a particular  value  of  the  relative 
density  difference.  In  each  of  these  graphs,  and  par- 
ticularly in  Figures  A. 6 to  A. 19  of  the  Appendix,  two 
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intersected  straight  lines  of  different  slopes  can  be 
distinguished  in  each  of  the  two  zones,  <$u  and  Su-p 
3.  In  the  zone  6u  the  intersection  of  the  straight  lines 
occurs  at  a relative  distance  (y  - yg)/Su  of  approxi- 
mately 0.15  (Figures  6.12  to  6.19  and  Figures  A. 6 to 
A. 12  of  the  Appendix).  The  first  straight  line  cor- 
responding to  the  logarithmic  distribution  immediately 
next  to  the  linear  zone  and  extending  until  the  point 
where  (y  - yg)/Su  = 0.15,  constitutes  the  "law  of  the 
interface"  of  the  upper  layer.  This  law  is  analogous 
to  the  wall  law  for  homogeneous  fluids  and  solid 
boundaries.  In  the  zone  6u^  the  intersection  of  the 
two  straight  line  segments  occurs  at  (yQ  - y)/56u-j 
=0.1.  The  straight  line  immediately  next  to  the 
linear  zone  and  up  to  the  previous  point  of  inter- 
section constitutes  the  law  of  the  interface  of  the 
lower  layer.  The  laws  of  the  interface  can  be 
expressed  by  the  following  approximations: 

a.  A first  approximation  is  given  by  the  equa- 
tion 


u - U/2  = 1 

u+  "0.37 


£n 


My  - y0) 


V 


1.61 


(7.14) 


for  the  zone  6u,  and  by  the  equation 
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u - U/2  = 
u* 


(7.15) 


for  the  zone  6u-| . Here  the  parameters  k^, 
c^j  and  c2  of  Equations  (7.11)  and  (7.13)  are 
taken  constant  and  equal  to  the  mean  of  their 
corresponding  values  from  Figures  6.36  to  6.43. 

b.  A second  approximation  is  given  by  the  graphs  of 
Figures  6.36  to  6.43,  where  the  previous 
parameters  appear  to  be  functions  of  the  rela- 
tive density  difference  Ap/(p  + Ap). 

c.  A third  approximation  can  be  obtained  from  the 
graphs  of  Figures  6.44  to  6.46  where  the 
parameters  k^,  k^,  c^,  and  are  given  as  follows: 


k 


(Su/h0)(ReFr2)0'4 


(7.16) 


k 


2 


and 


9.50  £n  |k3|  + 7.97 


(7.17) 


-6.37  £n  |k2|  - 4.71 
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It  follows  that  there  are  two  laws  for  the 
interface  for  each  layer,  respectively.  The 
parameters  and  are  not  constant,  as 
Csanady  [9]  claims,  but  functions  of  the  flow 
and  salinity  parameters. 

4.  The  straight  line  in  the  semi  log  plot,  which  is  defined 
between  the  points  of  (y  - yg)/6u  = 0.15  and 
(y  - yg)/5u  = 1 for  the  zone.du  and  between  the  points 
(yQ  - y)/56u-|  = 0.1  and  (yQ  - y)/56u-j  = 0.22  for  the 
zone  6u-j  defines  the  "outer  region  law"  for  each  layer 
(Figures  6.12  to  6.19,  6.36  to  6.43,  and  A. 6 to  A. 19). 
These  laws  are  given  as  functions  of  the  relative 
density  Ap/(p  + Ap)  by  the  corresponding  graphs.  A 
first  approximation  of  the  law  is  given  by  the  equation 


u - U/2 
u* 


0.45 


Jin 


u*(y  -y0F 


v 


+ 0.08 


(7.18) 


for  the  zone  6u  and  by  the  equation 

+ 4.04  (7.19) 

for  the  zone  6u-| . The  described  laws  define  completely 
the  dimensionless  velocity  distribution  in  the  shear 
flow  region  around  the  interface  of  an  arrested  saline 
wedge. 


u - U/2  . 


0.24 


£n 


u*(yn  - yJ" 
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When  the  theoretical  analysis  of  section  3.3  is  applied  to  the 

calculation  of  interfacial  and  bed  shear  stresses,  the  average  velocity 

— ~~2 

of  the  upper  layer,  , and  its  second  power,  V-j , over  the  length  of  the 
arrested  saline  wedge  must  be  taken  from  the  following  equations: 


V*  ■ Vo 


(A.l) 


and 


J-L 


V^dx  = l/Q 


(A. 2) 


where 


5[12(FrJ)2  £n  (Fr^)2/3  + 1 +8(Fr^)2  - 9(Fr' )8/3] 
3[1  - (Frg)2/3]3  • [6(Fr£)4/3  + 3(Fr£)2/3  + 1] 


and 


10[1  -9(Fr’)4/3  + 8(Fr^)2-4(Fr^)2  £n  (Fr^)] 

15  3[1  - (Frg)2/3]3  • [6(Fr^)4/3  + 3(Fr^)2/3  + 1] 

Figures  A.l  to  A. 5 give  the  dimensional  velocity  and  density  distribu- 
tion for  some  experiments  of  group  5 (Table  1),  where  characteristic 
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quantities  of  the  pertinent  parameters  are  defined.  The  data  of  these 
figures  are  summarized  in  Table  4.  Each  page  of  the  table  corresponds 
to  one  experiment  which  is  numbered  in  column  1.  In  this  table  the 
calculated  maximum  value  of  the  sum  + x^,  where  x^  are  the  viscous 
stresses  and  x-t  the  Reynolds  stresses,  is  tabulated.  The  interfacial 
coefficient,  f . , is  calculated  from  (x-  , + x..)  v [Equation  (6.4.8)] 
and  its  value  is  tabulated  in  column  26  of  Table  2. 

In  the  graphs  of  Figures  A. 6 to  A. 19  the  white  circles 
correspond  to  experimental  points  and  the  black  circles  to  interpolated 
points  obtained  from  Figures  A.l  to  A. 5.  From  Figures  A. 5 to  A. 12, 
the  coefficients  a^,  a^,  bg,  and  bg,  of  Equations  (6.4.3)  and  (6.4.4), 
were  calculated  and  tabulated  in  Table  3.  In  the  same  way,  the  coef- 
ficients kg,  kg,  Cg,  and  Cg  of  Equations  (6.5.2)  and  (6.5.3)  were 
calculated  and  tabulated  in  Table  3. 
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Table  4.  Geometrical  and  flow  parameters  of  each 
experiment 


Hxper im . 

u 

Vy 

u 

(u4j-)/u. 

u.(yQ-y ) 

e 

u 

'T 

cm 

cm/sec 

1 

2 

3 

4 

5 

6 

7 

8 

AS1’ 

a 

1.95 

0.700 

0. 194 

0.121 

-9.13 

19,52 

y0-y 

2.00 

0.850 

0.181 

0.147 

-8.51 

18,22 

-2—  = 0.031 
£ 

h =2.57 

2.10 

1.120 

0.153 

0.193 

-7.39 

15.40 

— — 

s 

n=o.09 

2.15 

1.220 

0.140 

0.210 

-6,97 

14.07 

3y^-=  0.21 
U.o 

2.20 

1.380 

0.126 

0,238 

-6,31 

12.62 

p-  = 5.844 
dy 

|hs+n=2.66 

1 

* 

2.25 

1.683 

0,112 

0.290 

-5,05 

11,28 

5u=4. 59 

2.30 

1.720 

0.099 

0.297 

-4,90 

9.96 

ul  = 0,0646 

£=3.65 

2.35 

1.850 

0.085 

0.319 

-4,36 

8.59 

2.40 

2.000 

0.071 

0.345 

-3.73 

7.14 

yo"hs+n 

2.45 

2,200 

0.058 

0.379 

-2.90 

5.80 

y0-y 

= 0.102 

2.50 

2.330 

0.044 

0.402 

-2.37 

4.45 

U=5.80 

2.55 

2.500 

0.030 

0.431 

-1  .66 

3.00 

= 0.20 

Vj=5.23 

y-y0 

u.(y-y0) 

u»0 

^ = 4.504 
dy 

v=0. 00874 

6u 

V 

2.85 

3.580 

0.041 

0.617 

2.82 

5.13 

ReFr  =83.2 

2.90 

3.700 

0.052 

0.638 

3.32 

6.56 

u;  = 0,0523 

— *-=24. 01. 
P+Ap 

* 

3.00 

3.05 

3.950 

4.081 

0.074 

0.085 

0.681 

0.704 

4.36 

4,90 

9.36 

10.76 

y-yn 

3,10 

4.120 

0,096 

0.710 

5,06 

12.18 

-*-2-  = 0.040 
ou 

3.20 

4.280 

0.118 

0.738 

5.73 

14.94 

-t/v  n 

2 

u*max= 

3,30 

4.400 

0.139 

0,759 

6,22 

17.57 

~2~  = 0.29 
u « o 

* 

3.35 

4.503 

0.150 

0,776 

6.65 

19.00 

^ = 3.840 
dy 

=0.0646 

u*=0,254 

3.40 

4.520 

0.161 

0.779 

6,72 

20.93 

3.50 

4.600 

0.183 

0.793 

7,05 

23.19 

u l = 0.0522 

f -=0.0189 

3,70 

4.710 

0,227 

0.814 

7.51 

28.71 

4.00 

4.900 

■0.292 

0.845 

8.30 

36.96 

u^0  = 0,0646 

ifig=0.0168 

::u,g=0,240 

* 

4.45 

5.080 

0.390 

0,876 

9.05 

49,36 

4.50 

5.120 

0.401 

0.883 

9.21 

50.72 

5,00 

5.220 

0,510 

0.900 

9,63 

64,59 

5.50 

5.320 

0.619 

0.917 

10,04 

78,36 

* 

5,95 

5.295 

0.717 

0,913 

9.94 

90,70 

6.00 

5.400 

0.728 

0.931 

10.37 

92.13 

6,50 

5.500 

0.837 

0,943 

10,79 

105.93 

7,00 

5,580 

0.946 

0.962 

11.12 

119.73 

★ 

7.45 

5,468 

1.044 

0.943 

10.66 

132.14 

7.50 

5,670 

1 .054 

0.978 

1 1 .49 

133.40 

* 

9.68 

5.853 

1.530 

1.009 

12.25 

193.64 

t Symbol  of  the  experiment 
a Measured  quantities 

- - f-o  value  of  f.  from  the  graph  and  ut  value  of  u4  from  f . 


Table  4— Continued 


1 

2 

3 

4 

5 

6 

7 

8 

A 

. * 

1.48 

0.469 

0.206 

0,099 

-9.25 

24.78 

1.60 

0.500 

0,182 

0.105 

-9,10 

21.90 

1.70 

0.620 

0.161 

0.131 

-8,52 

19,37 

y0-y 

-7—  = 0.048 

hs=2,03 

1 .80 

0.750 

0.141 

0.158 

-7,89 

16.96 

n=0,46 

* 

1 .88 

0.695 

0.124 

0.146 

-8.16 

14.92 

^ = 0.23 

1.95 

0.950 

0.110 

0.200 

-6.92 

13.23 

u*o 

hc  + n=2, 49 

2.00 

1 .050 

0,100 

0.221 

-6.43 

12.03 

j—  = 3. 50o 

5u=5. 91 

2.10 

1 .200 

0.080 

0.253 

-5.70 

9.62 

dy 

e=4 , 9 

2.15 

1.350 

0.069 

0.284 

-4,98 

8.30 

u 2 = 0.0386 

Vhs+n 

2.20 

1 .450 

0.059 

0.305 

-4.49 

7.10 

* 

2.28 

1.699 

0,043 

0.358 

-3.28 

5.17 

y-yn 

u„(y-yj 

y-yn 

U=4, 75 

0 - 0 040 

ou 

V 

5u 

Vj=4.25 
v=0. 00839 

2.70 

3.250 

0,035 

0.684 

4.25 

5.08 

— 77 

* 

2.78 

3.476 

0.049 

0,731 

5.34 

7.11 

.rUV  . 0 29 
2 

ReFr2=62.6 

2.80 

3.520 

0.052 

0.741 

5.56 

7.55 

u*o 

-on  iy 

2.90 

3.730 

0,069 

0,785 

6.58 

10.01 

h = 2.73C 
dy 

p-t/lp 

3.00 

3.880 

0,086 

0.817 

7.31 

12.48 

3,10 

3.970 

0.102 

0.836 

7.74 

14.80 

u.  = 0.0347 

3.20 

4.040 

0,120 

0.851 

8.08 

17,41 

* 

3,28 

4.  145 

0,132 

0.872 

8.59 

19,15 

3,30 

4.090 

0.135 

0,861 

8.33 

19.59 

u20  = 0.0403 

L 

u*max= 

3.40 

4.130 

0.152 

0.869 

8.52 

22,06 

=0.0403 

3.60 

4.180 

0.186 

0.880 

8.76 

26.99 

u„=0,201 

3.80 

4.200 

0.220 

0.885 

8.86 

31.92 

f .=0.0178 

4.00 

4.48 

4.220 

0.253 

0.887 

8.96 

36.71 

* 

4.298 

0.333 

0.905 

9.33 

48.32 

4.60 

4.300 

0.354 

0.905 

9.34 

51  .37 

5.00 

4.340 

0,421 

0.913 

9.54 

61 .09 

f =0.0190 

5.40 

lg  * 

4.390 

0,488 

0.925 

9.78 

70.81 

u „ = 0. 207 

* 

5.98 

4,414 

0.585 

0,930 

9.90 

84,89 

6.40 

4.470 

0.656 

0.941 

10.17 

95,19 

* 

6.98 

4.51 1 

0,754 

0.949 

10.37 

109,41 

7.40 

4.560 

0.823 

0.960 

10.61 

119.42 

8.00 

4.600 

0,924 

0.968 

10,80 

134.08 

* 

8.48 

4.733 

1,004 

0.997 

11.45 

145.69 

9.00 

4.660 

1,092 

0.981 

11.09 

158,46 

Table  4 — Continued 


1 

2 

3 

4 

5 

6 

7 

8 

A S4 

’1.92 

0.486 

0.234 

0.092 

-10.30 

29.25 

2.00 

0,450 

0.218 

0.085 

-10.48 

27,25 

y0-y 

0 = 0.048 

* 2.12 

0.478 

0.194 

0.090 

-10,34 

24.25 

£ 

O 

VO 

t\j 

ii 

1/ 

XI 

2.20 

0.660 

0.178 

0.125 

- 9,43 

22.25 

= 0.23 

n=0.49 

2,30 

0.820 

0.158 

0.155 

- 8.71 

19.75 

u*o 

h +n-3.09 

2.40 

C.920 

0.138 

0.174 

- 8.24 

17.25 

p-  = 4,179 
dy 

’ 2.52 

0.896 

0.114 

0.169 

- 8.35 

14.25 

6u=7 ,93 

2.60 

1.300 

0.098 

0.245 

- 6,43 

12,25 

u 2 = 0.0440 

£=5.00 

2.70 

1,480 

0,078 

0.279 

- 5.57 

9.75 

yo  = hs1'n 

2.75 

1.600 

0.068 

0.302 

- 5.00 

8.50 

* 2,82 

1.946 

0,054 

0.367 

- 3.35 

6.75 

y-y. 

0 = 0,023 

2.85 

1.950 

0.048 

0.368 

- 3,33 

6.00 

Ou 

U= 5 . 30 

y-y0 

u.(y-y0 ) 

= 0.20 
u.o 

Vj  =4 .88 

5u 

V 

v=0. 00840 

3.30 

3.200 

0.027 

0.604 

2.62 

5.35 

— = 4,315 

ReFr2=133.5 

3.40 

3.460 

0.039 

0.653 

3.86 

7.73 

dy 

p+Ap  ‘-3,8"> 

3.50 

3.650 

0.052 

0,689 

4.76 

10.31 

u,  = 0.0440 

’ 3.62 

3.951 

0.067 

0.745 

6,20 

13.28 

3.70 

3.960 

0.077 

0.744 

6.24 

15.27 

3.80 

4.100 

0.090 

0.774 

6.90 

17.84 

u20  = 0,0388 

4.00 

4.350 

0.115 

0,821 

8.10 

22.80 

2 

u*max= 

4.20 

4.480 

0.140 

0.845 

8.71 

27.76 

=0,0440 

4.30 

4.540 

0.153 

0.857 

9.00 

30.33 

’ 4.42 

4.574 

0.168 

0.863 

9.16 

33.31 

u. =0.210 

4.60 

4.650 

0.190 

0.877 

9.52 

37.67 

f i =0 .01 48 

5.00 

4.730 

0.241 

0.892 

9.90 

47.78 

5.50 

4.820 

0.304 

0.909 

10.33 

60.27 

6.00 

4.890 

0,367 

0.923 

10.67 

72.76 

f -=0.0143 

6.50 

85.25 

4,960 

0.430 

0,936 

1 1 .00 

u,„=  0,206 

’ 6.92 

4.960 

0.483 

0.936 

11,00 

95.75 

7.50 

5.040 

0.556 

0.950 

11.38 

110.23 

8.50 

5.090 

0.682 

0.960 

1 1 .62 

135,21 

9.50 

5.120 

0.808 

0.968 

11.76 

160.19 

10.50 

5.160 

0.934 

0.974 

11.95 

185.17 

11.50 

5.200 

1,061 

0.981 

12.14 

210.34 

’12.42 

5,167 

1.177 

0.975 

11  .99 

233.34 

’19.17 

5.296 

2.028 

0.999 

12.60 
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1 

2 

3 

4 

5 

6 

7 

8 

A U4: 

* 

1.42 

0.577 

0,254 

0.135 

-7.97 

24,37 

y -y 

0 -0.045 

1.50 

0.450 

0.236 

0,105 

-8.62 

22.64 

£ 

1,60 

0.550 

0.212 

0.129 

-8-.  1 1 

20,34 



h =2.35 

1.70 

0.650 

0.188 

0,152 

-7.60 

18,04 

^ = 0.23 

u.o 

0=0.14 

* 

1.82 

0,713 

0.160 

0.167 

-7.28 

15.35 

* 3.510 
dy 

hs+n=2.49 

1.90 

0.840 

0.140 

0,196 

-6.63 

13.43 

2.00 

0.950 

0.117 

0.222 

-6.07 

11.23 

u l = 0.0378 

5u  = 6 . 93 
£=4.20 

2.05 

2.10 

1.030 

0.105 

0.241 

-5.66 

10.07 

1.100 

0.093 

0.257 

-5,31 

8.92 

yo=hs+n 

* 

2.12 

1.087 

0.088 

0.254 

-5.37 

8.44 

2.20 

1.350 

0.069 

0.315 

-4a  03 

6.62 

y-yn 

2.25 

1 .470 

0,057 

0.343 

-3.42 

5.47 

, 0 = 0.023 
Ou 

U=4 ,28 

2.30 

1.620 

0.045 

0,378 

-2.65 

4.32 

”uv  o °n 

V1=3.78 

u.(y-yQ) 

y-y0 

u.o 

v=0. 00858 

6u 

V 

du  2.511 

ReFr^-44 ,4 

2.60 

2.500 

0.016 

0.584 

1 .84 

2.53 

dy 

-OC  1<V 

p+Ap  Z5, 

2.70 

2.720 

0.030 

0,636 

2,96 

4.75 

u l = 0.0282 

* 

2,72 

2.767 

0.033 

0.646 

3,20 

5.22 

2,80 

2. 850 

0.045 

0,666 

3.62 

7.12 

2.90 

2.950 

0.059 

0. 689 

4.13 

9.34 

u20  = 0.0329 

3.00 

3.030 

0,074 

0.708 

4.54 

11,71 

u„max=0,0378 

3.15 

3,170 

0.095 

0.749 

5.26 

15.04 

u* =0,194 

3,30 

3.270 

0.117 

0.764 

5.77 

18.52 

3.50 

3.380 

0.146 

0.790 

6.33 

23.11 

f =0.0212 

3.70 

3,450 

0.174 

0.806 

6.68 

27.54 

3.90 

3.500 

0.203 

0,818 

6.94 

32.14 

4.00 

3.520 

0.220 

0.822 

7.04 

34.83 

f .=0.0220 

4.22 

ig 

* 

3.633 

0.250 

0,849 

7,62 

39.58 

u. =0.198 

4.50 

3,630 

0.290 

0.848 

7.60 

45.91 

g 

5.00 

3.710 

0.362 

0.867 

8.01 

57,31 

5.50 

3,760 

0,433 

0,879 

8.27 

68,54 

6.00 

3.820 

0,506 

0.893 

8.57 

80.10 

7,00 

3,940 

0.651 

0.921 

9.18 

103.53 

* 

7.22 

3.965 

0.682 

0.926 

9.31 

107.97 

8.00 

4,050 

0,795 

0,946 

9.74 

125,85 

9,00 

4.150 

0.940 

0.970 

10.26 

148.81 

* 

9.48 

4.174 

1 .012 

0.975 

10.38 

160.21 

' 
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1 

2 

3 

4 

5 

6 

7 

8 

4 u5: 

* 1.56 

0,569 

0.192 

0.160 

-6.55 

19.04 

yQ-y 

1,60 

0.520 

0,183 

0,146 

-6,81 

18,15 

-y—  = 0.028 

1.70 

0.650 

0.161 

0.183 

-6,11 

15.97 

— 

n s =2 . 02 

1.75 

0.720 

0,150 

0.202 

-5.73 

14,88 

= 0.20 

n=  0.42 

1 .80 

0.750 

0.139 

0.211 

-5.57 

13.79 

u*o 

h -^1  = 2. 44 

’ 1.86 

0.614 

0,126 

0.172 

-6.30 

12.50 

du  . 

1 

-7-  = J.480 

1.90 

0.850 

0.117 

0.239 

-5.03 

11.60 

ay 

6u=5.91 

1.95 

0.920 

0,107 

0.258 

-4.65 

10,61 

u.  = 0.0356 

e=4 .60 

2.00 

1 .020 

0.096 

0.287 

-4,11 

9.52 

1 

Vhs+n 

2.05 

1 .070 

0.085 

0.301 

-3.84 

8.43 

1 

2.10 

1.150 

0.074 

0.323 

-3.41 

7.34 

1 

* 2,16 

1.223 

0,060 

0.344 

-3.01 

5.95 

y-y 

..  0 = 0.023 

u=3. 56 

2,20 

1.350 

0,052 

0.379 

-2.32 

5.16 

2.25 

1.430 

0,041 

0.402 

-1.89 

4,07 

^5  = 0.20 

VjO.os 

2.30 

1.520 

0.030 

0.427 

-1.41 

2.98 

2 

u*o 

v=0, 00858 

y-y0 

u, (y-yQ ) 

^ = 2.620 

ReFr2=28.7 

£ 

V 

dy 

AO  r„ 

2.55 

2.120 

0,019 

0,596 

1.84 

2,42 

u2  = 0.0263 

D+Ac  23,5,<" 

2.60 

2.220 

0,027 

0.624 

2.3S 

3,44 

2.70 

2.400 

0,044 

0,674 

3.35 

5,61 

2,80 

2,520 

0,061 

0.708 

4.00 

7.77 

u20  = 0.0290 

’ 2.86 

2.611 

0,071 

0.733 

4.49 

9,05 

u*max=0,0356 

2.90 

2,620 

0,078 

0.736 

4.54 

9.94 

u, =0.189 

3,00 

2,690 

0.095 

0.756 

4.92 

12.11 

3.20 

2,800 

0.129 

0.786 

5,51 

16.44 

•+> 

II 

O 

o 

OJ 

3.40 

2,860 

0.162 

0.803 

5.84 

20.64 

3,60 

2,920 

0.196 

0.820 

6.16 

24.98 

3.80 

2,950 

0.230 

0.829 

6.32 

29.31 

fi=0.0275 

* 3.86 

2,952 

0.240 

0.829 

6.34 

30,58 

4.00 

2,990 

0.264 

0,840 

6.54 

33,64 

u.|0. 181 

4.50 

3,070 

0.349 

0.862 

6,97 

44,47 

5.00 

3.130 

0.433 

0.879 

7,30 

55,17 

5.50 

3.200 

0.518 

0.899 

7.68 

66.01 

6,00 

3.230 

0.602 

0.907 

7.84 

76,71 

* 6,36 

3,257 

0.663 

0.915 

7.98 

84.49 

7.00 

3,320 

0.772 

0.933 

8.32 

98.38 

8.00 

3,380 

0.941 

0.949 

8.65 

119.91 

9.00 

3.460 

1.110 

0.972 

9.08 

141.45 

* 9.56 

3,534 

1.205 

0.993 

9.48 

53.55 
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1 

2 

3 

4 

5 

6 

7 

8 

A AA2: 

hs=3.00 

T=0.24 

hs+n=3.24 

6u=7,46 

e=6.35 

Vhs+n 

U=8.20 
V1=7.95 
v=0. 00860 
ReFr2=518.6 

^ P _ 0 O AO' 

d+Ap  23,4*‘ 

2 

u*max=0,0799 
u.=  0.283 
f .=0.0101 

f^O.0093 

u.f0.271 

2.10 

2,20 

2.30 

*2.36 

2.40 

2.50 

2.60 

2,70 

‘2,76 

2,80 

2,90 

*2,96 

3,00 

3.05 

0.450 

0,680 

0,900 

1.190 

1.100 

1.400 

1.700 

2.000 

2,127 

2.350 

2.650 

2,848 

3.080 

3.300 

0.180 

0.164 

0.148 

0.139 

0,132 

0,117 

0.101 

0.085 

0.076 

0.069 

0,054 

0.045 

0.038 

0.030 

0.055 

0,083 

0,110 

0,145 

0.134 

0,171 

0.207 

0,214 

0.259 

0.287 

0.323 

0,347 

0.376 

0.402 

-13,27 

-12.44 

-11.64 

-10,58 

-10.91 

- 9.82 

- 8.73 

- 7.64 

- 7.17 

- 6,36 

- 5.27 

- 4.55 

- 3.71 

- 2.91 

36.55 

33.30 

30.05 

28,22 

26.80 

23.76 

20.51 

17.26 

15.43 

14,-91 

10.96 

9.14 

7.72 

6.09 

y -y 

0 = 0,027 

£ 

0.19 

u.o 

^ = 7.700 
dy 

u2  = 0.0799 

6u°  = °*040 

= 0.29 
y.o 

~ = 3.448 
dy 

ul  - 0.0508 
u20  = 0,0730 

y-y0 

u„(y-yQ) 

5u 

V 

3,50 

3.60 
*3.86 

4.00 
4,40 

4.80 

5.20 

5.60 
*5.85 

7.20 

8.00 

8.80 
*9.36 

9.60 
10.40 

*29,86 

4.900 
5.150 
5.660 
5.850 
6.250 
6.650 

6.900 
7.180 
7.326 
7.700 

7.900 
8.000 
8.171 
8.020 
8.030 
8,000 

0.035 

0.048 

0.083 

0,102 

0.155 

0.209 

0.263 

0.316 

0.351 

0.531 

0.638 

0,74.5 

0.820 

0.853 

0,960 

0.570 

0,598 

0.628 

0.690 

0.713 

0.762 

0.811 

0.841 

0,876 

0,893 

0.940 

0,963 

0,976 

0,996 

0.978 

0.979 

0.977 

2.91 

3.82 
5.67 
6,36 

7.82 
9,27 

10.18 

11.20 

11.73 

13.09 

13,82 

14.18 
14,80 
14.25 
14.29 

14.18 

8.35 

11.45 

19.80 

24.33 

36.97 

49.86 

62.74 

75.38 

83.73 

126.67 

152.19 

177.72 

195.61 

203.48 

229.0 
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Figure  A.l  Dimensional  velocity  distribution  (experiments  S-j  and  S^) 
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Figure  A. 2 Dimensional  velocity  distribution  (experiment  S^) 
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Figure  A. 3 Dimensional  velocity  distribution  (experiments  and  U^) 


203 


Figure  A. 4. 


Dimensional  velocity  distribution  (experiment  AA-j ) 
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Figure  A. 5 Dimensional  velocity  distribution  (experiment  AA^) 
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Figure  A. 6 Dimensionless  velocity  distribution  outside  of  the  linear  zone  (experiment  S,) 
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Figure  A. 7 Dimensionless  velocity  distribution  outside  of  the  linear  zone  (experiment  S9) 
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Figure  A. 8 Dimensionless  velocity  distribution  outside  of  the  linear  zone  (experiment  S.) 
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Figure  A. 9 Dimensionless  velocity  distribution  outside  of  the  linear  zone  (experiment  U.) 


209 


Figure  A. 10  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (experiment  U ) 
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I 


Figure  A. 11  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (experiment  AA, ) 


EXPERIMENT  a A2 
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Figure  A. 12  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (experiment  AA9) 


EXPERIMENTS 
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Figure  A. 13  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (experiment  S, ) 
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Figure  A. 14  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (experiment  S?) 


EXPERIMENT  S 
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Figure  A. 15  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (experiment  Sj 
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Figure  A. 16  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (experiment  IL) 
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Figure  A. 17  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (experiment  IL) 
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Figure  A. 18  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (experiment  AA, ) 
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Figure  A. 19  Dimensionless  velocity  distribution  outside  of  the  linear  zone  (experiment  AA9) 
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